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I. Noţiuni de teoria grafurilor 


1.1. Grafuri neorientate şi grafuri orientate 


Pentru a explica ce este un graf vom da trei exemple în care 
apare această noţiune. 
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1. În figura 1.1.1 este reprezentată o parte a schemei căilor 
ferate din ţara noastră. 


Această, schemă ne indică legăturile directe dintre unele locali- 
tăți ale ţării pe calea, ferată. 


i Un asemenea desen este util atunci cînd vrem să ne deplasăm 
dintr-un oraș în altul urmînd calea cea mai scurtă, dacă cunoaștem 
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distanţa în kilometri dintre fiecare două, orașe unite printr-o linie. 


ie în circui E i i le ţării, întor- 
O excursie în circuit, care trece prin anumite orașe a i, T 
cîndu-ne în localitatea ae unde am plecat, trebuie concepută urmă- 
rindu-se pe o astfel de schemă, diferitele trasee posibile. Se spune că 


desenul din figura I.1.1 reprezintă un graf. Nodurile sau vîrfurile 


grafului sînt notate cu nume de localităţi, iar muchiile sau liniile 


graiului sînt liniile care reprezintă, legăturile directe pe calea ferată 


dintre două localităţi. - 


2. Să considerăm acum. un -exemplu din fizică, și anume calculul 
intensităţilor curenților care trec prin ramurile unei rețele electrice 
(fig. I.1.2). 


7 . tem.2 
Fig. 1.1.2 


Pentru a rezolva această problemă, cunoscînd schema, rețelei, 
tensiunile electromotoare şi valorile rezistențelor, vom serle legile 
lui Kirchhoff relative la noduri şi la ochiuri de rețea. Făcînd abstracție 
de elementele de circuit care se găsesc pe ramurile schemei putem 
desena, schema în felul următor (fig. 1.1.3): Pentru aplicarea celor 
două legi trebuie să ne fixăm un sens de parcurgere a ochiurilor 
de reţea din figura I.1.2. Obţinem astfel desenul din figura TA. 
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Fig. 1.1.3 Fig. 1.1.4 


Atât desenul din figura I.1.3 cît şi desenul din figura I.1.4 este 
un graf. Primul graf se numește graf neorientat, iar al doilea se 
numeşte graf orieniat, legăturile dintre noduri avînd o anumită 
orientare. i 


3. Desenul din figura 11.5 reprezintă o moleculă de tricvi- 
nacenă, care este un izomer al cielopolienei CH, 

Acest desen este un graf de un tip deosebit, care admite exis- 
tența unor muchii multiple între vîrturile sale, corespunzătoare ` 
valenței trei a grupării CH. 

Un graf cu muchii multiple se numeşte multigraf. 

Un graf G este definit printr-o mulțime X de vârfuri sau noduri 
și printr-o mulțime U de perechi de vîrturi din X. Mulțimea U ne 
arată ce legături există între vîrfurile grafului. Graful G cu X 
mulțimea yîrfurilor şi U mulţimea legăturilor dintre virturi se 
notează G = (X, U). De exemplu, un graf G are ca mulţime de 


` virfuri X = {2 a, g a sl, iar mulţimea U este formată din 


perechile [£n 3], [2u fs [tz al, [za Za, La £5] 
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Fig. 1.1.5 Fig. 1.1.6 


Acest graf G poate fi desenat în plan reprezentînd vîrfurile 
prin puncte, iar perechile din U prin linii (drepte sau curbe) care 
unesc vîrfurile care apar în perechile respective (fig. 1.1.6). 

Într-o astfel de reprezentare nu contează poziţia relativă, a, 
vîrfurilor grafului şi nici modul cum sînt desenate legăturile definite 
de mulțimea, U de perechi de vîrfuri. 

Astfel, toate desenele din figura 1.1.7 reprezintă unul şi acelaşi 
graf cu 14 virfuri, notate cu numerele de la 1 la 14. 

Mulțimea U pentru acest graf este formată din perechile [1, 2], 
[1 3], [ls 4], [1, 14), [2, 5], [2, 6], [2, 1], (2, 8], [3, 9], [3, 10], 
[3,11], [4, 12], 14, 13]. 
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Fig. 1.1.7 


Un astfel de graf se numește neorientai, deoarece cele două 
vîrfuri care apar într-o pereche de viriuri din U nu formează o 
pereche ordonată, nici unul nefiind privilegiat; în raport cu celălalt. 

De exemplu perechea [1, 2] poate fi citită la fel de bine [2, 1]. 
Pentru un graf neorientat o astfel de pereche se numeşte muchie, 


iar mulţimea, U se numeşte mulţimea muchiilor. O muchie se notează 


folosind paranteze drepte. 

Un triunghi poate fi privit ca un graf în care virturile sint tocmai 
virfurile triunghiului, iar muchiile sînt laturile triunghiului (fig. 1.1.8). 

Sint însă, situaţii cînd legăturile dintre virfurile unui graf sînt 
orientate de la un vîrf către alt virf, cum este graful din figura 1.1.4. 

Pentru a da un alt exemplu, să considerăm următoarea, situație : 
Academicianul A îl cunoaşte pe profesorul P (care i-a fost student), 
profesorul P îi arè ca elevi pe E, şi Fa, care sînt elevi la școli dife- 
rite şi nu se cunosc între ei. Æ, și A se cunosc reciproc, dar E, îl. 
cunoaște pe A numai de la televizor. Putem reprezenta relaţiile 
dintre A, P, E, E, printr-un graf cu 4 vixturi, ducînd o săgeată 
de la virful a către virful y, dacă w îl cunoaşte pe y. Se obține 
astfel un graf desenat în figura 1.1.9. Acest graf se numeşte graf 
orientat. El este caracterizat de mulţimea, de virfuri X = (A, P, 
E, E} şi de mulţimea, U, ale cărei elemente se numesc arce : 


U = (4, P), (P, A), (4, E), (En A), (P, E), (Bu P) (P, E») 
(En P), (Bu A). 


Pentru un are (z, y), virful œ se numeşte ewtremitate inițială, 
iar y se numeşte ewtremitate finală. Un arc va fi notat prin cele 
două extremități ale sale în paranteze rotunde. 

Arcul (z, y) se reprezintă grafic ca în figura I.1.10, ducînd o 
săgeată orientată de la æ la y. O astfel de pereche (w, y) se numește 
pereche ordonată (în care ordinea vîrturilor are importanță). 

Deci în cazul unui graf neorientat legăturile dintre vîrfuri se 
numese muchii, iar în cazul unui graf orientat legăturile dintre 
virfuri se numesc arce, care sînt perechi ordonate de viriuri. Să 
observăm că arcele sînt elemente ale produsului cartezian X x X, 
deoarece sînt perechi ordonate (%, y) cu æ, y €X, mulțimea virfu- 
rilor grafului. Atît în cazul unui graf neorientat, cât şi în cazul unui 
graf orientat, legătura dintre un virf şi el însuşi se numeşte buclă 
(fig. 11.11). 


Fig. 1.1.8 
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Fig. 1.1.9 
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Fig. 1.1.10 
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Pentru a înțelege mai bine semnificația 
unei bucle, să construim un graf în felul 


următor : 
Pornim de la numerele 2, 3, 4, 6, 12, pe ` 


care le reprezentăm ca vîrturile unui graf, du- 
cînd o săgeată de la virful la virful j dacă 
i este un divizor al lui j. Se obţine astfel 
graful din figura 11.12. . 

Pentru acest graf există bucle în fiecare 
virf, deoarece orice număr este un divizor 
al lui însuși. i 

Un subgraf al unui grat G = (X, U), se 
defineşte în mod asemănător cu noţiunea; de 
submulțime : el este un graf G, = (£n U) 
unde Xc X, iar U, reprezintă toate legătu- 
rile dintre vîrturile din X, care existau în 
graful G. Menţionăm că 'am notat prin e 
inclusiunea, nestrictă. De exemplu, pentru 
graful orientat din figura 1.1.12, subgraful 
cu mulțimea, de vîrturi (2, 3, 6) este desenat 
în figura I.1.13, iar subgraful cu mulțimea; 
de virfuri (2, 4, 12) este desenat în figura 1.1.14. 
Pentru graful orientat din figura 1.1.9 sub- 
graful care are ca mulțime de virfuri mulți- 
mea (A, P, E.) este desenat în figura, 1.1.15, 
iar subgraful care are ca mulţime de vîrturi 
mulţimea (P, E, E) este desenat în figura 
I.1.16. x 

Pentru graful neorientat din figura 11.17 
subgraful cu mulțimea de vîrfuri (2, 3, 4, 5} 
este desenat în figura I.1.18, iar subgraful cu 
mulțimea de vîrfuri (1, 2, 3, 5} este desenat 
în figura I.1.19. 

Un graf parțial al unui graf G = (X, U), 
este graful G, = (X, U.) unde V,c U, deci 
un graf parțial al unui graf G se obține din 
G prin suprimarea anumitor legături dintre 
vîrfurile lui G. De exemplu, graiul din figura 


1.1.20 este un graf parţial al grafului din 


figura, 11.17, iar graful din figura I.1.21, 
este un graf parţial al grafului orientat din 
figura I.1.12. 


Detiniţie: Într-un graf neorientat, gradul unui vîrf æ, notat die), 
este numărul muchiilor care se leagă de virful w. 


De exemplu, pentru graful din figura I.1.17 obţinem d) = 2, 
3, d(3) = 3, d(4)= 3, (5) 
Prin definiţie, un vîrf izolat este un virf care are gradul zero, 


a2) 


adică, nu se leagă prin muchii de nici un alt virf al grafului. 


A 
A 
A 
2 
PS E 
Fig. 11.12 Fig. 1.1.13 Fig. 1.1.44 Fig. 1.1.15 
7 
2 5 
2 5 
E. 
< SL >E 
5 P 3 4 3 + 
Fig. 1.1.16 Fig! 1.1.17 
i 
š : 
2 { N 5 
l a 
` 3 3 4 
Fig. 1.1.19 Fig. 1.1.20 Fig. 1.1.21 


Teorema 1.1.1. 


Pentru un graf neorientat G = (X, Uj cu m muchii, suma gra- 
delor vîrfurilor sale este egală cu pă i i : i 


De exemplu, graful din figura 1.1.17 are 7 muchii și &(1 
-+ (2) + d(8) + da) + d5) = 2 4+3 +3 +343 = w RE 
strația teoremei se face ţinind seama că, în. suma gradelor xîrfurilor 
graiului fiecare muchie [a, y] este numărată exact de 2 ori, ea 
contribuind cu o unitate la d(x) şi cu o unitate la (y). 


Aplicaţie. 
Să arătăm că un poligon convex cu n virfuri are E i mad 
a ci 2 
diagonale. i i 

„__ Un poligon convex cu n virfuri, pentru care s-au dus toate 
diagonalele, reprezintă un graf neorientat cu n virfuri care are 


proprietatea că oricare două virfuri sînt unite printr- i 
a LI aa printr-o: muchie 


Ke 
Fig. 1.1.22 


Un astfel de graf se numeşte graf complet cu n virfuri şi se 
notează K,. % 

__ Pentru a calcula numărul muchiilor acestui graf, să observăm 
că gradul fiecărui virf z este egal cu n — 1, fiecare vîrf fiind legat 
prin muchii de toate celelalte n — 1 vîrfuri ale grafului. Deci suma, 
gradelor vîrfurilor grafului este egală cu n(n — 1), care cu teorema 
anterioară este egală cu 2m, dublul numărului de muchii. 

f n(n — 1) i y 
Deci m = ———— = C} muchii. Pentru a afla numărul 
diagonalelor, trebuie să scădem din acest număr numărul laturilor, 
care este n, deci numărul diagonalelor unui poligon convex cu 2 
n(n — 1) __nn—3) 
2 


virfuri este egal cu 
2 
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Grafurile complete au aplicații în studiul grupurilor sociale. 
Astfel, dacă reprezentăm membrii unei colectivități umane prin 
Virfurile unui graf, legătura dintre două, vîrfuri reprezentind o 
relaţie de prietenie, de simpatie sau faptul că, doi membri ai colecti- 
vităţii au interese saú preferințe comune, subgraturile complete 
ale graiului care reprezintă grupul social studiat au semnificația 
unor subgrupuri sociale. Într-un astfel de subgrup oricare doi membri 
sînt; prieteni sau au preferințe şi interese comune. 

Astfel, să considerăm trei persoane A, B, C, pe care le repre- 
zentăr ca virturile unui graf şi să ducem o muchie între două vârfuri 


` ale gratului numai dacă persoanele respective sînt prietene. 


Se pot; obţine diferite situaţii, unele dintre ele fiind reprezentate 
în figura I.1.23. 


A 4 A „A 

A i i o o 

& Ci Bă e 8 e 8 c 
a) 3 é} c) d: 


Fig. 1.1.23 


Astfel în figura I.1.23a) graful este complet, în figura, 11.23 e) graul 
se descompune în două subgrafuri complete (A, B} şi (0), iar în 
figura, I.1.234) cele trei virfuri A, B, O formează fiecare cîte un 
subgraf complet, fiind virfuri izolate. În figura I.1.23 b) putem forma 
cu virfurile A, B un subgrat complet şi cu virful C un alt subgraf 
complet; sau cu vârfurile A şi C un subgraf complet şi cu B un alt 
subgraf complet. Observăm că în acest caz mai există, câte o muchie 
fie (4, C], fie (A, B], care uneşte virturi din subgrafuri complete 
diferite. În sociologie se spune că structura, reprezentată, de graful 


` din figura, 1.1.23b) nu este echilibrată. Într-adevăr, A se găsește 


într-o situaţie delicată, el fiind prieten cu B şi C, care nu sînt 
prieteni între ei. 

Studiul echilibrului unui graf care reprezintă un grup social 
se poate face şi în cazul unui număr mai mare de virturi, studiind 


subgrafurile sale complete. 


+ Un termen sinonim este acela de grafe. 
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În cazul grafurilor orient zü 
i rail ate, un graf complet graf 
pei a că orbire două, virfuri ale hd g e prea re nd 
au printr-un ar ; i È și 
7 PA Ta TC (7, 2) sau prin ambele arce (s, y) şi 


Fig. 1.1.24 


i Astfel, graful orientat din fig 
piara viris fa ăi T -ain unite printr-un arc, dar subgraful 
u mu rfuri P} sau i 
iu A, E B} sînt fiecare el alu. ic n en it ii d 
ormulele de structură ale substanțel imi 

de si i J or chim î i 
pentru care. legăturile dintre vîrfuri sînt tocmai este al ee 
rilor sau atomilor care compun molecula. e i 


f Fast apa are molecula, reprezentată în figura 1.1.25a, acetilena 
ue Stiri sr ae via a gaus I.1.26a, molecula de benzen 
e igura I.1.27a, cea d ă i & 

iar cea de eloropren în figura, Tia A i 


CH 
TAS cu cH 
HC == CH CH > 
Sy Ă CH / 


CH 


Fig. 1.1.25a. Fig. 1.1.26a. Fig. 1.1.27a 


cho 

be 
CH—on 
“CH—OH 
CH—0H CH, = CH — C = Ch, 
Chi OH |. 


ci 


Fig. Î.1.28a. Fig. 1.1.29a 
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ura 1.1.9 nu este graf complet,- s. () $ E 


„Fig. 11.25. Fig. I-1.26b. Fig. 1.1.27b. Fig. L4.28b. Fig. 1.1.29p. 


Acestea sînt grafuri pentru care vîrturile sînt atomii (respectiv 
grupările) din moleculă, legăturile dintre ele fiind muchiile grafutui 


(fig. 1.1.25b — I.1.29b). 


Graturile din figurile : I.1.26b, I.1.27b, 1.1.29b, nu sînt grafuri 
obişnuite în sensul definiţiei date, deoarece între anumite perechi 
de vîrfuri există mai multe muchii. 

Un astfel de graf cu muchii multiple se numeşte multigraj. 
Într-un graf (sau multigraf) care este desenul moleculei unei substanţe, 
gradul unui virf este tocmai valența atomului (grupării) respective. . 


Scurt istoric. 


Teoria grafurilor este o ramură a matematicii, relativ 
tînără : prima lucrare de teoria grafurilor a fost serisă de 
Euler acum două secole. În secolul trecut multe rezultate 
în teoria grafurilor aw fost date de Cayley, în timp ce 
Kirchhoff a aplicat graiurile la studiul circuitelor electrice. 

Studiul teoriei grafurilor a fost stimulat de apariția 
în 1936, la Leipzig, a cărții lui D. König de teoria grafu- 
rilor orientate şi neorientate, prima de acest gen din 
lume. Interesul pentru teoria grafurilor a crescut mult în 
ultimele decenii, odată cu apariția economiei matematice. 
În vreme ce la început teoria grafurilor părea destinată 
să dea soluții doar unor jocuri și amuzamente matema- 
tice, ea se aplică astăzi la elaborarea deciziilor optime în 
multe probleme economice, în studiul izomerilor substan- 
telor, în sociologie și psihologic, în lingvistica matematică 
şi inginerie. Teoria grafurilor poate fi considerată ca o 
ramură a teoriei mulțimilor, însă ea are de asemenea 
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multe interferenţe cu geometria (topologia), cu algebra, 
şi calculul probabilităților. În ultimul timp, teoria grafu- 
rilor se consideră ca făcînd parte din domeniul mai larg 
al matematicii combinaitorii (combinatoricii). Printre mate- 
maţicienii care au adus contribuţii importante în dome- 
» niul teoriei grafurilor în ultimele decenii, cițăm pe 0. Ore, 
-0. Berge, PF. Harary, P. Brdös, W. Tutte. ` 
În această ramură nouă a matematicii şi-au adus 
contribuția şi matematicienii români, autori ai unor rezul- 
tate care s-au impus deja pe plan mondial în acest domeniu, 
cît şi ai unor monografii publicate în limba română, 
care tratează despre teoria grafurilor și aplicațiile acesteia, 
De remarcat că în ţara noastră teoria grafurilor a 
primit un aport însemnat şi din-pariea celor care o aplică 
în economie, chimie, inginerie. 


[.2. Drumuri şi circuite în grafuri : 


O noţiune importantă în teoria grafurilor este aceea de arin. 
Pentru un graf G = (X, U) un drum de la virful v la vîrful y este o 
succesiune de vîrfuri pe care o notăm (2, Lı, a: - :5%p, Y) În cazul 
grafului orientat și [®, d, 8a. - -2p Y] în. cazul grafului neorientat, 
astfel încât : 

— dacă graful G este neorientat, perechile [e, mh. [t a]l <- 
wp y] sînt muchii ale grafului ; | 

— dacă graful G este orientat, perechile ordonate (2, 2), (n 22), 
(8p Y) sînt arce ale gratului. 

Dacă toate viriurile w, 2 Sa.. .;8p y Sint distincte, drumul 
(&, Lu- -8p Y} se numeşte drum elementar, iar dacă s = y, adică, 
drumul se întoarce în punctul de plecare, acest drum, se numește 
circuit. : : 


În figura 1.2.1 vom folosi următoarele notații: 


B = Bucureşti C = Cluj 

P = Ploieşti O = Orşova 
BV = Braşov CV = Craiova 
S = Sibiu BZ = Buzău 
T = Timişoara BC = Bacău 
A = Arad SV = Suceava 


D = Deva I = aşi 
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aful din figura 1.2.1. reprezintă o parte a reţelei C.F.R. din 
tara Pee. Deoarece legăturile pe calea ferată se fac i ambao 
sensuri; am presupus că acest graf este neorientat (în cazul ciro 
tației rutiere într-un oraş pot exista străzi cu sens unic, deci grafu 
trebuie considerat orientat). . i 


SV 


Fig. 12.1 


i: ; (BY, 8, T]; 
Exemple de drumuri: [B, P, BZ, BO, SV]; LBV, Ñ, i 
i0, D, AT Toate aceste drumuri sînt elementare. Un drum neele- 
mentar este de exemplu [BV, 8, D, A, T, 8, OV] pentru care se 
epetă, nodul S. 
SE Exemple de circuite : [B, P, BY, S, CV, B]; (8, Z, 0, OV, 8]; 
LB, P, BZ, BO, SV, 0, D, A, T, O, CV, B]; [BZ, I, BZ]. i 4 
Un drum (circuit) care trece o singură dată prin toate vârfurile unui 
graf se numeşte drum (eirouit) hamiltonian. Deci un drum hamiltonian 
este un drum elementar care trece prin toate virfurile unui graf. 
Pentru graful din figura I.2.1 nu există un circuit hamiltonian, 
deoarece-un astfel de circuit venind şi plecînd din virful 1 ar vali 
să treacă de două ori prin virful BZ, ceea ce contrazice definiţia 
unui circuit hamiltonian. Dacă suprimăm însă virful I, împreună 
cu muchia [I, BZ], subgraful obținut conține ùn circuit hamiltonian, 
i anume : an 
ip P, BZ, BO, SV, 0, BV, 8, D, 4, T, 0, OV, B), circuit care 
pliere de d hamiltonian își are istoria 
`- Denumirea, de circuit sau de drum hamiltonia i are istori 
sa. Astfel, în anul 1859 matematicianul Hamilton a pus în circulație 
eor dr lat făcut din lemn 
Partea sa principală era un dodecaedru regulat tâcu 
(fig. 1.2.2). Aces a un poliedru cu 12 fețe care sint toate penta 
goane regulate, iar în fiecare din cele 20 de virfuri se întilnese cîte 
3 muchii. - 


15. 


Fiecare vif al dodecaedrului lui Hamilton era marcat cu numele 


unui oraș important: Bruxelles, Canton, Delhi, Frankfurt ete. 


Jocul consta în a găsi un drum de-a, lungul muchiilor dodecaedrului 
care să treacă prin fiecare oraș exact o dată şi să se întoarcă, în 


- Fig. Ts 


oraşul din care a plecat. Ordinea trecerii prin primele oraşe era 
stabilită de la bun început. i 

Pentru a ușura memorarea, treceriloy» efectuate, în fiecare virf 
al dodecaedrului era cîte un cui cu o floare mare, astfel încît în 
jurul acestor cuie putea să se întindă, un fir care să indice drumul 
parcurs în această călătorie imaginară în jurul lumii, 

Un circuit hamiltonian în reprezentarea, plană a grafului dode- 
caedrului este trasat cu linii mai groase în figura, 1.2.3. 


Fig. 1.2.3 


Graiurile orientate și complete au o proprietate importanță, 
relativă la existenţa unui drum hamiltonian, dată de teorema urmă- 
toare : 


Teorema [.2.1. Orice graf orientat şi complet are un drum hamiltonian. 
Demonstrația acestei proprietăți se face prin inducţie după 
numărul virfurilor grafului, n.. Pentru n = 2 graful are două virturi 


æ şi y. Deoarece graful este complet, conform definiţiei, există un. 
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arc între aceste două viriuri, de exemplu (v, y), care este tocmai 
drumul hamiltonian căutat. Presupunind proprietatea adevărată 
pentru n — 1 vîrfuri, s-o demonstrăm pentru n virfuri. 

Dacă graful G are virfurile z,, So.. pt, Subgraful format cu 
vîrfurile a. .. 8n; este complet (conform definiţiei), deci el conține 


-un drum hamiltonian (conform ipotezei inducției). Fie acest drum 


(Lis Loro- -5a a). 


Xa 


Fig. 1.2.4 


Deoarece graful este. complet, există, cel puţin un arc între m, 
Şi 2. Dacă acest arc este (Sn, 2), atunci am obţinut un drum hamil- 
tonian cu n vârfuri, şi anume : (Lay Lys - - -3%a-a). Să presupunem con- 
trariul şi anume că există numai arcul (2, Va) 

Repetind acelaşi raționament pentru perechea de virfuri z, 
Și Zau dacă există arcul {2 y Sa) atunci putem forma drumul 
hamiltonian (2, o-- -3u-ay La). În caz contrar există numai arcul 
(Las n-a). Să presupunem deci că sintem în cazul cel mai nefavo- 
rabil, adică există numai arcul (e, 2) Şi arcul (£a aa) (fig. 
1.2.4). Deoarece z, se leagă prin arce cu vîrfurile &,,.. noa VOr 
exista pintre acestea două vîrfuri vecine, Sp şi Zu pentru care 
există arcele (Sr, Vn) Şi (La, Dia). Într-adeyăr, din a pleacă un 
arc ÎN Sp dar în Vp ultimul virf al drumului, vine un are din 
La, deci la un moment dat, plecînd din æ, către £s., sensul arcelor 
către z, trebuie să se schimbe. Dar şi în acest caz am găsit un 
drum hamiltonian cu n virfuri şi anume: (8, aye- cy m Dei 
---38a1) (fig. 1.2.4) şi teorema, este demonstrată. Un tip special de 
grai orientat şi complet este așa numitul graf-turneu, care are pro- 
prietatea, că între oricare două virturi ale sale există un are și numai 
anul. Denumirea, de graf-turneu provine din aceea că dacă, reprezen- 
tăm prin vârfuri jucătorii care iau parte la un turneu, astfel încât 
fiecare jucător să joace cu toţi ceilalţi (turneu neeliminator) şi să 
nu există meciuri nule, ducînăd un arc de la virful care reprezintă, 
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jucătorul i la vîrful care reprezintă jucătorul j dacă î l-a învins 

în meci direct pe j, graful astfel obținut este un graf-turneu (fig. 1.2.5). 
Din teorema, 1.2.1, rezultă că un astfel de graf admite un drum 

hamiltonian, de exemplu: (2, 5, 3, 4, 1), adică jucătorii au fost - 


1 


Fig. 1.2.5 


aranjați astfel încît 2 l-a învins pe 5, 5 l-a învins pe 3, 3 l-a întrecut 
pe 4 şi 4 pe 1. Deoarece jucătorul 1 a înregistrat numai înfrin- 


. geri, toate arcele grafului complet care ajung la virful 1 îl au pe 1 


ca extremitate terminală, deci 1 este nodul terminal al drumului 
hamiltonian, drum care poate fi găsit prin recurenţă, aşa cum se 
face în demonstrația teoremei 1.2.1. 

Jucătorii 3 şi -5 au cele mai multe victorii (cite 3), însă 5 poate 
fi considerat cîştigător al turneului, deoarece are o victorie directă 
asupra lui 3. ` A 

O aplicație a determinării drumurilor hamiltoniene într-un graf o 
constituie succesiunea optimă de faze în procesul prelucrării unor 


-piese la o maşină. Să presupunem că trei piese, notate Pi, Pa, Pa 


trebuie prelucrate pe o aceeaşi maşină. La trecerea de la prelucrarea 
unei piese la prelucrarea altei piese se pierde un anumit timp, deoarece 
la maşină treþuie demontate anumite dispozitive şi montate alte 
dispozitive, specifice prelucrării piesei respective. Acest timp de 
trecere poate fi interpretat şi ca un timp de pregătire a mașinii 
pentru prelucrarea altei piese (curățire, ungere ete.), sau pentru 
spălare şi schimbarea de coloranți la imprimarea unor ţesături. 

'Timpii de aşteptare între prelucrarea a două, piese sînt asociați 
arcelor grafului complet din figura 1.2.6. Astfel, vîrfurile grafului 
reprezintă, piesele, iar arcul (P, Pp) cu î,j=1, 2; 3 gi t #j are 


asociat un timp egal cu: timpul de aşteptare între prelucrarea lui 
P, şi prelucrarea lui P; pe maşină. Dacă definim „timpul“ asociat ` 


unui drum ca fiind suma timpurilor asociate arcelor care compun 
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acest drum, atunci drumul hamiltonian de timp minim va reprezenta 
toemai succesiunea, de prelucrare a celor trei piese pe mașină pentru 
care timpul total de aşteptare între prelucrări este minim (se negli- 
jează timpul de pregătire a maşinii pentru prelucrarea, primei piese), 


Fig. 1.2.6 


Suecesiunile posibile de piese pe maşină sint date în continuare, 
fiecare reprezentind cîte un drum hamiltonian în graful din figura, 
I.2.6: - 


Timp de așteptare 


(Po Pa Ps) 34+383=6 
(Po Pa Pa) - 0+6=6 
(Pa Pe Ps) ' 2+0=2 
(Pa Pa P) 3+l=4 
(Ps P; P,) 1+3 =4 
(Pa Pa P) 6+2=8 


Drumul hamiltonian de timp minim este (Pa, Pu P3) şi el cores- 
punde succesiunii optime de prelucrare P}, Pu Ps. Drumul hamil- 
tonian de timp maxim este (Py, Pa, P4) și el corespunde succesiunii 
celei mai nefavorabile de prelucrare şi anume Ps, Pa Pi. 

Dacă ţinem seama, și de timpul de pregătire a maşinii pentru 
prelucrarea primei piese, problema se reduce la problema găsirii 
drumului de timp minim din graful desenat în fig. 1.2.7, care pleacă 
din Pi, Pz sau P; şi trece o singură dată prin Pi, Pa Pa 

Am presupus că timpul necesar pentru pregătirea mașinii la, 
prelucrarea, primei piese P, este de 2 unități, timp asociat arcului 
(Pi, Pı), timpul necesar pentru pregătirea mașinii la prelucrarea, 


piesei P, este de 4 unităţi, timp asociat arcului (P2, P), iar timpul 
necesar pentru pregătirea mașinii la prelucrarea, piesei P, este de o 
unitate, timp asociat arcului (Ps, Pg). = 


Drumurile care pleacă din Pi, P; sau Ps şi care trec o singură - 


dată prin Py, Pa, P, în acest graf sînt următoarele : 


Timp de aşteptare 


(Pi Pı, Pa Pa) 64+2=8 
(Pi, Pu Pa Pa) A 6+2=8 
(Po Pa Py P;) 2+4=6 
(Pa Pa Pa Pi)’ 2 4+4=8 
(Ps, Pa Pr Pa) 4+i=5 
(Ps, Pa, Pa Pi) 8+1i=9 


Drumul de timp minim tiind (Ps, Ps, Pis Pa); rezultă că succe- 
siunea optimă de piese pe maşină s-a schimbat, devenind P}, Py 
P, în timp ce succesiunea cea mai nefavorabilă a rămas aceeaşi. 

Numărul drumurilor hamiltoniene ale grafului neorientat şi 


! 
complet cu w virfuri, graf notat prin E, este egal cu. În- 


tx-adevăr, primul vîrf al drumului poate fi ales în n moduri, al doilea 
în m — 1 moduri,..., ultimul într-un singur mod, rezultind în 
final n! posibilităţi de alegere. Însă, fiecare drum poate fi pareurs 
în două sensuri diferite, rezultind două permutări distincte, deci 
numărul n! al permutărilor. celor n vîrfuri trebuie împărțit la doi. 


oare o az 


Să observăm că acest număr creşte foarte repede odată, cu 4. În 


timp ce pentru graful K; din figura I.2.8 ave =60 drumuri ha- 


3 ! 
miltoniene, pentru graful K,, avem m 2 5-10% astfel de drumuri! 


Fig. 1.2.8 


Numărul circuitelor hamiltoniene ale grafului complet cu 
i —i)! 
n virfuri K, este mai'mic, fiind egal cu a= pentru n >3 (pro- 


blema 6). 7 

Dacă fixăm extremitățile drumurilor, de exemplu «, şi 4 pentru 
graful din figura 1.2.8, atunci există (n — 2)! drumuri hamiltoniene 
cu extremitățile fixate. Într-adevăr, fiind extremităţile drumului, 
celelalte n — 2 virfuri formează o permutare cu n — 2 elemente, 
numărul acestora fiind egal cu (n — 2)!. De exemplu, pentru graful 
din figura 1.2.8, drumurile hamiltoniene care au extremităţile mi 
şi æ; Sînt în număr de (5 — 2)! = 3! = 6 şi anume: 


[Zu Bo Pa Pa Bel 
[zu Lar Bos Za Loly 
[las Ly Pa Lay 4) ȘI 


[Za a Vas Bar 25], 
[i Pa a a Vsh 
[Ei Das a Boy Vgl- 


În teoria grafurilor există numeroase rezultate consacrate unor 
condiții suficiente de existență a drumurilor sau circuitelor hamil- 
toniene pentru grafurile orientate sau neorientate. Cea mai mare 
parte a acestora se reduc la condiţii impuse gradelor vîrfurilor gra- 
fului, căci este clar că dacă un graf G conţine suficient de multe 
muchii, el poate conține şi un drum sau un circuit hamiltonian. O 
condiție simplă este dată de problema 30. Noi nu ne vom ocupa 
însă aici cu probleme de existență, ci cu probleme de optimizare. 
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adică lungimea drumului hamiltonian (7, îs. - 


„Aceste probleme se referă la graful complet, care nu numai că 
conține un drum sau un circuit hamiltonian, ci conține, aşa cum 
am văzut, prea multe şi problema alegerii unuia, optimal poate îi 
dificilă. 


Un exemplu de optimizare a fost acela al prelucrării unor piese . . 


la o maşină. În cazul general, această problemă se enunță astfel : 
Pe o maşină M trebuie prelucrate n piese: P, Pa.. 3 Pa 


La, trecerea de la prelucrarea, piesei. P, la prelucrarea piesei È; se ` 


consumă un anumit timp, notat cu 6, care poate fi un număr real 
pozitiv sau zero. 


Am obținut astfel un graf orientat și complet G = (X, U) ale 


cărui virturi reprezintă cele n piese, avînd asociată o matrice pătrată, 


O = (Ciuma... de ordinul n, pe care o vom numi matrice a 
costurilor sau a "lungimilor arcelor graiului. 

Presupunînd că această matrice este simetrică, adică, ĉu = i SA 
pentru orice î,j = 1,.. :,n (vom defini e, = 0 pentru i = 1, 2, . . n); 
graful G poate fi considerat neorientat. 

n(n — 1) 
2 i 
muchii, fiecare muchie [5, 3] avînd asociat costul (lungimea) e. 

Se pune problema găsirii acelei succesiuni de prelucrare a celor 
n piese Pi, Pa...,P, la maşina M, astfel încît timpul total de 
nefuncţionare al maşinii M să fie minim.. 

Presupunînd că această succesiune ar fi : [Pa; Paseo 
suma timpilor de pregătire. a maşinii este- tocmai : 


El este un graf complet, deci mulțimea U conţine 02 = 


Giria F Gigi, H 00 F Cinane 


astfel definit. Vedem deci că modelul matematice al acestei probleme 
de minimizare a timpului de inactivitate -a utilajului M este o pro- 
blemă de drum hamilţonian minim în graful complet cu matricea 
costurilor 0 asociată acestei probleme. Dacă matricea C nu este 
simetrică, atunci graful orientat G are proprietatea că între două, 
vîrfuri oarecare există două arce de sensuri opuse, care în general 
am. costuri diferite, adică, Ce; É Cu. acest caz trebuie căntat un 
drum hamiltonian minim într-un graf orientat. 

Să considerării un alt exemplu de problemă înrudită. Fie un 
oficiu poştal care dispune de o maşină, pentru a colecta corespon: 
dența de la n cutii poştale aflate în raza sa. 
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Pinas Pal, , 


-sîa] în graful complet . 


AT 


inau umagen aen me aaa TIR sa 


Se cunoaşte matricea distanțelor minime D = (a)r ja1, =, n, unde 
dy reprezintă distanța minimă (în km) de lẹ cutia cu numărul 
i la:cutia cu numărul j. Dacă în oraş există străzi cu sens unic, 
se poate ca matricea D- să nu fie simetrică, deoarece pot exista, doi 
indiei ¿ şi j astfel ca dus du. 

Se pune problema, găsirii traseului optim al mașinii care trebuie 
să colecteze corespondența de la toate cele n cutii poştale, întor- 
cîndu-se la oficiu. Această problemă este o problemă de circuit 
hamiltonian minim într-un grai neorientat (dacă presupunem matricea 
D simetrică) sau orientat (cînd matricea D nu este simetrică). 

Această problemă mai este cunoscută şi sub denumirea de 
„problema voiajorului comercial“, voiajor care trebuie să viziteze 
n oraşe cu mostre din producţia fabricii pe care o reprezintă, întor- 
cîndu-se de unde a plecat, astfel încît distanța parcursă să fie minimă 
(ceea ce corespunde unui cost minim al deplasării). O altă variantă 
posibilă, a problemei cutiilor poștale este următoarea : Maşina care 
colectează corespondența nu pleacă în cursă de la oficiul poştal, ci 
de la' un garaj, situat undeva în oraș. Ea trebuie să se întoarcă 
însă a oficiu pentru a preda corespondenţa, strinsă de la cutiile 
poştale. . 

În acest caz avem de-a-face cu o problemă de drum hamiltonian 
minim cu extremităţile fixate, respectiv garajul şi oficiul poștal. 

Dacă corespondenţa, se ridică de mai multe ori pe zi, a doua 
cursă va îi un circuit hamiltonian, mașina plecînd de la oficiu și: 
înapoindu-se tot la oficiu. i 

La fel se pot pune problemele de distribuire a unui produs la 
anumiți consumatori cu o maşină care execută un circuit. hamil- 
tonian. În acest caz apare şi condiţia, suplimentară ca să nu se depă- 
şească capacitatea de transport a maşinii. 

Dacă această capacitate nu este suficientă, se vor afecta mai 
multe mașini, fiecare urmînd un traseu de deplasare care este un 
circuit hamiltonian minim în subgraful format din centrele pe care 
le deserveşte maşina respectivă. 

Din analiza acestor probleme practice distingem trei tipuri de 
probleme de trasee hamiltoniene minime în grafuri : 

P1) drum hamiltonian minim ; 

P2) drum hamiltonian minim "cu extremităţile fixate, î, şi ia; 

P3) circuit hamiltonian minim. 

Pentru prescurtare, vom nota un drum hamiltonian minim 
prin DHM. Dacă mulțimea virfurilor grafului complet este X, iar 
matricea costurilor este 0, vom mai nota prin DHM (X, C) ur 
drum hamiltonian minim relativ la graful cu mulțimea virturilor 
X şi matricea, costurilor 0. 
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Pentru problema P2 vom nota prin DH Mis, sau DEM, (X, 0) 
un drum hamiltonian minim cu extremitățile fixate, å şi å în 
graful complet cu mulțimea vîrturilor X şi matricea costurilor C. 
Pentru problema P3 vom nota prin CHM sau CHM(X, C) un 
circuit hamiltonian minim în graful complet cù mulțimea, vîirfurilor 
X şi matricea costurilor C. 

Vom arăta că aceste trei probleme sînt echivalente între ele, 
adică un algoritm pentru rezolvarea uneia oarecare dintre ele este- 
în acelaşi timp un algoritm pentru rezolvarea celorlalte două, cu 
foarte mici modificări ale mulţimii X.a vîrturilor gratului sau ale 
matricei costurilor ©. . - 

Să presupunem că ştim să, rezolvăm problema, P1, adică știm 
să găsim un DHM. Să vedem cum putem rezolva problema P2. 

Vom modifica datele pentru problema P2, lăsînd neschimbată 
mulţimea X, însă schimbînd matricea costurilor C astfel : 

Adunăm la elementele liniilor i, și î2, ca şi la elementele coloa- 
nelor å, şi în o 'constantă K, mai mare ca, lungimea oricărui drum 
hamiltonian din graful G = (X, U) cu matricea costurilor C. De 
exemplu, K se poate alege suma, elementelor matricei 0. 

Elementele de pe diagonala principală a matricei O vor rămîne 
însă neschimbate, egale cu zero, deoarece ele nu au nici o semni- 
ficaţie pentru problema considerată, 


Să observăm că prin această transformare elementele situate 


pe poziţiile (i, ia) și (t, î,) se vor mări cu 2K. Deci matricea costurilor 
C va deveni: 


i PRI On HE ceea, tr Eee. Czyn 


Gna +E. 5 -à REFET ET Cii + 2K- dra + K 
ra= a : ; f 

Gia PE: e tipa IE e Deea... Cin HE 

Engst Cna tK. * Cai +E. 0 


Toate elementele care nu se găsesc pe liniile į sau fî sau pe 
coloanele î, sau i, rămîn neschimbate. 

Vom arăta acum că un DHM pentru graful 6 = (X, U) cu 
matricea, costurilor modificată, C.i este un DH Mai, pentru graful 
G cu matricea costurilor O. 

Trebuie deci să demonstrăm că orice DHM în graful G cu 
matricea costurilor modificată în modul indicat este un drum hamil- 


fi 
k 
E 
h 
d 


Ei ra roi ta 7 Deer aan Ca 


tonian minim cu extremităţile fixate, î, şi î. în graful @ cu matricea 


costurilor C. , X 
Pot apărea 3 situații : 


1) DEM (X, Caa) are ca extremități pe în şi iaj 


2) DHM (X, Can) are una dintre extremități pe i sau i, şi 
cealaltă extremitate diferită de î, şi îs; 


3) DHM (X, Cia) nu are nici o extremitate din mulțimea {iy ia}. 


Vom arăta că situațiile 2) şi 3) nu sînt posibile, ele contrazi- 
cînd ipoteza de minimalitate a drumului hamiltonian, din modul de 
alegere a constantei K. 

în cazul 1), dacă notăm lungimea unui DHM (X, Cii) prin A 
acelaşi drum hamiltonian 'va avea în graful inițial cu matricea 
costurilor C o. lungime }, astfel încît : 


l=} +2E 


din definiția matricei Cis, (fig. 1.2.9). 


Fig. 1.2.9 
- Într-adevăr, acest drum conţine numai două, muchii care se 
leagă de î, şi ip deci a căror lungime a crescut cu câte K unităţi. 
În cazul 2) avem situația : 
l=LĠ + 83K, 


deoarece singurele cazuri posibile sînt cele din figurile [.2.10 şi 
1.2.11, care corespund cazurilor cînd î, şi î2 sînt virfuri vecine 


B 


. drumul hamiltonian DHM (X, Ci) Sau. nu. 


Fig. 1.2.10 


În cazul 3) avem: 
i l=L+4K 


aşa cum se constată din, figura I.2.12 şi ime 

um se co t . fi .2.12 şi I.2.13, deoar i 

PAER REA, din matricea Ca crescut cu. DE iar ei 
a avînd o singură extremitate 4, sau i, ai crescut cu K 


D 


n N 


Fig. 1.2.12 


Să arătăm că numai cazul 1) e posibi 

2 m c ste posibil. 
Într-adevăr, presupunind prin absurd că DAM (Z, 0, ) 
Li TEZI 


are lungimea ! şi sîntem în i i lui 
a m pen cazul 2), atunci lungimea, lui pentru matri- 


l= l + 3K. 


` 


O ET 


fa N 


Fig. 1.3.13- 


Să, considerăm acum un drum hamiltonian oarecare, care are 


ca extremităţi pe i, și i de lungime L Cesare 
Atunci am văzut că CĂ gime pentru ;matricea Oiri, 
egolitaten:, ăzut că lungimea sa L, pentru matricea 0 verifică 


L= +2K 


ii 
i 
ji 


i Ne liane aa aut A 


De aici se vede că L <1, sau: 
In +2K <4 +38 K, adică: 
L—k<kK, 
inegalitate adevărată, deoarece K a fost ales mai mare ca lungimea 
oricărui drum hamiltonian în graful G cu matricea 0, deci E > L, > 


> Ii.— h, lungimile muchiilor fiind numere reale pozitive sau zero. 
- Am ajuns deci la o contradicţie, şi anume că DHM (X, Ciia) 


`- nu este minim, deoarece am găsit un alt drum hamiltonian, cu 


extremități 4 şi î2 de lungime strict mai mică. 
“La, fel se tratează. şi cazul 3), deci am demonstrat că pentru 
graful G cu matricea lungimilor (costurilor) modificată, Ci,s drumul 


- hamiltonian minim are în mod obligatoriu extremitățile 1, şi îe 


fixate anterior. 

Să arătăm acum .că dacă ştim să rezolvăm problema P2, atunci 
ştim să rezolvăm şi problema Pi. E af 

într-adevăr, să adăugăm la graful Œ două noi viriuri, Baa Și 
Bayo, deci Xa = X U {Lari Bus) Și să definim lungimea muchiilor 
[usa d] cu î <n ca fiind egală cu a, jar lungimea, muchiilor 
[tn Sa] cu j < n fiind egală cu b, cu a, b>0 (fig. 12.14). 

- Muchiei [ua p+a] îi vom da de asemenea o lungime arbi- 
trară e > 0. 

Dacă determinăm în graful cu mulţimea, de virfuri X, un drum 
hamiltonian minim cu extremităţile Za Şi aa atunci suprimind 
extremităţile se obţine un drum hamiltonian minim în graful G 
inițial cu mulţimea, de virturi X. 

Într-adevăr, presupunînd prin absurd că drumul hamiltonian 
astfel obţinut în graful 6 nu ar îi minim, atunci l-am putea înlocui 
printr-un drum hamiltonian mai scurt, obţinînd un drum bamilto- 
nian cu extremităţile ua Şi Sata în graful cu mulţimea, de vîrturi A 


- care ar fi mai scurt decît cel minim, ceea ce este absurd. 


Aceasta, deoarece toate muchiile [tara il Și [ra 2] CU 
i, j < n au respectiv aceleaşi lungimi (fig. I.2.15). 

Deci problemele P1 şi P2 sînt echivalente din punctul de 
vedere al rezolvärii, în sensul dat anterior. Vom arăta acum că ştiind 
să rezolvăm P2 putem rezolva problema P3 astfel : 

Adăugăm un nou virf, notat #a+ı la mulţimea X de virturi a 


graului iniţial eu matricea costurilor (lungimilor) 0. 
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Vom alege un vîrf oarecare g; în mulţi i fin 
Yo! k i, mulţimea X şi vom ini 
seu O LOT [lny 4] astfel: lungimea, muehiei ei z] 
i 3 wi iei f : i Ă ] 4 

ed ii gimea muchiei [&, z;,] pentru orice i = 1, 2,... 


o E unra vine la completarea matricei costurilor sau lungimilor 
S ă o linie și coloană cu Cosa; = Cangi = Gia pentru i=l 
2, sn $ : tao te DEGI A la 


Pentru acest nou graf să determinăm un DHM cu extremităție - 


fixate, a, şi æ Presu î 
vat i n+a. Presipunînd că în acest drum hami i 
7 AE j amiltoni: 
ra yael anterior, lui Bagi este m, să suprimăm muchia 
Eain s] pi saro Diloeniin prin muchia [;, a, ] de aceeaşi lungime 
torm definiției date. Astfel se obține un circuit hamiltonian în 
gra; a inițial 6. Să demonstrăm că acesta este un CHM (X, 0) 
P E e = 7, j 
SER eri amp prin absurd că acela nu ar fi un circuit hamil- 
onig „ R înioeuim prin altul de. lungime strict i mică, 
Notînd cu æ, virful anteri i î r irent, k suprimat 
€ A rior lui æ, în acest circuit, să imăm 
muchia [a &;, ] şi să o înlocui ri i e FA si 
i a A ocuim prin muchia (4, æ i 
lungime, în graful cu mulţimea, vîrturilor X- LPa Terah de areas 


Am obținut astfel un drum hamiltonian cu extremitățile a, 


ȘI S, mai scurt decit cel minim, ceea ce este absurd. Deci problema 


P3 se rezolvă ca o problemă P2, însă, într-un alt graf. - 


Să arătăm acum că ştiind să rezolvăm problema P3, ştim să 


rezolvăm problem: t rg i 
apă DP ri pu ceea ce va demonstra echivalenţa proble- 


Într-adevăr, să, adăugăm la graful Œ cu mulțimea, vîrfurior X. 
> 


un nou virf 4, pe care îl legămcu vi ile di i 
Ș à 7 rfuril îi 
avînd toate o aceeași mngime i > 0 (fig. Laan TRE pi Aa 


Găsind un circuit hamiltonian minim în graful cu mulțimea, de 
vârturi Xa = X U (asa), prin suprimarea muchiilor care îl unese 
pe Sn CU cele două viriuri din X se obține un DHM în graful 
iniţial. Demonstrarea  minimalițăţii drumului hamiltonian astfel 

Xn 


Xatt 


xX. 
Fig. 1.2.16 Fig. 1.2.17 
. obținut se face tot prin reducere la ábsurd, analog cu demonstratiile 
de mai înainte, deoarece din orice drum hamiltonian în X se obține 
un circuit hamiltonian în X,, de lungime mai mare cu 2a, prin 
legarea, extremităților sale cu noul virf aia, prin muchii de 
lungime a. A 

Vedem deci că cele trei probleme P1, P2 și P3 sînt echivalente 
din punctul de vedere al rezolvării. Vom prezenta mai departe un 
algoritm pentru rezolvarea problemei P1, în cazul grafurilor neorien- 
tate. Problemele P1, P2, P3 se pot enunţa şi în cazul grafurilor 
orientate, demonstrarea echivalenței lor fiind analoagă cu cazul 
neorientat, dînd noilor arce introduse în graful cu mulțimea vîriu- 
rilor X, o orientare convenabilă în toate cazurile. 

Un exemplu de determinarea drumului de timp mazim într-un 
graf îl constituie găsirea drumului critic într-un graf de activități 
(fig. 12.18). 

În acest; graf arcele reprezintă, etapele (activităţile) unei lucrări 
complexe, compusă din mai multe activităţi, iar numerele încercuite 
asociate arcelor reprezintă timpul de desfăşurare a activităţilor 
asociate arcelor gratului, momentul începerii activităţii fiind repre- 
żentat de extremitatea inițială, şi momentul terminării activităţii 
fiind reprezentat de extremitatea finală a arcului respectiv. i 

Începerea, întregii lucrări este reprezentată prin vîrful A, iar 
terminarea, ei prin vîrful I. | 
: În această lucrare anumite activități se pot desfășura în acelaşi 
timp. Arcul dintre B și G este desenat cu linie întreruptă, deoarece 
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operaţia, reprezentată, prin: arcul (G, I), nu poate începe înainte de 
terminarea, operaţiei reprezentate prin arcul (4, B). Timpul asociat; 
acestui ard se consideră, că, este egal cu zero. 


Drumul de timp mawim dintre A şi I se numeşie drum critic în ` 


graful activităților din figura 1.2.18 şi el este desenat cu linie îngro- 
şată. Timpul asociat acestui drum, egal cu suma, timpurilor asociate 
arcelor sale, este de 17 unităţi (de exemplu zile, săptămîni, luni) și 
el reprezintă timpul de realizare a, întregii lucrări. . 


Fig. 1.2.18 


Drumul critic reunește activităţi de a căror realizare la timpul 
prevăzut depinde realizarea la timp a, întregii lucrări. Orice mărire 


a timpului asociat arcelor drumului critic conduce la mărirea, timpului ` 


de realizare a întregii lucrări, de aceea activităţile situate pe drumul 
critic din graful activităților trebuie supravegheate cu deosebită 
grijă pentru æ nu depăși timpul planificat; de realizare a lucrării, 

Pentru a preciza mai bine cele spuse, vom numi program o 
mulțime de operaţii care conduc la realizarea unui obiectiv (proiect, 
ansamblu de construcții etc.), astfel încît se cunose pentru fiecare 
dintre operaţii durata, sa, și relațiile de ordine temporală, privitoare 
la aceasta, (antecedente obligatorii). 


Un program poate: fi reprezentat printr-un graf orientat, numit ` 


graf de program, ale cărui arce reprezintă, operaţiile, cum este graful 
din figura I:2.18. Duratele operaţiilor sau timpii operatori fiind 
presupuși cunoscuți, se va atașa arcelor grafului o valoare ne-negativă 
reprezentind fiecare din aceşti timpi operatori. 

Virturile grafului, numite evenimente, pot fi interpretate ca 
indicînd realizarea, obiectivelor parţiale (de exemplu, montarea unui 
subansamblu). 
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CR TRAI ca EAI FEED Aa STITI 


Durata operaţiei P,, care este una dintre operaţiile care con- 


` tribuie la realizarea evenimentului H, va fi notată cu t}, valoare 


asociată arcului (£;, E,) în graful de program (fig. I.2.19). 


Fig. 1.2.19 


ili lui necesită pentru fiecare operație cunoașterea 
N eco ea d nemijlocit, graful traducind. relațiile de 
ordine în timp existente în mulțimea operațiilor. E 
Graful unui program are o proprietate importantă, și anume 
aceea că este fără circuite, deoarece în caz contrar o aceeaşi operație 
ar putea, fi făcută să urmeze ei însăşi. E af, i Sp: 
Primul lucru de îndeplinit în vederea, aplicării metodei drumului 
critic va fi întocmirea unei liste de operaţii pentru care timpii opera- 
tori să fie de valori comparabile. ROR - 
Apoi se vor determina, antecedentele fiecăreia dintre operaţii. 
Această parte, numită analiză de lucru, este cea mai importantă, și 
esea cea mai lungă. sa A 
i Vom da dude de situații care pot fi întilnite în alcă- 
tuirea grafului unui program. i | 
ui copiile că între două evenimente E și T, se pianonun 
două operații diferite b şi ¢ care urmează unei operații æ şi sin! 
urmate de o operație d. i e 
Putem reprezenta, această situație ca în figura 1.2.20, gia 
care reprezintă un multigraf. Pentru a lucra cu grafuri obişnuite, 


RUR : ea 
. vom introduce un eveniment fictiv Ez şi o operaţie) virtuală (de 


timp zero) între F; şi E; (fig. 1.2.21). 


Fig. 1.2.21 
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Să considerăm acum situația din figura I.2.22 ii i 

d urmează am ele lui a şi d. Să Prenume Ri ae ai 
„urmează, lui a, d îi urmează, lui b, dar e nu îi urmează lui b. Nu se 
mai poate folosi același graf, trebuie introdus un eveniment; fictiv 
s ȘI o operaţie fictivă, desenată cu linie întreruptă în figura I.2.23 


b . 
&) EG) 


“Fig. 1.2.23 


Fig. 1.2.22 


aa ali graful unui program, cu începerea lucrărilor 
prezenta ao S ieste E, şi încheierea lucrărilor reprezentată 
Sa en m VOM căuta, data realizării ansamblului de lucrări 
adică durata programului de realizat. ? 


ipali ge Saol zi ai ea a a aO oni 
gare i To peost doug vîrfuri o sumă maximă de timpi operatori. 
n i um (care poate să nu fie unic) se numeşte drum 
Luînd pentru dur i ări impi 
operatori luaţi Bi ram ca ra bak ie i A "a PA ne-am REEN 
EA rat a e pot fi efectiv realizate, ţinînd seama, 
Calculul datei de realizare a evenimentului E, revine deci la 


căutarea în graf a drumului cel mai lung S drumului critie 
e 
al lung sau a wu: 


Pentru graful din figura 1.2.18 vom atribui evenimentului inițial - 


A termenul 0.'Cel mai lan; 

l 0. g drum de la A la Z sau drumul criti 
sli „desenat îngroșat pe figură şi el are lungimea, de 17. ei 
umăr reprezintă timpul de realizare pornind de la data zero. 


Deci dacă programul se desfășoară fără incidente, durata sa 


va fi de.17 unităţi de timp (de exemplu luni jii 
va 7 ; plu luni). Operaţiile reprezentate 
aa y A Oh (0, F), (F, 6), (6, ) ale drumului critic se numesc 
Lungimile maxime ale drumurilor de la vi te 
A i ime i a virful A la; celelalte 
pui B, 0, D, F, F, G, H, I ale grafului din figura I.2.18 au 
a recute pe desen în pătrăţele. Ele reprezintă datele aşteptate de 
realizare a evenimentelor reprezentate de vîrfurile graiului. 
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ta 


Scan orienta a e 


Să, observăm că, operaţiile critice trebuie să, înceapă la datele 
evenimentelor origine ale arcelor. Astfel, operaţia reprezentată de 
arcul (0, F) trebuie să înceapă la data evenimentului 0, àdicà la 
data de 4. Dacă pornirea unei operații critice este întirziată, atunci 
tot programul va întîrzia. De exemplu, dacă operaţia, reprezentată 
de arcul (F, G) nu începe la 10, ci la 13, lucrarea în intregimea ei 
va, întârzia cu 3 luni, fiind terminată cel mai devreme la daia de 20. 

Operaţiile critice sînt, deci, acele operaţii ale programului asupra 
cărora trebuie să ne îndreptăm toată atenţia, căci de buna lor 
desfăşurare depinde data de realizare ansamblului lucrărilor. 

Vom vedea, în continuare, că operațiile necritice admit unele 
întârzieri în începerea lucrărilor. 

Datele aşteptate ale diferitelor evenimente F, care compun 
programul, notate cu t; se calculează evaluind cele mai lungi distanțe 
de la evenimentul origine E, la celelalte evenimente reprezentate 
ca vîrturi ale grafului programului. Este acum interesant de cunoscut, 
pentru fiecare eveniment E, neeritie data sa limită de realizare, 
dată de la depăşirea căreia tot programul va fi întirziat. 

Să considerăm evenimentul Æ;. Timpul minim necesar pentru 
realizarea 6peraţiilor situate între F; şi Æ, se obţine căutînd în 
graf drumul cel mai lung de la F; la Fn, deoarece trebuie să fim 
asiguraţi că operațiile care urmează lui E, pot fi toate realizate, 
tinind seama de duratele lor. 

Data limită tf a evenimentului F; se va obţine deci scăzind 
din data t, a evenimentului final E, lungimea celui mai lung drum 
de la E; la En i 

Să, observăm că pentru evenimentele critice, fiind situate pe 
cel mai lung drum de la E, la En această diferență este egală chiar 
cu t, lungimea celui mai lung drum de la E, la E. s 

Într-adevăr, dacă E, se găseşte pe drumul cel mai lung de 
la E, la E, atunci lungimea maximă a ărumurilor de la Æ, la E, 
este egală cu suma dintre lungimile maxime ale drumurilor de la 
DB, la E, şi de la E, la E. 

Am obținut; astfel, pentru fiecare eveniment, două date : 

— data t; este data așteptată a realizării evenimentului £,; 

— data t* este data limită de realizare a evenimentului E, 
după depăşirea căreia timpul total de execuție a ansamblului lueră- 
rilor este majorat. . 

Data î, este numită şi data cea mai apropiată, iar data limită 
tf este numită şi data cea mai târzie. 

Am văzut că pentru evenimentele critice data limită 4% se 
confundă cu data aşteptată t. 
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Aceste evenimente critice nu admit nici o întârziere în ceea ce - 


priveşte începerea lucrărilor. Ă 


Intervalul (i, 1] se numeşte interval de flucluaţie. El este . 


deci intervalul în care poate începe realizarea evenimentului necritic 
Ea fără a modifica timpul total de execuţie a ansamblului lucrărilor, 
În cazul evenimentelor critice, acest interval se va reduce la un 
singur număr, deoarece t; = t*. i 

Să reluăäm exemplul grřafului din figura I.2.18. Pentru calcu- 
larea datelor limită ł¥ putem proceda astfel: Din virful H pleacă 
un singur arc la T, care este un eveniment critic. Lungimea, acestui 
arc fiind 2, rezultă că tă = 17—2 = 15, intervalul de fluctuație 
al lui E find [12, 15]. i 

Din E pleacă un singur drum pînă la I, de lungime 6, deci 
tă =17—6 = 11, intervalul lui Ẹ fiind [7, 11]. Pentrù D obţinem 
i = 10, iar t$ = 14. În. final, datele evenimentelor A — I se vor 
putea plasa în intervalele următoare, aceste intervale reducîndu-se 


la o singură dată obligatorie în cazul evenimentelor situate pe. 


drumul critic : 


A:0 P.: 10 
B: [3,14] . . Gil 
i C:4 H : [12,15] 
D.: [6,10] 1:17 
B: [7,11] 


Vom fi interesaţi de asemenea, să cunoaştem, pentru fiecare operaţie 
O; reprezentată de arcul (E, B,) în graful programului, întîrzierea 
care poate fi admisă la începerea sa, fără a perturba data așteptată, 
de realizare a evenimentului £,. Această, întîrziere, numită marginea. 
liberă a operaţiei 0, este egală cu: . - 


îs — bn — tys 


unde 4 şi î, sînt datele așteptate ale evenimentelor B, și £,, care 
încadrează opeiaţia Oy. i S 
Întivzierea care poate fi admisă la începerea operaţiei O, 
fără a perturba data limită de realizare a evenimentului E, se 
numește marginea totală a operaţiei O; şi ea este egală cu: 


E i A E 


Aceasta este întîrzierea maximă, care: poate fi admisă la începerea 
operaţiei O; fără a. mări durata, întregului proiect. 
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.  citatea unui program. 


i 


ERN 


Marginile libere și marginile totale ale operațiilor critice sînt , 
nule, deoarece începerea lor trebuie să se facă fără nici o întârziere. 
a Deoarece iř > 4, marginea liberă este mai mică sau egală cu 

i totală ricărei atii a programului. 
marginea, totală a oricărei operaţii a prog aulai, _ 

Sntervalele đe fluctuație şi marginile operațiilor măsoară elasti 

Cu cît acestea sînt mai reduse, cu atât „programul este ma 
rigid. Un program total rigid este un program în care toate drumurile 


i sînt critice : nu va exista nici un interval de fluctuaţie, iar marginile 
- libere şi foţale- vor fi nule. 


lt aspect; al noţiunii de margine liberă poate fi examinat | 
în Elea a și care timpii: operatori ai operațiilor necritice pot îi 
iaca AnS liberă corespunde atunci măririi posibile 2 durate 
operației necritice Oy, lăsînd neschimbate datele t şi 4, ale eveni- 
mentelor E, şi E; care încadrează această operație. Ne iaa 
De exemplu, pentru graful din figura 1.2.18, marginea liberă 

-a operaţiei reprezentate de arcul (D, E) este: 


te —tp— tip =T —6—1=0, 


în timp ce marginea sa totală este: 


H tor =11—6—1=4. 


E to 


i rația (D, E) poate admite o întârziere: maximă de 4 luni, 
mea eg întivztea execuţiei întregului proiect, reprezentat 

g din figura I.2.18. , 
z E aaa intervalele de fluctuaţie ale evenimentelor, precum 
şi marginile totale ale operaţiilor în paranteze, graful din figura 
12.18 devine graful din figura [.2.24. e a PI 

Importanţa notiunilor de intervale do uctiakie şi de margini 

are în cursul realizării programului. , MECN 

vale să toate evenimentele încep înaintea datei. lor limită Li 
programul se va desfăşura normal şi data realizării finale în n E i 
respectată. Dacă însă se întîmplă ca, dintr-un motiv oarecare, data 
de pornire a unui eveniment să depăşească data „limită i, atiunei 
trebuie să se accelereze sarcinile situate pe cel mai lung drum w la 
E, la E, Dacă nu există această posibilitate, vom putea evalua 
întîrzierea datei î, de realizare a ansamblului lucrărilor. A 

Datele limită, sînt deci limite de securitate şi joacă rolul de 
sonerii de alarmă pentru conducătorul lucrărilor. în timpul punerii 
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în execuţie a programului. Un exemplu practic, Gin lucrarea [6], 
este tratat de problema 26. . 

Metoda drumului critic este mult utilizată în urmărirea tucră- 
rilor de construcții-montaj. i i 


e BA 


4], po 


E 
a C em F © 


Fig. 1.2.24 


Pentru găsirea drumului sau a drumurilor critice într-un graf 
există algoritmi care reduc această operaţie la un calcul sistematic, 


‘care evită calculul timpurilor așociate tuturor drumurilor din graf. 


Dacă interpretăm numerele asociate acelor unui graf ca, nişte distanțe 
dintre două virturi, atunci drumul de timp maxim reprezintă drumul 
de lungime maximă dintre cele două vîrfuri, lungimea unui drum 
fiind egală prin definiţie cu suma lungimilor arcelor (muchiilor) 
care compun acel drum. Numărul asociat unui are (unei muchii) 
mai poate avea şi semnificaţia unui cost al transportului dintre 
donă virfuri. 

În probleme practice sîntem interesaţi în determinarea unui 
drum de lungime minimă sau a celui mai scurt drum dintre două 
localităţi, care pot; fi reprezentate ca vîrfuri ale unui graf, numărul 
asociat muchiei care le unește fiind lungimea în km a drumului 
auto direct dintre cele două localităţi. De exemplu* pentru graful 


din figura I.2.25 în care virfarile reprezintă oraşe, iar numerele . 


asociate muchiilor reprezintă distanţele directe (în sute km) dintre 
aceste oraşe, drumul minim dintre A şi H este wrmătorul : i 


[4, 0, E, F, H] 


care are o lungime egală cu 10 unități (1000 km). 


` * vezi și problema 11 
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eane oa Di 3 manba i EE 


Dacă îi asociem muchiei [D, F] numărul 4, adică între D şi P 
se construiește un drum mai scurt, de 400 km, atunci între A și 
H mai există un drum minim de aceeaşi lungime, și anume: 


[4, B, D, F, H] 


OR 


Ð 
(0) 


Pentru calculul lungimii minime a dramurilor de la un sia 
notat æ, la toate celelalte vîriuri da, Bay- +,%u ale ni em j 
vom prezenta un algoritm cuie il IA anume algoritmu 

i Dantaig.care constă în următoarele: 
e e notează cu i funcția distanța minimă de la a la celelate 
vîturi, adică Ka) Zi (RR minimă a drumurilor de la virtu 
îrful a, al graiului. i 1 , 
“ zi Xa se nateeză o submulțime formată din q virfuri (ara) 
ale grafului G = (X, U) cu funcția de distanță l: UR, care 
asociază fiecărui are (muchii) lungimea sa. 

Algoritmul are următorii paşi : 

Pasul 1. Se notează X, = (a) şi se defineşte iz) =0; A 
Pasul 2. Dacă, s-a determinat mulțimea Xg =n, iy one Ciad 
pentru fiecare virf zip din X, se determină acel virf notat Yip care nu 
aparține lui X, şi în plus este cel mai apropiat de zip „ adică lungimea 
arcului (Sip, Yip ) este minimă în mulțimea, lungimilor arcelor care 
pleacă. din Sip către vifuri care nu aparțin mulțimii X, A | 

Dacă, există mai multe vîriuri care se leagă de Xip prin arce de 
lungime minimă, se va alege unul oarecare dintre ele. 

Se determină acel virf y; cu proprietatea că : 


(i) + Via Yie) = min (Vy ) + Uz Yi )) şi se definește X; = 
= X,U {Yu}, definindu-se funcţia t în Yi; prin : 


UYis) = HDi) T Uliss Yis) 
Se repetă pasul 2 pentru q = 1, 2,.. 9 — L 
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i La sfirşitul aplicării algoritmului se -obține valoarea funcției 
distanță minimă pentru toate virfurile Va, g.. 2, adică distanțele 
minime de la virful a, la celelalte virturi ale prafului. i K 

„Dacă funcția i nu poate fi definită într-un virf &,, înseamnă 
că nu există nici un drum de la v, la æ, şi vom lua prin definiție 
Hep) = 0. > ; i 

Se observă că acest algoritm construieşte prin recurenţă, în 


mod sistematic, distanțele minime de. la æ, la celelalte vîrturi ale ` i 


graiului. - 
f Acest algoritm se poate aplica şi grafurilor neorientate, schim- 
bind în descrierea sa cuvintul arc prin cuvîntul niuchie. 
Să, considerăm exemplul grafului din figura I.2.26. 


Fig. 1.2.26 


Aplicînd algoritmul Di obținem : 


X, = {m}, Wa) = 0. Yy cel mai apropiat vîrt de m, este unic deter- 


minat și el este 2. Vom defini deci X, = (au, 2a} şi (£a) = (a) + 


Ho 2) =0 4+1 =.. 
Găsim acum : Y; = da; Va = Vz SAU dy. 
Luăm de exemplu yı = £z, Ya = Wz- 


Deoarece (2) + Han 83) = 0 +3 >l) + Ula 2) = 1 +1=2 
alegem ca. nou vîrf pe pă și definim î(23) = 2. Ca A 

Deci Xa = (up. V} Acum cele mai apropiate vîirfuri din 
afara mulțimii X, vor fi: 


Yı = a Ya = Vz, Ys Nu există, 


Deci avem obligatoriu X, =, o By & 
= min (H0) + Wep w) Wa) + kas a) = min (2,5 


l% 
] 


j 
) 
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„Acum găsim Yı = 84; Ya Ya, Ys DU există, deoarece nu există 


arce de la Dzn Lz ȘI V; către aa. 


Deci Xy = {2u Vo gy Va as} Şi aa) = ia) + 4 = 4. Am 
obținut astfel distanțele minime : E 


ia) =0; Ue) =1; Hs) =2; 


Ha) = d; Wa) = 2. 


„Dacă notăm ay = Uz æ) pentru î, j = 1... distanţele 


” directe dintre vîrfurile grafului G, matricea : 


A = (yhet... 


i - se numeşte matricea distanțelor directe a grafului G. În această 


matrice, a,=— co prin definiție, dacă nu există arcul (i, j) în graful G, - 
ceea ce se mai poate serie (i, j) ¢ U. 

Vom nota cu ağ lungimea minimă a drumurilor de la virful 
e, la virful x al grafului, cînd există măcar un drum de la s la 
a, în graf. Dacă nu există nici un drum de la la 2 vom pune 
prin definiţie e = co. Matricea A 


A* = (CAE AP 


se numeşte matricea distanțelor minime ale grafului G. Prin. definiție 
avem Gu = ağ = 0 pentru i = 1, 2j... m% . 

Algoritmul D1 al lui Dantzig determină deci toate distanțele 
minime de forma ağ pentru i= 1,...,, deci o singură linie a 
matricei A*. Pentru găsirea întregii matrice A*, adică pentru deter- 
minarea distanțelor minime: dintre oricare două vîrturi àle grafului G, 
se poate aplica de ori algoritmul Dl. Un alt algoritm pentru 
găsirea tuturor distanțelor minime eğ este algoritmul D2 al lui 
Dantzig, care are o rază mai mare de acţiune, putînd detecta și 
eventualele cirenite negative ale gratului, presupunînd că funcţia 
distanţă 1: U —> R poate lua și valori negative. 

Este clar că dacă un graf admite un circuit; de lungime totală 
negativă, atunci distanţa dintre oricare două virfuri ale circuitului 


poate fi făcută negativă şi oricît de mică parcurgînd de un număr 
` suficient de mare de ori circuitul, deci nu există o cea mai mică 


distanţă între vîrturile circuitului. 

Se fac mai întîi următoarele notații : 

Se notează Xn = {8u ae) Și cu ag- distanța minimă, 
între a, și æ (i, j < m— 1) în subgraful cu mulţimea de virturi 
Za; iar cu d? distanţa minimă între æ; şi æ; (i, į < m) în sub- 
graful cu mulțimea de vîrturi Xn. 
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[i 


j 


i : aii imă ai ă există un drum de la 
Pasul 1. Se defineşte d? = aa; dă, = az, unde ay sînt distanţele Inegalitatea 3 exprimă tocmâi faptul că exis 


; - î patul cu mulţimea, 
i directe dintre vîrturile æ; şi æ; ale grafului dat. "a, la Em, urmat'de un drum de la 7 la æ, în subi i) 


; Š i " jl imilor este negativă, adică 
aip Dacă a + au < 0, stop. S-a detectat un circuit negativ. „de virfuri Am D a îi apa ome cui d i 
; i i oemai un circui 1 Ş „Pi A 
„Pasul 2. Dacă s-au determinat distanțele minime d, distanţele 5 Egalitatea 4 folosește proprietatea că drumul ia e e n = 
minime d} se determină astfel : i - : în subgraful cu vifurile din Xn fie că trece p m 


i ns caz lungimea sa este egală cu die + d ue tit sirere 
ii 2 (diam LESMI prin virîul 2, cînd lungimea, sa este dy”, el fiind un 
see m=1 - în subgraful cu multimea đe virfuri Ama: ee 
i s i. Să considerăm exemplul grafului din figura I.2.27. 
2. du =min -o (Om HT sism=l. ; ? 
i d= s‘ ; 
ł 


8. Se verifică dacă există un indice i cu proprietatea că 1 sis 
<mn-—l şi ; i 

| di, + dm, < 0. 

j „ Dacă da, stop. Există un circuit negativ în subgraful cu mulțimea, 
© de vîrfuri Xn care trece prin vîrfurile Bi ŞI Bm 


p X. 
3 Fig, 1.2.27 


4. Se calculează : 
| l d% = min (434, da + a) 
pentru i, j = 1, 2,...m—lşiigj. 


de =1, d = ©; dis = 


finim là început X, = {2p #2}; 
n L2)=2; = 


= min (d) + da Be + az) = min (0 -+ 2, 14 

min (d, + da da + as) = min (%0 +2, 0+2) = 2. 

Da fel se găseşte că d, =co; d =c0; 

die = min (Bo, da + d3) = min (1, 00) = 1; 

Bi = min (dy, d + dă) = 00. | îm 
Pentru graful din figura 1.2.27, găsim distanţele minime : | 

l- ag) = min(o0, 2, 3) =2; 


Se repetă pasul 2 pentru m = 3,.. Re . : 
| La sfârşitul aplicării algoritmului distanţele minime d, vor fi f 
| distanţele minime în graful cu mulțimea de virfuri X = (ay ae 
li «--3%a}, deci vor fi tocmai ağ. ` A 
| iii Dacă graful nu conţine arce de lungime negativă, punctul 3 
i] de la; pasul 2 va îi suprimat. Dacă graful conţine arce de lungime 
| negativă, la 3 găsind un circuit negativ ne oprim, deoarece am văzut; 
i că în acest caz nu există, distanțe minime, pentru unele perechi de - 
“| vîriuri putînd avea distanțe oricât de mici (nemărginite inferior). 
Se presupurie de asemenea că d” = 0 pentru orice è și m. 
Algoritmul D2 determină recursiv distanțele minime în subgra- 
îurile eu mulțimile de vituri (e, Sa}, {2u a Pola {ln Las Cay da), | 
ee {2u Dope ct = X. ada, 
Justificarea acestui algoritm se poate face astfel : Egalitatea 1 fi 
de la pasul 2 exprimă faptul că distanța minimă de la m; la £, în. 
subgraful cu mulțimea de virfuri X, este egală cu minimul dintre | 
distanța directă a, (pentru j = i) şi sumele distanțelor minime 
de la v, la un nod arbitrar w; în subgraful cu mulţimea de virfuri $! 
Ama cu distanțele directe am. Egalitatea 2 se justifică în mod 
analog. i È 


až, = č, = min (É, + ta die + aa dis - 


Snin 3)=1; 
aă = di, = min (d, Mia dz + aa das + asa) = min (0, 1, ) 7, 


- aa) = min (co, œ% 1)=1; 


Á 3 
až, = di = min (dh + a, Be + Oza dă ṣ 


x 43 3 J 
a, = dh = min (aa + fh, 42 t dzis, da3 


- 3) = min (2, 00, œ) =2; 


f Symi œo) =3:; 
ah, = dh = min (an + des aa + Bas ass + die) = min (3, co, %) d 


+ 48) = min (4, o, %0) =£; 


a% = Qi; = MİN (aa -+ dis, aa + dis, ao 7 
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at; = din = min (dhe, dh + dł) = min (1, 2 + 3) = 1; 


ah = dis = min (dis, di + dh) = min (2, 2+ 4) = 2; 


ah = dis = min (dh, di + dh) = min È, 1+4) =2; 


af = i, = min (Aa, d + di) = min (o, 1 +2) =3; 


ai = da = min (dhs, di + dh) = min (%0, 1 +3) = 4; 


ai = di = min (dhi; dh + dh) = min (00, 1+2)=3. 


Pentru exemplul dat, care este un graf cu numai patru virfuri, 
aceste distanțe minime se pot găsi mai uşor urmărind desenul gra- 
fului, ceea ce nu mai este cazul pentru grafuri cu un număr mare 
de vîrturi şi de arc. i : - i 

Rezultă că pentru graful din figura, 1.2.27, matricea A a distan- 
telor directe, respectiv matricea A*.a, distanțelor minime stai 
următoarele : a 


01 2 œ 0 12 2 
a=|9 0 2 1] 4s _ [3 0 2 ı 
co co 0 1 3 4 0 1 
2 © co 0 2 3 4 0 


Pentru determinarea drumurilor elementare de lungime minimă, 
în cazul în care cunoaştem matricea; distanțelor minime A* sau 
numai prima linie a acesteia, se poate folosi un algoritm simplu, 
care constă în următoarele : 


Presupunem că vrem să găsim toate drumurile elementare de 
lungime minimă de la vîrful d la viitul a, (i > 2) al gratului cu 
matricea distanțelor directe A. i 

Dacă. df co vor exista drumuri, deci şi drumuri minime 
între a, şi a. i 


Dacă ai; = au, arcul (a, z,) este un drum de lungime minimă 
de la a, la m. 


Dacă există un indice k diferit de 1 şi î care verifică ecuaţia. : 
x k 
aii = ir F ago 


atunci un drum de lungime minimă de la a, la æ, este de forma 


(Ham 20, unde u = (22: - pp) este un drum de lungime minimă 
de la a, la a. , 5 
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iames onne 


PETE, 


E 


i iile, diferite dè 1 şi î, ale 
acest caz, fie ky ko... soluţiile, c d i 
Sie scrise anterior. în această sinapis, sm, EE îi e o l 
awilor minime de la a, la & sint respectiv : (Sry i), (d 
iih æ). Repetînd procedeul indicat, găsim arcele rupt 
sior minime care ajung respectiv În Pry Dra seo Gks JAMA, 
ă, câ jungem în m, A NI d. 
pei loa că aceste drumuri minime vor F reconnu n 
i 's ori vii lor, găsind prin recurenţă toate ar ile 
sensul invers orientării arcelor, nT enta toate aronimte 
re minime de la virful a, la virful æ; E 
era ES 1 şi d care să verifice . ecuația: an = a%, + asi 
iri singurul drum de lungime minimă ya la i la îi Sate arcul 
{ emplu, să considerăm graful din figura 1.2.28. P 
ia 2), altaQ disantele minime de la vîrful æ, la celelalte virfuri, 
găsim că ele sînt egale respectiv cu: 


2, 0% 4 dia 5. 


* LA =, a 
ah = 1, a = 2, al = 3, ais = 3, ie 


Pentru a găsi drumurile minime de la a, la ~, aplicăm algo- 
ritmul deseris anterior. , . 

Găsim mai întîi că ecuaţia : 

a = ak, F arg, SAU: 
5 = alr F ara f 

l ii, şi = 6 şi k = i ultimele arce 
sle ela Aa eu ea sai r a (Bes e) şi 
(7 e) (fig. I.2.29). ; 


DAEA 
" Dho A E 


Vom determina drumurile minime care ajung în a. Rezolvăm 
ecuaţia : , | 
ae = Air Que Sau: 


2 = aă + are 
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şi găsim soluţia unică (diferită de 1 şi 6): k = 
Rezolvînd oa ic di 


E casa f 
als = Git + arg, Bau: 
5 
— g* 
3 = a + aus 


găsim k = 3 gi rezolvînd ecuaţia: : 
+ E 
alr = aik F ar Sau: 
4 = &i + a 


găsim k = 4 şi k = ð. Am ajuns deci la gr in fi 
> i n graful din figura I.2.30: 
pop sina aceeaşi operație pentru virturile 42, a şi &, găsim : 
FE entru v, : nu există k # 1 și 2, deci arcul (a, &2) este singurul 
w zoini de la m, la a; i 
entru 4: nu există k și i an i 
ürun A A air ia 1 şi 8, deci arcul (z,, #3) este singurul 
entru z,: nu există k și i ar ste si 
dr minm iși ata ist 1 şi 4, deci arcul (a, 24) este singurul 
Ă m obţinut astfel graful Qin figura, 1.2.31,; 
- i astfel „2. care est 
porpal ralu inițial, ci proprietatea că orice dam. are 
la iriuri este un drum minim, în particular t rumu- 
rile de la a, la æ, vor fi drumuri minime, Şi Amine 3 Gi ii 


E. . ȘI . i & 
(Ur, Bor Ze 22); (Lry Bar Biy e); (Lair Vas d Lr Ba) Și (Lu Lay Ey Ba). 


Fig. 1.2.30 


Fig. I.2.31 


În general, graful obținut în acest mod es 
i ţ te un i 
al graiului de la care am plecat, deoarece el se Saua Ie 
inițial prin suprimarea anumitor arce şi a anumitor vârfuri. 
„9 proprietate importantă a acestui subgraf parțial este aceea, 
că el nu conţine circuite, existenţa, circuitelor ducînd la mărirea, 
lungimii drumurilor din graf, care astfel nu mai sînt minime. 
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Justificarea acestui algoritm de determinare a drumurilor 
minime se bazează pe proprietatea următoare : Dacă există relaţia : 


at, = ale + au 
atunci există un drum minim de la v, la æ; al cărui ultim are este 
arcul {Sy &;). Într-adevăr, această egalitate ne spune că distanța 
minimă de la æ, la æ; este egală cu distanța minimă de la s, la 
æ, plus lungimea arcului (Sr, £;). i 


1.3. Grafuri conexe 


În acesi paragraf ne vom ocupa numai de graturile neorientate. 
Gratul din figura 1.2.1, care reprezintă o reţea de comunicaţii, 
are o proprietate importantă, şi anume aceea că din oricare virf 
al său se poate ajunge în alt virf, deplasindu-ne de-a lungul unui 
drum al grafului. Un astfel de graf se numeşte graf conex. 


Definiție: Un graf G care are proprietatea că oricare două 
vârfuri ale sale sînt legate printr-un drum din G se numeşte graf cones. 


Noţiunea de conezxitate are importanţă în aplicațiile practice. 
Astfel, graful care reprezintă rețeaua de comunicaţii auto sau tele- 
fonice dintr-un oraş trebuie să fié un grai conex, asigurîndu-se 
astfel legătura (directă sau indirectă, prin mai multe vîrfuri inter- 
mediare), dintre oricare două puncte ale sistemului. 

În cadrul grafurilor conexe un loc aparte îl ocupă arborii. 
Graful desenat în figura 1.1.7 este un arbore. 

Prin definiție, un arbore este un graf neorientat, conex și fără 
circuite. ; 

Un arbore mai poate fi caracterizat şi prin proprietatea că este 
un graf neorientat, conex, avind un număr minim de muchii, 
adică dacă îi suprimăm o singură muchie, se obţine un graf care 
nu mai este conex. ` , A 

Astfel, pentru arborele din figura I.3.1a, suprimarea unei 
muchii îl face neconex (fig. I.3.1b). Arborele din figura I.3.2 se 
numește stea. Pentru arbori, între numărul vârfurilor n și numărul 
muchiilor m există o relație simplă, dată de teorema următoare : 


Teorema 1.3.1. 


Dacă G este um arbore cu n %rfuri şi m muchii, atunci m =n — 1. 
Pentru a demonstra acest rezultat, să arătăm mai întâi că 
pentru un arbore A există cel puţin un vîrt de grad 1 (numit și 
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virf terminal), deci un vîrf la eare sose: i că i a 
și un y şte o singură much: 

Ba alt ea pauu Să presuponen contrariul, S anume Că tarare 
ğ ratului are gradul d(z) > 2. Fie æ, un virf al i 

El se uneşte prin muchii cu cel puțin alte două virfuri, æg a DS 


a) &) 
Fig, 13.1. 


Fig. 13.2 


alwa) 22, deci 4, se uneşte printr-o muchie cu cel puțin un alë 
v 29. a a = o, atunci am obţinut că arborele ar conţine 
circuit cu 3 vîrfuri și anume Sp Sy d. Procedind astfel din 
aproape în aproape găsim un drum (Vo, Ly.. .; U) Virtul a, are şi 
el Uas) 22. Ştim că 4, se leagă printr-o muchie cu ®, K Dacă 
2 s-ar lega printr-o a II-a muchie cu unul din virfurile #o, Dis. -bp 
ar rezulta că arborele conţine un circuit, ceea ce contrazice definiţia 
unui arbore. Deci a, se leagă cu un nou vîrf ær, Deoarece numărul 
vârfurilor este finit; vom ajunge la situaţia cînd k =n —1 (fig 
L 3. 3), Dar virful g,_. are gradul A(önı) 7> 2, deci Sa_ı-.Se leagă cu 
Pag şi cu cel puţin unul din vîrfurile rămase, co; By... Zaa, deti 
graful conține un circuit, ceea ce contrazice definiția unui” arbore. 


Am demonstrat deci prin reducere la absurd că orice arbore 
are un virf de grad unu. Suprimînd acest; virf și muchia, care se 
leagă cu el, obținem un nou arbore cu m’ = m — 1 muchii şi n’ = 
=y — 1 virfuri, deci p—m=w'—m. -- ` BrT 
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Continuind acest procedeu obținem numai arbori pentru care 
diferența dintre numărul vîrfurilor şi numărul muchiilor este aceeaşi. 
În final, vom obține arborele cu un virf şi fără nici o muchie, deci 
pentru care diferența dintre numărul virfurilor și numărul muchiilor 
este egală cu 1. Rezultă că n — m = m RR ae 
... = 1, deci pentru arborele de la care am plecat n — m = 1 sau 
m =n — 1, eat. 

În figura 1.3.4 este ilustrată aplicarea acestui procedeu de 
eliminare a virfurilor și muchiilor în ordinea : 2, 6, 1, 7, 4, 5. 


7 
5 
Nge ON; Ri 


Fig. 1.3.4 


Demonstrația acestei teoreme se poate face şi prin inducție 
după numărul vîrfurilor, n (problema 18. 

O proprietate importantă a unui arbore este aceea, că adăugînă 
o muchie între două virfuri æ; şi œ; care nu erau legate printr-o 
muchie se creează un circuit unic. Într-adevăr, deoarece a; şi &, sînt 
yirfurile unui arbore, rezultă că există un drum care unește pe 4, cu Xy. 
Acest drum, împreună cu muchia [e 2;] formează un circuit. 

Să presupunem, prin reducere la absurd, că s-ar forma două 
circuite prin adăugarea muchiei [4 2]. Aceste circuite trebuie să 
conțină ambele muchia [s æ] Prin suprimarea acestei muchii ar 
rezulta că arborele conţine cel puţin un circuit (fig. 1.3.5), ceea ce 
este absurd. 

Faptul că arborii sînt grafuri conexe cu număr minim de muchii 
axe importanță în anumite probleme de tipul următor (problema 
arborelui minim) : 


Fig. 13.5 


__ O stație de pompare S trebuie să alimenteze cu apă consuma- 
torii A, B, 0, D; E, F. Costurile de construcție ale. conductei de 
legătură dintre stație și consumatori, cât; şi între anumiţi consumatori 
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sînt date în unităţi băneşti (de exemplu zeci de mii de lei), asociate 
muchiilor grafului din figura I.3.6. Să se determine care dintre legă- 
“turi trebuie realizate, pentru a se putea alimenta, cu apă toți consu- 
matorii şi astfel încît; cheltuielile de construcție ale rețelei de alimen- 
tare să fie minime. 

Din cerințele problemei se deduce că graful care reprezintă, 
legăturile dintre staţie și consumatori trebuie să fie un graf conex 
(pentru a se putea alimenta toţi consumatorii) şi fără circuite (pentru 
a avea un cost cît mai mic), deci el este un arbore. 

Într-adevăr, dacă graful alimentării cu apă ar conţine un circuit, 
de exemplu [B, D, F, B] (fig. I.3.7a;), costul construcţiei ar include 


Fig. 1.3.7 


suma costurilor asociate muchiilor [B, D], [D, P] și [B, F], adică 
3 +3 + 4 = 10 unităţi. Dacă suprimăm însă muchia, [D, F], costul 
se reduce cu 4 unităţi (fig. I.3.7b). Din acest exemplu vedem că 
pentru a avea un cost de construcţie cît mai mic, gratul-soluţie 
trebuie să nu conţină circuite. Pentru determinarea unui arbore de 
cost minim se poate aplica următorul algoritm ` (algoritmul lui 
Kruskal) : 


1. Se alege o muchie a grafului cu costul cel mai mie. 


2. Dintre muchiile nealese, se alege una care nu formează circuite 
cu muchiile deja alese şi care are costul cel mai mic. 2 


Se aplică pasul 2 al algoritmului pînă la obţinerea unui graf 
conex și fără circuite cu aceeaşi mulțime de virfuri ca şi graful conex 
de la care am plecat. i 

În general, dacă există mai multe muchii cu acelaşi cost, se 
poate ca problema, să aibă mai multe soluţii cu costul minim. 

. Pentru exeinplul din figura 1.3.6, se alege mai întîi muchia 
[B, E] de cost minim (egal cu 1), apoi muchia (4, E] (sau (0, D]) 
de cost 2 etc., pînă la obţinerea arborelui de cost minim (fig. 1.3.8). 
Costul asociat; acestei soluţii optime este 14. O altă soluţie optimă 
este reprezentată în figura 1.3.9, corespunzătoare. alegerii muchiei 
[0, F] în locul muchiei [B, F). 
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Dintre aceste două soluții de construcţie a rețelei de alimentare cu 
apă, putem prefera pe una dintre ele, după alte criterii. De exemplu, 
dacă F consumă mai multă apă decît C, vom prefera prima soluție, 


8 8 


è 
(O) 
© 


A F 
Fig. 1.3.8 


deoarece pentru arborele din figura 1.3.9, care este un drum hamil- 
tonian, F se găseşte la capătul reţelei, deci apa va avea în F o pre- 
siune mai mică decît în O. i 


Fig. 1.3.9 


Să arătăm acum că algoritmul lui Kruskal ne magie 
într-adevăr un arbore parţial de cost minim. Dacă graiul conex 
de la care plecăm are n virfuri, toţi arborii parţiali ai săi vor avea 
n — 1 muchii. i 
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„4 — 0.2161 


Să presupunem, prin reducere la absurd, că arborele parțial A. 
obţinut cu algoritmul lui Kruskal nu este minim, ci există un 
alt arbore A, de lungime totală minimă, deci 1(4,) < UA). Dacă 
muchiile lui A sînt Uy, a.. -șta_a Să presupunem că numai primele 
It muchii sînt comune cu arborele A, deci arborele A, are muchiile 
Uyy Map: -p Ury Die eo Daca UNDAE Urti. y Un Sînt muchii dife- 
rite de Vry Urto- -s Un-ı (putem avea întotdeauna o astfel de 
situație, eventual prin renumerotarea muchiilor). 

Să adăugăm muchia Ur, = [2 2] la arborele A, Deoarece 
tr, este o muchie care nu aparţine arborelui A, rezultă că în 
graful astfel format G, apare un cirenit unie, conform unei proprie- 
tăţi a arborilor. Acest circuit conține măcar o muchie v(i > k +1) 
care nu aparține arborelui A, deoarece în caz contrar ar rezulta că 
arborele A conţine un circuit, ceea ce contrazice definiţia unui arbore. 
Dacă suprimăm acum muchia v; obţinem tot un arbore A». În- 
tr-adevăr, graful obţinut din G, este tot conex şi nu conţine cireuite, 
deoarece am suprimat circuitul unice al grafului G}. Lungimea; acestui 
arbore este egală cu : 


UA) == UA) + Htr) — Ho). 


Însă, din descrierea algoritmului lui Kruskal, ttz, este cea mai 
scurtă muchie care nu formează circuite cu muchiile tyy... Ur 

Cum nici muchia v; nu formează circuite cu aceste muchii, 
deoarece A, este arbore, rezultă că avemi următoarea inegalitate 
între lungimile celor două muchii: Hur) SHu), deci (Aa) < 
< (4). Continuînd acest procedeu înlocuim toate muchiile Vr+1,- . 
«+ 39a-a ale arborelui A, cu muchiile p+r. -ttn ale arborelui A, 
obţinînd un şir de arbori: - 


Aa As. x «Ån =A 


astfel încît lungimile lor să verifice inegalitățile : ; 
UA, > (Aa) > UA) > ... > UAn) = UA) deci (4) > UA). Am 
ajuns astfel la o contradicție, deoarece am presupus că (4A) < HA) şi 
deci algoritmul lui Kruskal este complet justificat. 
vederea programării algoritmului la un calculator, este bine 

ca operațiile cerute să fie aranjate într-un anumit mod, în special 
pentru a verifica condiţia ca muchia nou aleasă să nu formeze 
circuite cu muchiile deja alese. i - 

Se observă din descrierea algoritmului că grafurile parțiale 
care se formează sînt; toate arbori, care la sfirșit se leagă între ele 
într-un singur arbore parțial care uneşte toate viriurile graiului. 
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Din acest motiv este bine ca componentele conexe (subgraiturile 
conexe maximale) obţinute pe parcurs să fie reprezentate prin 
indicii virfurilor din acea componentă, fiecare componeată fiind 
scrisă pe câte un nivel. Dacă alegem o muchie de lungime minimă 
dintre cele nealese, pot apărea următoarele situații : 

a) Ambele ei extremităţi se găsesc pe un nivel, adică fac parte 
dintr-o componentă conexă a graului obținut pînă în acel moment 
(componentă care este un arbore). 

În acest caz, legînd două viriuri ale arborelui printr-o nouă 


. muchie, se formează un circuit, deci vom alege o altă muchie în 


ordinea crescătoare a lungimilor, 

b) O extremitate a muchiei se găsește pe un nivel, cealaltă 
extremitate se găseşte pe un alt nivel. 

Alegerea, este bună, cei doi arbori se unesc într-un singur arbore, 
deci vom unifica cele două nivele pe care se găsesc cele două extre- 
mităţi ale muchiei. | 

c) O singură extremitate a mmchiei se găsește pe un nivel, iar 
cealaltă nu apare pe nici un nivel. i 

În acest caz alegerea este bună, vom adăuga la acel nivel şi 
extremitatea a doua a muchiei selectate. i 

d) Nici o extremitate a muchiei alese nu. se află pe nici un 
nivel. Alegerea este bună, vom adăuga un nou nivel, format; din 
cele” două extremităţi ale muchiei selectate. i 

La sfîrşitul aplicării algoritmului vom avea un sińgur nivel, 
graful parţial obţinut va fi un graf conex (avind o singură compo- 
nentă eonexă) și fără circuite, deci va fi un arbore. Deci sfîrşitul 
aplicării algoritmului va fi indicat de obţinerea unui singur -nivel 
care conţine toate vîrfurile graiului. 

Muchiile alese-le vom marea, într-un fel oarecare, de exemplu 
prin subliniere. 7 

Reluînd exemplul grafului din figura 1.3.6, vom avea urma- 
toarea listă de muchii, care au ataşat în paranteză costul : 


IB, 8] (3); [4, E] (2; [0,_D] (2); [B,_D] (3); [B,_F] (3); 


(6, PA; (4 815; B, 0l (9; 5, 0 (8; D, 7) (8; 
[4, B] (5); LE, F] 6); [S, B] (7). S 
Această listă reprezintă datele problemei. 
Aplicînd algoritmul descris cu ajutorul nivelelor, obținem : 
1) Nivelul 1: {B, E) 
2) Nivelul 1: (4, B, E) 
3) Nivelul 1: (4, B, B) 


DI 


Nivelul 2: (0, D} 

4) Muchia [B, D] produce unificarea celor două nivele: 

Nivelul 1: (4, B, C, D, E 

5) Nivelul 1: {4, B, 0, D, E, F} 

6) Muchia [C, F] nu poate fi aleasă, deoarece extremităţile 
sale aparţin nivelului 1. Este selecționată următoarea muchie (în 
ordinea crescătoare a lungimilor), (A, 8]: 

Nivelul 1: (4, B, 0, D, E, F, S) 

Deoarece s-a obținut un singur nivel, care conţine toate virturile 
gratului, aplicarea algoritmului este încheiată. i 
i Muchiile care formează arborele minim (desenat în fig. 1.3.8) 
sînt muchiile subliniate în lista, iniţială, pentru care suma lungimilor 
este egală cu 14. 

În momentul selectării unei muchii, ea a fost subliniată, în lista 
de muchii care constituie datele problemei. 

$ În acest fel pot fi tratate probleme oricît de mari, toate opera- 
pile reducîndu-se la crearea, sau unificarea unor nivele, în funcţie 
A roenan căutării în lista de muchii ale grafului, după regu- 

ş „Pentru a uşura, căutarea în lista de muchii, este bine ca acestea 
să fie scrise în listă în ordinea crescătoare a lungimilor. l 

i În continuare vom reveni asupra problemei drumului hamilto- 
nian minim într-un graf şi vom arăta legătura acestei probleme cu 
problema arborelui minim. 

„___ Fiind dat un graf neorientat G cu n viriuri, care are asociată 
fiecărei muchii fu, v] o lungime l[u, v] > 0, vom determina cu 
algoritmul lui Kruskal un arbore parţial de lungime minimă. Dacă 
acest arbore are numai vîrfuri de gradul unu şi de gradul doi, atunci 
el este un drum hamiltonian, deci este chiar drumul hamiltonian 
de lungime minimă, deoarece mulţimea drumurilor hamiltoniene ale 
grafului G este inclusă în mulțimea arborilor parţiali ai lui &. 

Să demonstrăm că dacă un arbore parțial al lui G are numai 


vîrturi de gradul unu şi de gradul doi, atunci el este un drum 


hamiltonian în graful G. 
ntr-adevăr, notind prin g numărul virfurilor arborelui de 
gradul unu, vor rezulta n — q vîrfuri de gradul doi. i 
Cum suma gradelor este egală cu dublul numărului de muchii, 
rezultă că : 


q+2(n—g9)=2(m—1), 


deoarece un arbore cu n vîirfuri are n — 1 muchii. Din această 
ecuație se deduce g = 2, deci arborele are exact două virturi de 
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a ner onm De Tarcea 


- gradul unu, iar n — 2 vîrfuri de gradul doi, deci el este un drum 


hamiltonian. 
Dacă arborele parțial minim nu este un drum hamiltonian, vom 


descrie un algoritm euristie de obţinere a unui drum hamiltonian 
minim pornind de la'arborele parțial minim obținut cu algoritmul lui 
Kruskal. Să lămurini mai întâi ce este un algoritm euristic. Există 
multe probleme de optimizare pentru care nu cunoaștem un algoritm 
de rezolvare, sau dacă un astfel de algoritm poate fi conceput, re- 
alizarea, lui la calculatorul electronic necesită un volum imens de 
calcule, care fac ca, aplicarea, acelui algoritm să fie practic imposibilă. 

De aceea, uneori sînt preferați algoritmii numiţi „euristicitt, 
care nu produc în mod sigur soluția cea mai bună din punct de 
vedere matematic, ci o soluţie care este destul de bună pentru 
practică într-un timp şi cu un volum de calcule rezonabil. De obicei 
aceşti algoritmi pleacă de la o soluţie posibilă a problemei (nu 
neapărat optimă) şi o îmbunătăţesc succesiv prin aplicarea unuia 
sau mai multor criterii, specifice problemei respective. 

Cînd procesul de îmbunătăţire a soluţiei nu mai poate avea 
loc, aplicarea algoritmului se opreşte. Nimeni nu ne asigură însă că 
dacă plecăm de la o altă soluţie inițială sau dacă aplicăm un alt 
criteriu de optimizare, soluţia obţinută pe această cale nu este 
mai bună decit soluţia obţinută anterior. De multe ori în practică 
rapiditatea luării unei decizii este mai importantă decît micul cîştig 
pe care îl poate aduce soluţia, cea mai bună față de o soluţie destui 
de bună, dacă, soluţia, cea mai bună se obţine într-un timp mult 
ma îndelungat şi cu capacităţi de calcul mult mai puternice decit 
soluția euristică, i 

De aceea în Cercetarea Operajională întîlnim. la ora actuală, 
alături de algoritmi exacţi, mulţi algoritmi euristici sau aproximativi, 
care sint foarte utili în abordarea problemelor complexe pentru 
care găsirea unei soluţii analitice exacte este sau imposibilă sau 
foarte greoaie. ă 

Algoritmul euristic pe care îl vom prezenta pleacă de la arborele 

arţial minim al unui graf G, transformîndu-l într-un drum hamil- 
tonian al lui 4 (despre care nu putem îi însă siguri că este un drum 
hamiltonian minim), cu minimizarea la fiecare etapă a creșterii 
lungimii totale. i 

Pentru justificarea algoritmului, să observăm mai întîi că un 
arbore parţial al unui graf G care nu este dum hamiltonian se poate 
transforma, într-un arbore parţial care să conțină virfuri cu grade 
mai mici, prin următoarele operații: 

01) Se unese două vîrfuri terminale e şi y printr-o muchie [z, yl. 
În acèst. mod apare un circuit unic în graf. Pe acest circuit va exista 
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cel puţin un virf z de grad d(2)>3, deoarece în caz contrar arborele 
de la care am plecat era un drum hamiltonian şi problema era rezolvată. 

g Vom suprima din acest circuit o muchie [v, z] adiacentă cu un 
virf a de grad mai mare ca doi, graful obținut fiind tot arbore şi 
avind un vîrf sau două virfuri de grad mai mare ca doi pentru care 
gradul s-a micşorat cu cîte o unitate, fără să mai apară alte viriuri 
de grad mai mare ca doi (fig. 1.3.10). . 
02) Un vîrf terminal y, legat de un vîrf œ cu d(z) > 3 se intercalează 
între alte două virturi, z și t, legate printr-o muchie, astfel încît 
gradele virfurilor z şi t să fie cel mult doi. 

Pentru aceasta, se suprimă muchiile [, y] și [2, t] şi se adaugă 

muchiile [y, 2] şi fy, î] (fig. I.3.11). 


Fig. 1.8.10 Fig. 13.11 


93) Un vîrf neterminal, v, legat de un virf œ cu d(z) > 3, se leagă 
de un vîrf terminal y care nu face parte din aceeași componentă 
conexă cu v, obţinută prin suprimarea muchiei [v, 2].. Apoi se 
suprimă muchia [&, v], obţinîndu-se tot un arbore (fig. [.3.12). 


Fig. 1.3.12 


În primul caz creşterea costului arborelui este egală cu : [4] — 
— liv, 2], în al doilea caz ea este egală cu : liy, 2] + liy, t] — îl, y] — 
— ilz, î], iar în al treilea caz ea este egală cu: liv, y] — lis; o]. 
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Să observăm că în toate cele trei cazuri printre virfurile de grad 
mai mare sau egal cu trei ale arborelui, cel puţin unul își micşorează 
gradul cu o unitate, iar nici un astfel de virf nu își mărește gradul. 

Evident că după aplicarea repetată a acestui procedeu pentru 
un graf complet; vom obține un arbore care conține numai viriuri 
de gradul unu sau doi, deci care este un drum hamiltonian. 

Pentru a avea o lungime cît mai mică a drumului hamiltonian, 
este bine ca transformările succesive făcute asupra arborelui să 
minimizeze, în fiecare etapă, creşterea costului (lungimii) arborelui. 

Pentru aceasta vom proceda astfel : 

1) Vom alcătui lista vîrfurilor terminale Yı, Yz- -Yy Ale arbo- 
relui, adică a viîrturilor cu gradul unu. 

Pentru fiecare virf y; vom considera pe rind toate virfurile 
y ca j>i 

Pentru un j ales, vom determina cea măi lungă muchie de pe 
drumul care îl uneşte pe y; cu Yp adiacentă cu un virf de grad 
mai mare ca doi, apoi vom determina, creșterea lungimii corespun- 
zătoare eliminării acestei muchii cu operația ol. Pentru fiecare virt 
terminal vom calcula apoi creşterea minim a lungimilor, corespun- 
zătoare legării sale cu celelalte vîrfuri terminale, de indice mai mare. 

2) Pentru toate vîrturile terminale y legate de vîrfuri de grad 
mai mare ca doi, vom determina minimul creșterii lungimilor prin 
operația de intercalare, descrisă de 02. 

3) Pentru toate virturile neterminale v, legate de virfuri œ de 
grad mai mare ca doi, vom determina minimul creşterii lungimilor 
prin legarea lor cu virfuri terminale y (care nu aparţin aceleiași 
componente cu v prin suprimarea viriului œ), prin operația 03. 

4) Dintre minimele astfel calculate la 1), 2) şi 3), vom, alege 
cea, mai mică valoare. Dacă există mai multe vîrfuri care realizează 
aceeaşi valoare minimă, vom alege acea transformare care scade 
gradul cu cîte o unitate la două virturi de grad mai mare ca doi 
(numărul maxim). i 

Aplicăm transformarea care realizează minimul creşterii costului 
arborelui, de tipul o1, sau 02, sau 03. 

. Dacă arborele obţinut are numai vîrfuri de gradul unu sau doi, 
adică este un drum hamiltonian, ne' oprim, Dacă nu, se reia algo- 
ritmul de la punctul 1). i 
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et Sapca acest algoritm euristic pentru graful complet G- eu 
virfuri 25 Zo, ay Ca sp Lg, avînd următoarea matrice a lungimilor 
muchiilor : 


Xi - ` Ta ta Ta Ts Tg 
o j u | 27 15 2 | 3 
14 |. o | 18 | 3 6 | 2 
C E E | o j 35 17 | 18 
15 ! 3 j 5 o |] 1 | 
20 |e | 17 1 | o j 
15 | 2 | 18 4 | 4 - 


Aplicind algoritmul lui Kruskal, găsim arborele minim din figura 
I.3.13, pentru care vîrturile terminale sînt £j, La; La, g 


LA Xa Xy 


% 


Xp Si Xs 


Fig. 1.3.13 


Prin operaţia ol, ealculăm astfel minimul creşterii lungimii î 
Pentru s: se adaugă muchia [æ 8] şi se elimină muchia 
[2i 42]. Creşterea lungimii : 1. 


MU [%3, 24]. Creşterea lungimii : 2. 
| | Pentru &,: se adaugă muchia [5 &] și se elimină muchia 


| [2 a]. Creșterea lungimii: 1. Operația 02 nu se poate aplica, 
Mi deoarece nu există muchii cu ambele extremităţi de gradul unu 
| Ru da Pentri operaţia 03 trebuie să cercetăm numai virfurile 
| 2 Și a. 
Il Pentru w, minimul creşterilor se realizează adăugînd muchia 
| (2, 4] şi eliminînd muchia (2, 24], cu creşterea lungimii : 3 unități. 
(i Pentru 4, minimul creșterilor se realizează adăugîind muchia 
| [24 2] şi eliminind muchia [Se 24], cu creşterea lungimii de o 
unitate, Vedem deci că în trei cazuri creşterea, este minimă şi anume 
l egală cu o unitate. i 
INI În două cazuri din trei cite un singur virf de grad mai mare 
ca doi şi-a micşorat gradul cu o unitate, iar în cel de-al treilea caz 
două virfuri și-au micșorat gradul cu cite o unitate. i 
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Pentru aş: se adaugă muchia [a %;] şi se elimină muchia 


Vom alege deci acest caz (adăugarea muchiei [£s æ] si elimi- 
narea muchiei (7, t4]), pentru care se obţine chiar drumul hamilto- 
nian minim din figura 1.3.14, de lungime egală cu M4+2+4+ 
+14 15 = 36 unităţi. 


>, -X Xe Xs a 3 


Fig. 1.3.14 


Pentru cercetarea directă a tuturor cazurilor posibile (în număr 
de 61/2 = 360) trebuia să efectuäm 1440 adunări şi 359 comparații. 

Este clar că drumul hamiltonian obținut este minim, deoarece 
arborele minim (care este unic) are o lungime de 35 unități gi nu 
este drum hamiltonian, deci orice drum hamiltonian are o lungime 
de cel puţin 36 unități, lungimile muchiilor fiind numere întregi. 

Se pot construi şi exemple în care acest algoritm. nu ne conduce 
la drumul hamiltonian minim, avantajul său constind în simpli- 
tatea aplicării din cauza numărului redus de operații. 

Desigur că aici nu am făcut decît o scurtă introducere în teoria 
grafurilor, am prezentat unele noţiuni și concepte, algoritmi de 
calcul şi aplicaţii. 

Pentru cititorul dornic să cunoască mai multe despre teoria 
grafurilor şi aplicaţiile ei, recomandăm următoarele cărți apărute 
în limba română : : 

1. 0. Ore, Grafele și aplicațiile lor, Editura ştiinţifică, Bucureşti, 
1968. 

2. O. Berge, Teoria grafurilor şi aplicaţiile ei, Editura tehnică, 
Bucureşti, 1969. a 

3. H. Ionescu, V. Burlacu, 0. Dinescu, Teoria grafelor cu unele 
aplicaţii în economie, Editura ştiinţitică, Bucureşti, 1969. 

4. Al. Roșu, Teoria grafelor, Editura militară, Bucureşti, 1969. 

5. I. Aurian, Gh. Boldur, S. Lazăr, Cercetarea operaţională în 
constructii, Editura ştiinţifică, Bucureşti, 1967. 

6. A. Kaufmann, G. Desbazeille, Metoda drumului critic, Editura 
tehnică, București, 1971: 

7: I. Tomescu, Introducere în combinatorică, Editura tehnică, 
Bucureşti, 1972. i 

9. T. Ionescu, Grafuri, Aplicaţii. vol. I—II, Editura didactică 
şi pedagogică, București, 1974. d 
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|.4. Probleme 


1. Să se arate că cele două desene din figura 1.4.1 reprezintă, 
unul şi același graf. Aceeaşi problemă pentru figura 1.4.2. 
Două astfel de grafuri se numesc grafuri izomorfe. 


ALO D 


Fig. LA Fig. 1.4.2 
R: Se vor numerota convenabil virturile. 
2. Să se arate că pentru un graf G numărul virfurilor sale de 
grad impar este par (Ind: Se va utiliza teorema 1.1.1.). 


3. Să se arate că într-un graf fără circuite nu pot exista drumuri 
neelementare. 


4. Să se determine drumurile elementare dintre virfurile 1 şi 
7 ale graiului din figura 1.4.3. 


2 4 6 
1 

3 5 7 

Fig. 1.4.3 


R: 13 drumuri elementare, 7 care incep cu muchia [1, 2] 
și 6 care încep cu muchia [1, 3]. 


5. Să se găsească un drum şi un circuit hamiltonian pentru 
graful din figura 1.4.3. 

6. Să se arate că numărul circuitelor hamiltoniene ale graului 

5 i PERE n —1)! - 
neorientat şi complet cu n > 3 vîrfuri este egal cu a - (Ind: 


Se va folosi inducția completă după n). 
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7. Pentru graful din figura 1.4.4 să se găsească drumul hamil- 
tonian de timp minim şi drumul hamiltonian de timp maxim. 


Fig. 1.4.4 


R: Prin enumerare directă se găseşte că din cele 24 permutări ale 
virturilor 1, 2, 3, 4 numai 14 corespund unor drumuri hamiltoniene in 
graiul dat. Drumurile hamiltoniene minime sint (3, 2, 4, 1) şi (4, 1, 3, 2) 
de timp 4, iar drumul hamiltonian de timp maxim, egal cu 10, este 
(Q, 1,4,3). 


8. Un delegat al unei instituţii din oraşul 1 trebuie să meargă 
în orașele 2, 3, 4, 5, apoi să se întoarcă în oraşul 1. 


Fig. 1.4.5 


Cunoscînd costurile deplasării dintr-un oraș în altul, care sint 
numerele asociate muchiilor grafului din figura 1.4.5, să se determine 
succesiunea de oraşe pentru care costul total al deplasării să fie 
minim şi ştiind că sosirea în oraşul 4 este mai urgentă decît sosirea 
în oraşul 2. 


R: Dintre cele 8 circuite hamiltoniene ale graiului dat există 2 cu 
cost minim (egal cu 12 unităţi) şi anume : [£, 5, 2, 3, 4,1] şi [1,4,3, 2, 5,1]. 
Răspunde cerinţelor problemei circuitul (1, 4, 3, 2, 5, 1). 


9. Să se găsească drumul critic (de timp maxim) dintre vîrturile 
A şi G ale grafului din figura 1.4.6, precum şi intervalele de fiuc- 
tuaţie ale evenimentelor şi marginile totale ale operaţiilor. 
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10. Aceeaşi problemă pentru graful din figura 1.4.7 și perechea 
de noduri A şi G. 


a 
o 
X 
3 
14. Folosind algoritmul lui Dantzig, D1, să se găsească distan- 


tele minime (în km) pe calea ferată de la Bucureşti la celelalte 
oraşe ale ţării, care apar ca vîrfuri ale graului din figura 1.4.8. 


Suceava 
sesti p 


7 


Fig. 1.4.6 


ar 


Ma 


Distanţele directe în km, asociate fiecărei muchii, sînt scrise lingă h. £ 
muchia respectivă. i S 3 
T S 
12. Pentru graful orientat din figura I.4.9, să se calculeze distan- Să 
tele minime de la virful A Ja celelalte virturi, folosind algoritmul Dl 
al lui Dantzig, precum şi drumurile minime care pleacă din A. 


Baia Mare 
Roşiori 


13. Să se arate că orice graf-turneu cu n > 4 virfuri conține > 
un subgraf cu 3 vîrfuri de forma următoare indicată în figura 1.4.10 
Und: Este suficient să se demonstreze pentru n = 4, analizind toate 
situaţiile posibile). 


„Satu Mare 


14. Pentru un graf-turneu cu n viriuri Sy, Vz.. - Bu SE notează 
cu 7, numărul arcelor care intră în æ; și cu s numărul arcelor care 
ses din a. Să se arate că: 


Craiova 


I)rxrs=n-—l 


Caransebeș $ 


e 
îm 
în. 
| 
iM 
M 


D 
L 
N 
Q 
e 
Š 
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Fig. I.4,8 


Inä: Relațiile 1) şi 2) se deduc pornind de la definiția unui graf- 
turneu, tinind seama că numărul muchiilor grafului complet cu n 
vârfuri este egal cu C2, iar 3) se deduce din 1) prin înmulţire cu F; 
respectiv s; şi prin însumare, ținînd seama de 2). 


Fig. 1.4.10 


15. Dacă numărul m al muchiilor unui graf G cu n virfwri 
verifică inegalitatea : 


m > Oaa 


să se arate că graful G nu are virfuri izolate. 
Ind: Se va arăta că dacă G are un virf izolat, atunci m <Q. 


16. Pentru un graf conex cu m muchii şi n virfuri, să se arate 
că numărul minim de muchii care trebuie înlăturate astfel încit 
graful obținut să nu conţină circuite este egal cu m — n + 1 (acest 
număr se numeşte numărul ciclomatic al gratului). 


Ind: Se va ţine seama de teorema I.3.1. 


17. Să se arate că un arbore cu n > 2 virfuri are cel puţin 


două virfuri de gradul 1. | 
Ind: Se va ţine seama de teorema I.3.1. şi de teorema I.1.1. 


18. Să se demonstreze prin inducţie teorema 1.3.1. 


19. Să se găsească arborele de cost minim pentru graful conex 
din figura I.4.11 şi pentru graful conex din figura 1.4.12. 


20. Un oraş este așezat pe malurile unui rîu, peste care există 
cinci poduri, aşezate ca în figura alăturată (fig. 14.13). 

Este posibil, ca plecind dintr-un anumit; loc, să trecem o singură 
dată peste cele cinci poduri, întorcîndu-ne de unde am plecat? 


$2 


„__ Rs Considerind multigraful din figura 1.4.14, unde virfurile Tepre- 
zintă Tespectiv malul drept al rîului, malul sting al rîului şi insula, iar 
muchiile reprezintă legăturile acestora prin poduri, Euler a arătat că o 
problemă de acest fel nu are soluţie, deoarece orice trecere printr-un virf 


Fig. LA Fig. 1.4.12 
MD 
Fig. 1.4.13 Fig. [4.14 


consumă două muchii, pe care nu le mai putem utiliza, iar multigratul are 
virîurile MD şi MS de grad impar, egal cu trei. Deci nu ne putem deplasa - 
astfel încît să parcurgem toate muchiile multigrafului din figura 1.4.14 © 
singură dată. 


nai 21. Să se coloreze regiunile din figura 1.4.15 cu 4 culori, astfel 
încît fiecare două regiuni care au o frontieră comună să fie colorate 
cu culori diferite. 


"Fig: 1.4.15" 
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Este posibil să realizăm.-o astfel de colorare cu numai trei culori? 


R: Nu este posibil, deoarece A, B, C, D au frontiere comune două 


cite două. 
(S-a demonstrat că pentru o hartă de acest fel cu cel mult 40 regiuni 


sint suficiente pentru colorare cel mult 4 culori. Problema colorării hărților 
cu cel mult 4 culori rezistă de aproape un secol eforturilor matematicie- 
nilor de a o rezolva în cazul general.) 


29. Să se deseneze graful relaţiei de divizibilitate avînd ca 
virfuri numerele 2, 3, 4, 5, 12, 15, 18, 20. 


23. Să se arate că numărul circuitelor elementare ale 
graului neorientat și complet cu 2 virfuri (n > 3) este egal cu 


1 2 nn — (n k ti), 


PX 


2 rs k 


24. Folosind algoritmul euristic, să se determine un drum 
hamiltonian de lungime minimă în graful neorientat, cu opt virfuri, 
care are asociată următoarea matrice a costurilor : 


1 2 3 4 5 6 7 8 
1 0 2 7 | 1 |: 10 8 9 21 
2 2 o 23 a5 | 35 14 | 7 5 
3 7 23 0 _8 1i it |243 22 
4115 45 8 0 12 | 32 a | 5 
5 | 10 | 35 11 12 o j 4&2 | 38 19 
6| 8 14 17 32 | 42 o j 3 36 
7 9 7 13 21 j 38 | 3 o | 18 
sian |53| 2 15 19 36 18 0 


R: Se obţine chiar drumul hamiltonian minim (3, 5, 4, 8, 7,2, 6, 1] 
de lungime 85. 


25. Să se demonstreze că problemele Pi, P2 și P3 de la drumuri 
şi circuite hamiltoniene sînt echivalente şi eu problema P4: 
P4) Găsirea unui drum hamiltonian minim avind una dintre 
extremităţi într-un virf dat al grafului. 
R: Se adaugă grafului dat un nou virf, legat de celelalte viriuri prin 
muchii de lungime nulă. 
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E =7); 


26. Să se construiască graful programului, să i 
drumul critic, intervalele de fuctuațis le ernim loa pi arnie 
totale ale operațiilor în următoarea situație de construcție a aii 
ansamblu hidroelectric, cu unitatea de timp luna : Operațiile care 
trebuie să fie realizate sînt următoarele : 


a) construcţia drumurilor de acces la uzină şi 1 i 
obținerea materialelor (durata t, = 4); ER piru 


b) pregătirea carierelor de exploatare şi a fundațiilor (t, = 6); 


& Bica taiata unui centru pentru personal şi administraţie 
AA ; i 


d) comanda şi construcţia materialului ie şi hi i 
5 ţi ui electrice hidr: 
(generatoare, turbine, conductă forțată etc.) (t4= 12); a 
e) construcția uzinei (t = 10); i 


construcții iului isuri : ; 
tai A DU. ţia barajului, a digurilor şi a deversorului de supra- 


g) construcţia, galeriilor de fugă şi a conductelor de aducțiune 


h) montajul uzinei şi al conductelor (t, = 10); 

i) verificarea funcţionării (t = 3). 

Evenimentele sau etapele programului sînt; următoarele : 
1) pornirea lucrărilor ; 

2) terminarea drumurilor đe acces ; 


3) terminar i ; PI sea ăi 
PR oaia rea centrului pentru personal şi administrație şi a 


Fig. 1.4.16 


5 — e. 2161 


4) terminarea uzinei şi æ conductelor de aducțiune ; 
5) terminarea ansamblului energetic ; 
6) livrarea, lucrărilor către beneficiar. 


R: Se obţine graful din figura 1.4.16. 


27. Să se arate că pentru orice graf-iurneu G există un virf æ 
astfel încât toate virturile y # æ ale grafului pot fi atinse plecind 
din æ pe drumuri care au un arc sau două arce. 


R: Se va lua pentru æ virful din care pleacă un număr maxim de 
arce în G. Dacă un virf z nu poate fi atins din x pe drumuri de lungime 
unu sau doi, se deduce că există arcul (z, <£), iar din z pleacă arce către 
toate virfurile la care sosesc arce din z, ceea ce contrazice definiţia lui x. 


28. Să se calculeze numărul grafurilor cu n virfuri date a 
2 ** "0 Și fără bucle, avind proprietăţile următoare : 
a) neorientate ; 
b) orientate ; 
c) grafuri-turneu ; 
d) orientate şi complete. 
R: a) 2%; b) 42; c) 2%; d) 3%, unde «= GG. 


29. Să se arate că un graf neorientat cu n virfuri și n muchii 
conţine cel puțin un circuit. 
R: Dacă graful nu conţine circuite, fiecare componentă conexă a sa 
este un arbore cu g viriuri și g — 1 muchii. 


30. Să- se demonstreze că un graf neorientat cu n virfuri şi 
m muchii, m > C, + 2, conține cel puţin un circuit hamiltonian. 
R: Se va arăta că dacă G nu este graful complet Fa» atunci există 
un virf æ cu d(x) = n — 2, apoi se va utiliza inducția după n, considerind 
subgraful obţinut din G prin suprimarea lui g. 


II. Noţiuni de programare liniară 


1.1. Formularea problemelor de programare liniară 


Pentru a vedea ce sînt problemele de pro rana iniană că gX 
citeva exemple concrete de astfel de Aroblenes, aie limata Bidai 


1. Un prim exemplu este pri i i 
! X problema transportului unui prod 
de la anumite centre de produeţi i rare 
ei e Aole producție la anumite puncte de desfacere 
Să considerăm că întreprinderea 7, î ămi 
onsiderăm ı produce într-o săptămină 
T ooo t cărămizi, iar I, produce în acelaşi interval de timp 650 t. 
ceste cărămizi trebuie transportate la consumatorii C, şi 0,. Se 
ştie că 0, poate primi cel mult 700 t și 0z poate primi cel mult 


-950 t, iar costul transportului unei tone de produs este dată de 


mairicea. următoare : 


Cc, GA 
L 40 leijt 60 leijt 
I, | 50 lejt 70 leijt 
Disponibil Capacitate ge 
40 lei fe depozitare 
1000. L, / G 700 


Fig. 1.1.1. 
Interpretarea acestei matrice a costurilor este următoarea : 


La intersecția liniei ¿ cu coloana j se găsește costul (în lei/tonă) 
al transportului unei tone de cărămizi de la întreprinderea cu numărul ¿ 
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la consumatorul cu numărul j. Astfel transportul unei tone de la 
Ila 0, costă 40 lei, transportul unei tone de la 4 la 0» costă 60 lei ete. 

Problema care se pune este aceea de a determina un plan de 
transport care să ţină seama de cerințele problemei şi care să ducă 
la cheltuieli de transport minime. 

Din punct de vedere matematic, aceasta, revine la determinarea, 
numerelor v, y, Z, t, care au următoarea semnificație: œ reprezintă 
numărul de tone transportate de la Z, la Ou, y reprezintă numărul 
de tone transportate de la I, la 02, 2 reprezintă numărul de tone 
transportate de la T, la O, şi i reprezintă numărul de tone transportate 
de la I, la C, Ținind seama de această, interpretare, obținem că 
numerele æ, y, 2, t trebuie să fie nenegative, adică g, y, Z t20. 


Restricţiile problemei se scriu astfel : 
1) Cantitatea transportată din I, trebuie să fie egală cu dispo- 
nibilul din J, adică: 
æ+ y= 1000 
2) Cantitatea transportată din I, trebuie să fie egală cu disponi- 
bilal din 7, adică: 
z + t= 650 
3) Cantitatea transportată în 0, trebuie să nu depăşească capa- 
citatea de depozitare în O, adică : 
A a +z < 100 
4) Cantitatea transportată în 02 trebuie să nu depăşească, capa- 
citatea de depozitare în Oz adică : , 
y + 6 950 
5) Costul total al transportului, adică, : 
40% + 60y + 502 + 70t 
să fie minim; adică să aibä cea mai mică valoare posibilă, ceea ce 
vom nota simbolic : 
min (40 + 60y +- 502 + 701). 


'Prebuie de asemenea să fim atenţi dacă problema pusă este posibilă, 
adică poate avea o soluție. Într-adevăr, adunind cele două egalităţi 
de la punctele 1) şi 2) obţinem : 
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Prin adunarea membru cu membru a celor două i ități 
=: ereziei r două inegalităţi de la 


æ +y +e+t s 1650 


Aceste două relații nu vin în contradicție şi di 
u: t mt şi din compararea lor 

deducem că problema este compatibilă, adică kapaditates depozitelor 
este Aa puțin egală cu producția din J, şi Ip. 

acă de exemplu C, putea primi 600 t iar C, pute: imi Í 

NE caiac 1 Oti putea primi 900 i 
pro aie a n opani formulare era imposibilă, deoarece pe de o parte 
æ +y -+z + t= 1650, 


iar pe de altă parte am fi avut: 


æ +y +z- t= (s +g) + (y + î) s 600 + 900 = 1500, 


ceea ce era imposibil. 
Formularea matematică a i 
j a: acestei probleme d; 
deci următoarea, : ii e 
Să se determine numerele 
se € reale v > 0 
care verifică relațiile : pândă da d a dl 


æ + y = 1000 
z -+= 650 
æ+z< 700 


y+i< 950 
astfel încât : 


min (40 æ + 60 y + 50 z +70 ŝi). 


Uneori problemele de transport se scriu numai cu egalități ji 
în loc de @ + 2 < 700 şi y +t < 950 vom serie: 2 i ea 
y + t = 950, dacă interpretăm problema astfel: 0, este un consu: 
mator care are nevoie de 700 t cărămizi, iar O, este un consumator 
Ca ue deci t caramizi acest caz, pentru ca problema 
posibilă, trebuie ca di; ibi i i iar 
lut posibi pr e de la Z, şi I, să fie chiar egal 
Problema de transport enunțată se poate formul 
general (însă tot pentru două întreprinderi ti două ae ari af 
Se ştie că întreprinderea, J, produce într-o anumită perioadă 
de timp o cantitate de produse egală cu a, unităţi (a, > 0), între- 
prinderea I, produce în aceeaşi perioadă o cantitate de produse 
egală cu a2 unităţi (a. > 0), depozitul C, poate primi cel mult e, 


69 


unităţi (e, > 0), depozitul 0, poate primi cel mult c, unităţi (ea >0), 
iar ta Al aan pe “unitate de produs sînt date de matricea : 


G Ca 
h a kaet] 
h or Qaz 


Se pune problema determinării unui plan de transport optim, 
adică a unor numere æ, y, z, t care verifică următoarele restricţii : 


2>0,y>0,2>0,420 


æ + y= 
z+ t= m 
a+Hz<s& 
y+is ce 


min (a8 -+ aY + taZ F teat) 
Condiţia de compatibilitate se deduce ca mai înainte, scriind : 
a +a =g ty Hett lete) lyti S aHa. 


deci numerele a; de; Cy 62 trebuie să verifice inegalitatea : 


a F ag S G F Ca 


2. Un alt exemplu de problemă de programare liniară este problema 
dietei (a meniului optim). n i N Sc 
Sal ştie că viaţa unui om este condiționată, de anumite cantități 
minime de principii nutritive: glucide, proteine, grăsimi, albumine etc. 
Acestea se află în proporţii diferite în alimentele pe care le 
consumă, E ee 
Totalul acestor cantități este calculat pe baze ştiinţifice pentru 
un organism, sănătos în funcţie de sex şi vîrstă, ţinind seama și de 
eforturile la care este supus omul în cursul unei zile. ; 
Cantităţile minime din diferitele principii nutritive pe care le 
sumă un om reprezintă necesarul biologie. — i 
Sa Problema dietai cere să se determine cantităţile care trebuie 
consumate din fiecare aliment astfel încît, ținînd seama de proporţiile 
în care se găsesc principiile nutritive în fiecare din ele, să se asigure 
necesarul biologio la un preţ de cost minim. 
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Pentru exemplificare, se ştie că pentru un organism normal în 
condiţii de efort fizic obişnuit sînt necesare 3000 de. calorii şi 100 g 
proteine pe zi. 

Se ştie că 500 g piine conţin 1000 calorii și 25 g proteine, iar 
500 g brînză conţin 2000 calorii şi 100 g proteine. 

Preţul a 500 g pîine este de 2 lei, iar preţul a 500 g brinză este 
de 7 lei. Se cere să se determine cantitatea de piine şi cantitatea de 
brînză care trebuie consumată de un om într-o zi (în cazul unui 
„meniu compus numai din aceste produse), pentru a asigura necesarul 
biologie (care este 3000 calorii și 100 g proteine) la un preţ de cost 
minim. Pentru a formula matematic problema, să considerăni 1000 
calorii şi 25 g proteine ca unităţi de măsură a acestor principii nutri- 
tive, iar unitatea de măsură pentru pîine și brînză să fie de 500 g 

Să notăm cu æ numărul de unităţi de piine şi cu y numărul 
unităţilor de brînză ce urmează a fi consumate. Modelul matematic 
al acestei probleme este următorul : ` 


220, y>0 


s + 2y>3 (1) 


s+ 4y 4t 
min (2% + ty). 


Inegalitatea @ + 2y > 3 exprimă condiţia ca numărul de calorii 
rezultate prin consumarea a w unități de piine şi a y unități de 
brînză să fie cel puțin egal cu 3 unități de măsură pentru calorii, 
adică cu 3 :1000 = 3000 de calorii. 

Înegalitatea v + 4y > 4 exprimă condiţia ca numărul de grame 
de proteine existente în æ unităţi de piine şi în y unităţi de brînză; 
să fie cel puţin egal cu 4 unități de măsură pentru proteine, adică 
cu 4:25 = 100 g. 

Funcţia 2% -+ 7y reprezintă tocmai costul total în lei al meniului, 
al cărui minim se caută. 

O formulare mai generală a problemei este următoarea : 

Să presupunem că dispunem de trei alimente As Aa şi A cu 
care vrem să preparăm o mîncare care să satisfacă anumite cerinţe 
în substanţe nutritive (albumine, grăsimi, glucide) ale organismului 
și care să fie cît mai ieftină, 

Să presupunem că 1 kg din A, conţine au g albumine, a £ 
grăsimi, a g glucide, 1 kg din A, conține ap g albumine, as, € 
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răsimi, a: lucide, iar 1 kg din A conţine da $ albumine, @s g 
E alia, ași 5 Sucide. Mai ştim că necesarul organismului la o masă 
este de b, g albumine, b g grăsimi, bs g glucide, iar costul unui 
kg din Aa, As A; este respectiv de Cr Cz, Ca lei. 

Problema se retranscrie matematic sub forma, : 
Să se determine vectorul (2, y, 2)e R° cu componente nenega- 
tive : v, y, z > 0, care satisface restricțiile : : 


Gat + tY + az > bi 
ant + Goo] + A? > ba (2) 
Cast F Ggatj ++ das Z ba 


astfel încât ce -+ Cay + Cz Să aibă cea mai mică valoare posibilă 
pentru 4, Y, 2 verificînd (2), ceea ce vom scrie : 


min (8 -+ CY -t 032). 


Interpretarea variabilelor z, y, 2 este următoarea: „dieta asociată 
soluţiei (e, y, 2) conține æ kg din Ayn Y kg din A. şi 2 kg din Aa. 
Functia liniară cw -+ cy + 6a reprezintă tocmai costul acestei diete, 
iar restricțiile problemei ne anioni satisfacerea nevoilor organismului 
î ine, grăsimi și glucide. Eet des 
T iri clar o D, y ai o semnificație pentru problema dietei 
numai dacă sint nenegative. Condițiile w > 0, y > 0,2 > 0, împreună 
cu inegalităţile (2) formează restricțiile acestei probleme de progra- 
mare liniară, iar funcţia de minimiza 0? + Cay + 092 Se numeşte 
funcția obiectiv. în cazul problemei transporturilor şi a problemei 
dietei presupunem că produsul de transportat, respectiv alimentele 
As, Ap As pot fi divizate oricît de mult, deoarece căutăm soluții 
exprimate în numere reale nenegative. _ A 

Această ipoteză este realistă în multe cazuri, de exemplu 
cînd transportăm ciment în vrac, grîu, produse petrolifere etc. In 
alte cazuri însă produsele transportate pot fi piese, mașini, panouri 
şi în acest; caz cerem ca, soluţia să aibă componente întregi și pozitive. 

Astfel de probleme de programare se numesc probleme de pro- 

re “în numere întregi. DAR 
Erang Un alt ezemplu fo problemă de programare liniară este 
următorul: Într-o fabrică funcţionează trei procese de producţie 
notate Pa, Pa Ps care consumă materiile prime M, şi M, Și anume 
în unitatea de timp P, consumă @n unităţi din M, și d unităţi 
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Gin M, P, consumă a» unităţi din M, ṣi ap unități din Ma Pg 
consumă aa unități din M, şi ass unităţi din M,. În depozit există 
lun stoc format din b, unități din M, şib, unităţi din M,. Mai presu- 
punem că procesul de producție P, produce într-o oră o valoare 
egală cu c lei, procesul P, produce într-o oră o valoare egală cu c2 
lei, iar procesul P, produce pe oră o valoare egală cu cs lei. Prin 
aceasta se presupune că procesul P, (i = 1, 2, 3) funeţionind ¢ ore 
produce o valoare egală cu te, lei (ipoteza de liniaritate a modelului 
de producţie). i 

Se cere să se determine un plan de producție optim, adică un 
vector v = (m, Y, 2) e R3 cu æ, y, 2 > 0 astfel încât æ, y, 2 să repre- 
zinte respectiv timpul de funcționare în ore a proceselor Pi, Pa Ps 
care satisface inegalitățile : 


Ant F Aay F du < bu (3) 
Gap F- az] + Oa S bo 


astfel încât max (GE + Cy + 0a2), adică valoarea nou creată în pro- 
cesul de producţie (care constă din cele trei subprocese Py, Pa, P3) 
să fie maximă, cu condiția ca materia primă consumată să nu depă- 
şească rezervele întreprinderii la un moment dat. 
Kourt istoric. 

Programarea liniară datează din anii 1939—1941 cînd au apărut 
primele lucrări ale lui Z. Kaniorovici (1939) şi F. Hitchcock (1941). 
Din punct de vedere istoric primele probleme rezolvate au fost cele 


. đe transporturi, aplicabile la organizarea optimă a transporturilor 


maritime, la necesităţile de aprovizionare a frontului sau la plani- 
ficarea misiunilor aviației de bombardament. Metoda simplex a 
fost creată în 1947 de G. Dantzig şi J. von Neumann. Ulterior progra- 
marea, liniară a cunoscut un mare avînt prin lucrările unor matema- 
ticieni și economiști ca, T. Koopmans, L. Ford, D. Fulkerson, 
A. Charnes, W. Cooper, R. Gomory, H. Kuhn ete., găsindu-şi un 
cîmp foarte larg de aplicații în economie. 

Necesităjile reale au condus la apariţia şi dezvoltarea, altor 
tipuri de programări, cum ar fi programarea pätratică, programarea 
convexă, programarea în numere întregi, programarea stohastică, 
programarea dinamică etc., toate aceste domenii fiind înglobate în 
termenul generic de programare matematică. În domeniul programării 
matematice o serie de contribuţii importante au fost aduse de mate- 
majicieni şi economiști români. 

Deci o problemă de programare liniară constă în găsirea unor 
numere reale care verifică anumite inegalităţi, cum sînt (1), (2) 
sau (3), numite restricțiile problemei și care în plus fac maximă (sau 
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minimă, în funcţie de problemă) o anumită funcție liniară, adică o 
sumă de variabile înmulţite cu coeficienți numere reale. 

Un vector v ale cărui componente satisfac restricțiile unei astfel 
de probleme de programare liniară, se numeşte soluție posibilă sau 
soluție admisă. Dacă în plus o soluţie admisă maximizează (sau 
minimizează, în funcţie de problemă) funcţia liniară asociată acelei 
probleme, ea se numește soluție sau soluție optimă a problemei de 
programare liniară. i 

Am văzut din aceste exemple că mulţimea soluţiilor posibile 
pentru o problemă de programare liniară, este caracterizată printr-un 
sistem de inegalităţi liniare. Vom studia aceste mulțimi de puncte 
M(x, y) din plan care satisfac restricţii de forma inegalităţilor liniare 
şi vom arăta că aceste mulțimi sînt mulțimi conveae. 


11.2. Mulţimi convexe determinate 
de inegalităţi liniare 


În planul Oy orice punct M are două coordonate, abscisa æ l 


şi ordonata y. O ecuație de forma : 


az + by = c (1) 


reprezintă o dreaptă, ale cărei intersecții cu axa absciselor şi axa 


ordonatelor au respectiv coordonatele 
e i c 
(2-0) şi (0.7) pentru a, b #0. 
a b 


Dacă a = 0 atunci dreapta A de ecuație (1) este paralelă cu Os, iar 
dacă b = 0 atunci A este paralelă cu Oy. 


Dreapta A împarte planul în două semiplane 5, şi 82, situate 
de o parte şi de alta a dreptei A. Dacă pentru un punct M e8, 
coordonatele sale verifică inegalitatea 


ax + by <e (2) 
atunci coordonatele tuturor punctelor diù S, verifică aceeași inega- 
litate (2), iar coordonatele tuturor punctelor M e 8, verifică inega- 
litatea contrară : 

aw + by >c. (3) 


T4 


Deci soluţiile unei inecuaţii de forma 
az + by se (4) 


vor reprezenta coordonatele unor puncte din planul Oy situate 

într-unul din semiplanele determinate de dreapta A de ecuație 

(1) sau pe dreapta A. i 
Semiplanul pentru care există inegalitatea, 


az + by <e 


se determină verificînd sensul inegalităţi pentru un. punct M(z, y) 
aparținînd unuia din semiplane (de exemplu originea, dacă e + 0, 
adică dreapta A nu trece prin origine). În celălalt; semiplan vom avea 
verificată inegalitatea de sens contrar. 

Pentru un sistem de inegalităţi : 


aw + by sa . 
at + bY S Ca (5) 
Aa -+ bY < Cs 


reprezentarea grafică a mulțimii punctelor ale căror coordonate 
satisfac toate aceste inegalităţi se obține prin intersecția semipla- 
nelor (împreună cu frontierele respective) deduse din rezolvarea 
fiecărei inegalităţi în parte. 

De exemplu, să considerăm sistemul de inegalităţi : 


22 — 4y 


(6) 
4a + 


< 8 

—3% + 2y < 6 
8 

3 


6a -+ 5y 0 

Acest sistem poate fi adus la forma (5) cu toate inegalitățile de 
același sens prin schimbarea semnului în ambii membri ai inega- 
lității a treia, care devine 


— 4s —2y s—38 


Dacă unele din inegalităţi ar fi stricte, am considera semiplanul 
respectiv fără dreapta frontieră, ale cărei puncte verifică egalitatea 
dintre cei doi membri ai inegalităţii. 

Reprezentând grafic mulţimea punctelor din planul sOy ale 
căror coordonate verifică (6), se obţine un patrulater convex împreună 
cu interiorul acestuia (fig. 112.1), ale cărui laturi au fost îngroșate. 


y} 
ë Qi 
Eya 
$ 
ir 
1 (| 
7 NE p= x 
E = TeX : 
m Sa 
EA O 
E 
S = i 
Fig. II.2.1 Fig, II.2.2 


Dacă suprimăm însă din (4) inegalitatea 6% +- 5y < 30, se va 
obține un domeniu nemărginit* desenat în figura 1[.2.2. Vom arăta 
mai departe că mulțimea punctelor din plan ale căror coordonate 
verifică, sistemul de inegalităţi (5) are o proprietate caracteristică, 
şi anume aceea de a fi mulţime conveză,. - 

Prin definiţie, o mulțime S din plan sau din spațiu se numește 
multime convexă, dacă odată cu oricare două puncte ue și wes 
ea contine toate punctele segmentului care uneşte pe u cu u. 

Mulţimile de puncte desenate în figura II.2.1 şi 11.2.2 sînt 
mulțimi convexe. Mulţimile de puncte din figura TI.2.8 æ sînt convexe, 
iar mulțimile de puncte din figura IL.2.3b nu sînt convewe (mulţimi 
reprezentate prin porțiuni nehaçurate ale planului). 

- A XII-a figură convexă desenată reprezintă mulțimea punctelor 
(x, y) pentru care æ? < 2y, adică mulţimea punctelor interioare 


parabolei de ecuație y = Zar Mulțimea formată dintr-un singur 


* Prin definiție o mulțime din plan este mărginită dacă ea se găseşte în interiorul 
unui cerc cu centrul în origine. 
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punct este convexă, la fel mulțimea punctelor unui segment de 
dreaptă este convexă. În spațiu exemple de mulţimi convexe sînt 
mulţimea formată dintr-un singur punct, mulțimea punctelor unui 
segment de dreaptă, mulţimea punctelor unei drepte, mulțimea 
punctelor unui plan, unei sfere sau tetraedru (incluzind si mulțime: 
punctelor interioare suprafeţei respective). s 


y; 


Fig. IT, 2.3a Fig. I1.2.3b 


Pentru a da o exprimare analitică noțiunii de convexitate, să 
vedem cum se exprimă coordonatele punctelor unui segment în 
funcție de coordonatele extremităților sale. Pentru aceasta, să consi- 
derăm în planul sOy un segment AB astfel încit A să aibă coordo- 
natele (7, Yı); iar B să aibă coordonatele (w, Y1). Fie M un punct 


pe segmentul AB şi să notăm A -i (fig. II.2.4). Este clar că 


0 < à < 1, pentru à = 0 punctul M coincide cu extremitatea A a 
segmentuluilA B, iar pentru A=} punctul M coincide cu extremitatea 


a B(X, Yd 


ARJ) S 


a| A M B X 


Fig. 11.2.4 


Ba segmentului AB. Notind cu A’, M’, B', proiecţiile punctelor A, M, 
B pe axa Oz şi cu 0, D intersecțiile paralelei din A la Oa cu perpendi- 


cularele din M şi B pe Oz, obţinem : AU, = ^, deci AC = 
AD AB 


a AD = M'B' = Ma, — F). Abscisa punctului M va fi deci egală 


cu OM = 04' + A'M' = g t an — Aa = An H (L — Aa. Ana- 


11 


sertare taierea te A 


log obținem că ordonata punctului M este egală cu MM'=MO0+ 
+ CUM = ABD + Ya = My — VF Ya = Mk (LL — Aa: 

Deci dacă notăm vectorii de poziţie corespunzători extremităților 
A şi B ale segmentului prin u = (3, ja) şi U = (2, Yı), atunci vectorul 
corespunzător punctului M de pe segmentul AB cu proprietatea că 


—— =) va fi M+(1—a. Se poate arăta un rezultat 


analog şi pentru punctele situate pe segmente în spaţiu, în care 
caz vectorii au trei componente. i 

Acest rezultat ne permite să transcriem prima definiție a unei 
mulțimi convexe în forma următoare : E 

O mulţime S (din R? sau R?) se numeşte conveză dacă peniru 
oricare două puncte ue, weg şi pentru orice 0 < A < 1, există 
relaţia du -+ (L — Mie S. Am văzut că mulțimea soluţiilor posibile 
pentru o problemă de programare liniară este formată din mulțimea 
vectorilor + ale căror coordonate verifică inegalităţi de forma (5) 
pentru u = (w, y) e R?, sau de forma următoare pentru u = 
= (7,9, 2jeBs: 


da (1) 


<s 

< 
az + bay + C < d 

<x 


Teorema JI.2.1. 


Mulţimile de puncte ale căror coordonate verifică sistemele de 
inegalităţi (5), respectiv (7), sînt mulţimi conveze de puncte din plan, 
respectiv din spatiu. 

Vom face demonstraţia pentru mulțimile S din plan caracte- 
rizate de inegalităţile (5), demonstraţia pentru mulțimile din spaţiu 
caracterizate de (7) fiind analoagă. 

Să presupunem că două puncte ue 9, îi e 8, deci coordonatele 
lor verifică (5) şi obținem deci pentru u = (2, Yı) & = (y J1): 


aia + bY < 6y añ, + bY S cu 


(8) 


aat + bDi S Ca QT, + bah S los 


* se identifică ùn punct din plan (spaţiu) prin vectorul de poziție corespunzător, 
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Trebuie să arătăm că şi vectorul Au + (1 — A)u verifică (5) 
pentru orice 0 < à < 1. Într-adevăr, componentele acestui vector 
fiind Aa + (1 — A}Ū, respectiv Ayı + (1 — Aa, prima inegalitate 
din (5) pentru acest vector se serie : 


alag + (1 — AE] H ba H A L A sa (9) 


Această inegalitate se demonstrează prin regruparea termenilor, 
ţinînd seama de (8): 


ala + (1 — NE] + Bala + (1 — YI = 
= Maat, + Bays) + (1 — Da, + bar) < ae + (Î — de = 


deoarece 0 < à < 1,0 < 1 — à < 1. La fel se demonstrează şi cele- 
lalte inegalităţi din (5) pentru vectorul du + (1 — Aju, deci Au + 
+ (1 — Y e S pentru orice 0 < à < 1 c.et.d. 

Mulțimea convexă S formată din soluţiile inegalitäților (5) 
sau (7) se numeşte tronson. 

Dacă unele din inegalităţile nestricte (5) se transformă în inega- 
lităţi stricte sau chiar egalități, se verifică prin acelaşi calcul că 
mulţimea, S caracterizată de un sistem de egalităţi și inegalităţi 
liniare (stricte sau nu), este convexă, 

În spaţiu, fiecare din inegalităţile (5) determină un semispaţiu 
împreună cu planul de separație, iar intersecţia acestor semispajii 
este o mulțime convexă, mărginită sau nu, limitată de plane. 

Fiind dată o mulțime oarecare S din plan sau din spaţiu, se 
numește acoperire convexă a mulţimii 3 și se notează prin <S> 
cea mai mică mulțime convexă care conține mulţimea 8. Dacă 
mulţimea S este chiar convexă, atunci este clar că <S) = 8. 

De exemplu în plan acoperirea convexă a mulţimii formate 
dintr-un segment AB şi un punct C este triunghiul ABC (împreună 
cu punctele interioare) sau segmentul AB, dacă O se găseşte pe 
segmentul AB sau segmentul AC dacă B se găsește pe segmentul 
AC sau segmentul CB, dacă A se găsește pe segmentul OB. Acope- 
rirea convexă a mulțimii formate din două segmente AB şi CD este 
desenată în figura II.2.5 pentru câteva poziţii relative ale punctelor 
A, B, 0, D. 

Pentru a construi acoperirea convexă a unei mulțimi § se adaugă 
la S toate punctele situate pe segmentele care unesc puncte din S, 
pentru mulţimea astfel obţinută se repetă operaţia etc. pînă cînd 
ajungem la o mulțime convexă, care odată cu oricare două puncte 
ale sale conţine mulțimea punctelor segmentului care le unește. 
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Să observăm că distanța OP de la originea O a axelor de 
coordonate la dreapta d este dată de: 


O mulţime convexă P se numeşte poliedru convex dacă există o 
mulţime finită M care generează pe P, adică, pentru care P = (M). 

Poliedrele regulate din spaţiu : tetraedrul, cubul, octaedrul, 
dodecaedrul, icosaedrul regulat, sînt poliedre convexe în sensul 


i 
aa d ip ; n | 
acestei definiţii, ele fiind generate de mulțimea, vîrfurilor lor. $ | AN 
h ; | op — _OM:0N OM-ON b a M 
| = A à 
| | MN VOM:+ONE a Paz F va K 
5 a 


unde K este constanta Vaz + d2. În cazul cînd dreapta đ se găseşte 
în semiplanul S, determinat de inegalitatea as + by < 0, atunci 


Fig. 11.2.5 i 


i O mulţime convexă care nu este poliedru convex este de exemplu 


mulţimea punctelor suprafeței şi interiorului unei sfere. Un tronson 
care nu este poliedru convex este desenat în figura, II.2.2. 

În continuare vom arăta că extremele unei funcţii liniare de 
două variabile se ating în vîrfurile unui tronson din plan și vom 
prezenta algoritmul simplex pentru rezolvarea problemelor de pro- 
gramare liniară. 


1.3. Interpretarea geometrică a problemelor 
de programare liniară : 


În problemele de programare liniară trebuie să maximizăm 
(minimizăm) o funcţie liniară în care apar coordonatele unui punct 
u care verifică restricţiile problemei, adică aparţine unei mulțimi 
convexe determinate de anumite inegalităţi liniare. 

Vom arăta pe un exemplu că extremele unei funcţii liniare se 
ating în virturile mulțimii convexe determinate de restricțiile pro- 
blemei. 

Să considerăm o mulţime convexă D mărginită de drepte din 
planul Oy (fig. II.3.1). ` 

Să căutăm punctele 4 = (&, y) aparținîind mulțimii D pentru 
care funcţia aw + by ia cea mai mare, respectiv cea mai mică 
valoare, unde a, b sint numere reale. 

Să desenăm dreapta d de ecuaţie ax + by = ^, unde à este 
un parametru real şi să presupunem că dreapta d se găseşte în 
semiplanul S$, determinat de inegalitatea az -}- by > 0, adică A > 0. 
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Fig. 1.3.1 


; A 
à < 0 ṣi 0P = — Fa Fäcîndu-l pe à să varieze, se obțin diferite 


drepte paralele cu d. Din calculul anterior rezultă că valoarea à a 
funcției obiectiv aw -+ by în punctul u de coordonate (æ, y) e8, este 
egală cu ìà = 0P Va 4-52 = K.OP, adică constanta K înmulțită 
cu distanţa de la origine la dreapta paralelă cu d care trece prin 
punctul respectiv. Dacă (w, y) e Sp atunci A = — E.0P. 

Deoarece u trebuie să aparțină mulțimii D a punctelor care 
verifică restricțiile problemei, rezultă că minimul funcției am -+ by 
se atinge în virful A al lui D, iar maximul funcției as + by se 
atinge în vîrful B al lui D. Într-adevăr, distanța de la origine la 
dreapta paralelă cu d, care trece prin A este minimă, iar distanța 
de la origine la dreapta paralelă cu d care trece prin B este maximă, 


„dintre toate dreptele paralele cu d car€ trec prin puncte ue D. 


Maximul, respectiv minimul funcției az + by se pot atinge 
pentru toate punctele aparținînd unei laturi a mulțimii convexe D, 
dacă acea latură este paralelă cu d, deci există situații cînd o pro- 
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plemă de programare liniară are o infinitate de soluţii optime 
(fig. IT.3.2). 

Astfel, deoarece laturile AB şi A'B’ ale tronsonului convex D 
sînt paralele cu dreapta d, rezultă, că funcţia az -+ by este minimă 
pentru orice punct 4 = (4, y) apartinind segmentului AB și este 
maximă pentru orice punct u aparținind segmentului A'B’. 


Fig. 11.3.2 Fig. 1.3.3 


O altă situaţie care poate apărea în problemele de programare 
liniară este aceea cînd funcţia obiectiv az + by nu este mărginită 
pentru u = (4, y) e D, mulțimea, soluţiilor posibile ale problemei de 
programare liniară. Acest lucru se întîmplă cînd tronsonul D nu 
este mărginit. 

Astfel pentru tronsonul D din figura I[.3.3 funcţia aw + by 
îşi atinge minimul în A, dar poate lua valori oricît de mari, deoarece 
pentru. puncte u e D suficient de depărtate de origine distanța de 
la, origine la dreapta paralelă cu d care trece prin u poate fi făcută 
oricât de mare. Pentru tronsonul D’ maximul funcţiei az + by se 
atinge în B. Funcția az + by nu are un minim pentru u = (v, y)eD'. 
Într-adevăr, deoarece D’ are puncte oricît de depärtate de origine, 
situate în semiplanul determinat de az + by < 0, rezultă că distanța 
OP de la origine la dreapta paralelă cu & care trece prin ue D 
poate fi tăcută oricât de mare, deci valoarea funcției obiectiv 
à = — K.OP poate fi făcută oricît de mică. 

Rezultă că putem determina mulțimea soluţiilor unei probleme 
de programare liniară în două variabile cercetînd în planul Oy 
distanţele de la origine la dreptele paralele cu dreapta d de ecuaţie 
ax -+ by = 0 care trec prin viriurile tronsonului convex al soluțiilor 
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posibile ale problemei de programare liniară, cît şi i 
> tg 
pote de dreapta d în care se găsesc aceste virtari i caise 
| stiel pentru tronsonul convex D al soluţiilor posi ile din fi 
el ; ons t posibile din figura 

TE minia] funcției aw + by este atins în A, jar maximul zi B 

e se găseşte în semiplanul caracterizat de inegalitatea 
as + by < 0, iar B se găseşte în semiplanul caracterizat de inega- 


Fig. IL3.4 


litatea aw -+ by > 0, iar aceste virturi reali ă î i 
s > a ealizează în semiplaneli 
i ri maximele distanțelor de la origine la dreptele paralele 
Su, reap a e ecuaţie az + by = 0, care trec prin vîrturile tron- 
Am mai văzut că funcţia obiectiv îşi i 
Ă n ma şi poate atinge extremele 
n a infinitate de puncte aparţinînd mulţimii soluțiilor posibile 
ae problemei de programare liniară sau poate să nu fie mărginită 
erior (sau superior) pe mulțimea soluțiilor posibile. 


Egvemplul 1. 


Să rezolvăm prin metoda geometrică probl i i 
POE gala unda g problema de meniu optim 
Am văzut că formularea matematică a pr i ă 
ă 1 problemei este urmă- 
toarea : Să se determine numerele reale z, y astfel încât : 4 


2>0,y>0 
æ + 2y >3 
s+ 4y >44 
min (22 + Ty). 
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Reprezentînd grafic în planul sOy mulţimea convexă D caracte- 
rizată de inegalităţile e > 0, y > 0, @ + 2 >3 şi z + 4y > i se 
obține tronsonul din figura II.3.5 cu vîrfurile A o$), B(?, =) 
C (0,4). 


Dind dreptei d de ecuaţie 2% + Ty = 0 o mişcare de translație 
(deplasare paralelă cu ea însăşi), primul virf atins dintre A, B, € 


Fig. 11.3.5 


este vîrful B, deci ducînd prin A, B, C drepte paralele cu d, dreapta 
care trece prin B are cea mai mică distanță la origine. În concluzie 


L ` o kg 
coordonatele lui B şi anume g = 2, y S realizează minimul 


- funcției liniare 22 + Ty pentru (æ, y)e D, mulțimea convexă a 


soluțiilor posibile ale problemei. 3 i 
Faptul că dintre cele trei puncte aog) B24) C(0,4), 


punctul B realizează minimul funcţiei 2% +- Ty se poate verifica și 
prin calcul direct, deoarece notind f(s, y) = 2% + Ty obținem : 


1[o3)=3 (2 z) F, (0,4) = 28. 


2 
Deci soluția problemei dietei este & = 2, y = 3 sau, ținînd seama 
de semnificația acestor valori, meniul optim va conţine 2 unități 
de piine și 2 unităţi de brînză, adică 1000 g piine şi 250 g brînză 
(L unitate de măsură pentru piine = 1 unitate de măsură pentru 
brînză = 500 g). 
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Egemplul 2. 


Să rezolvăm prin aceeași metodă problema de transporturi 


enunțată la pagina 67 
“Să se determine « > 0, y > 0,2 > 0, t> 0 astfel încit 


æ +y = 1000 
z+t= 650 
æ+z< 700 
y +t< 950 
min (40 æ + 60y + 507 + 70t). 

Pentru a rezolva această problemă să reducem mai întîi numărul 
de variabile, folosind primele două, egalități din care deducem : 
y = 1000 — a 

(1) 
t = 650 — z. 
Înlocuind aceste valori în ultima inegalitate óbținem : œ + z> 
> 700, care împreună cu v -+z < 700 ne dă: 
æ + 2 = 700 


Pentru a rezolva geometric problema în planul gOz să vedem între 
ce limite pot varia œ şi z. Din datele problemei obținem æ > 0, 
2 > 0. Din ultima inegalitate obţinem y < 950, deci e =1000 —y > 50. 
Din ultima inegalitate obţinem de asemenea t < 950, deci 
z = 650 — t > — 300, care nu ne aduce nimic nou față de z > 0. 
Deci în planul xOz vom considera mulțimea convexă determinată de : 


s> 50, z> 0 
æ + z = 700 
care este desenată în figura II.3.6 şi care se reduce la mulţimea 


punctelor segmentului cu capetele în A (50, 650) şi B (700, 0). 
Funcția obiectiv devine, ţinînd seama de substituţiile făcute : 


40 + 60 (1000 — æ) + 502 + 70 (650 — z) = 105 500 — 20g — 
— 202 = 105 500 — 20(2 + 2) = 91 500, 
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adică este constantă pentru orice punct u = (xv, 2) e AB, mulţimea 
soluţiilor posibile ale problemei de programare liniară. Deci orice 
punct de pe segmentul A.B reprezintă o soluție optimă a problemei 
de transport. Se observă de altfel că funcţia obiectiv este minimă 


24 44 (50.650) 


Fig. 11.3.6 


atunci cînd suma œ + 2 este maximă, iar segmentul AB este paralel 
cu dreapta d de ecuaţie z+2=0. Considerînd Mm = (2 2) cu 
(2, 2&4) € AB, rezultă că æ, + 2 = 100 şi din relaţiile (1) deducem 
valorile celorlalte variabile : y, = 1000 — a, i = 650 — z astfel 
încât costul transportului să fie minim, egal cu 91 500 lei. De exemplu, 
luînd æ, = 50 t obținem z, = 650 î, y, = 950 t, t, = 0 t sau luind 
æ, = 300 ż obținem z, = 400 £, yı = 700 t, î, = 2507. 

În cazurile cind obținem mai multe soluții ale -unei probleme 
de programare liniară, putem alege una dintre ele după alte criterii. 

De exemplu dacă în acest caz întreprinderea I; poate expedia 
mai uşor cărămizile în C, decit în C, (există mijloace de transport) 
în timp ce I, nu exprimă nici o preferință, vom alege pentru s, 
valoarea maximă şi anume &, = 700 t, rezultînd 2, = 0, yı = 800 t, 
i = 650 t. 


I.4. Transformarea unei matrice prin metoda substituţiei 


Regulile de calcul prezentate în acest paragraf sint necesare în 
fundamentarea algoritmului simplex, însă servesc şi la rezolvarea 
altor probleme de algebră liniară, aşa cum vom vedea mai departe. 

Să presupunem că vectorii a; a2 aa din R? sînt liniar indepen- 
denți, iar vectorii arbitrari bı, ba, bg din R3 se exprimă în funcție 
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de a, az; a, în modul următor : 


Di = ant + oala + auda 


| ba = rgy + aaa F Osea | G) 


bs = casa F da3da + agsia 


Folosind coeficienții « € R (i, j = 1, 2, 3) formăm tabelul : 


| b ` ba bz 

i R 
ty | an Xa žig 
ng t oan za oa 
ag | aa aa oa 


Din modul de aranjare a elementelor, în tabel se observă că 
coeficienții exprimării vectorului b; în funcţie de vectorii lr, Ca ag 


` sînt tocmai elementele care se găsesc pe coloana, i, situată sub litera 


b: din tabel pentru i = 1, 2, 3. Am încadrat în tabel coloana a 
Il-a şi în egalităţile (1) egalitatea a II-a pentru a pune mai bine 
în evidenţă această corespondență. De asemenea fiecare coeficient 
4; în (1) se înmulțește cu elementul a, situat în tabel în aceeaşi 
linie cu el. 

Să presupunem că vrem să introducem în tabel vectorul b, în 
locul vectorului a. Pentru a fixa, ideile, să considerăm de exemplu 
k = 3 şil = 2, adică să obținem un tabel care are pe prima coloană 
literele a, ba, a. iar prima linie să rămînă neschimbată, adică să 
conţină literele b, by, ba. 

Condiţia pe care vrem s-o impunem unei astfel de substituții 
este ca şi vectorii a,, ba, & să fie liniar independenţi, la fel cum erau 
i, day as. Problemele pe care vrem să le rezolvăm sînt deci urmă- 
toarele : 

-~ P1) Să găsim condiția pentru ca o astfel de substituție să fie 
posibilă. 

P2) Să găsim noii coeficienți a; ai tabelului obţinut; prin 
substituție în funcţie de vechii coeficienți «y care apăreau în tabelul 
iniţial. 

Coeficienții &,, fiind numere reale, vom arăta, că această substi- 
tuţie este posibilă dacă şi numai dacă coeficientul «n, situat în 
matricea (&y)ij=z,2,8 pe linia lui a, şi coloana lui b, este diferit 
de zero. 
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Teorema IF.4.1. 


Putem substitui vectorul b, în locul vectorului a, în sistemul de 
vectori (a, ax, aa}, oblinând un sistem de vectori liniar independenţi. 
dacă și numai dacă coeficientul an, + 0. Dacă această condiție este 
depitităi formulele de transformare a coeficienţilor ay sînt urmă- 
toarele : 


& ` : 
ay = «y —— ay pentru î Æl şi 
Cry 


(2) 


Înainte de a face demonstraţia, să observăm că formulele (2) 
au sens, deoarece ay, 0. Blementele de pe linia lui a, (vectorul 
care este substituit de b,) se împart cu coeficientul a, care se 
numeşte şi pivot al transformării, iar cele de pe celelalte linii se 
transformă cu regula dreptunghiului : 


j k 


Fig. 1.4.1 


Se scade din vechea valoare e, produsul elementelor om şi a, 
„proiecțiile“ lui «y pe coloana şi respectiv linia pivotului, împărțit 
la pivot. 

Să observăm că nu mai este nevoie să aplicăm această regulă 
pentru elementele situate pe coloana pivotului, deoarece din (2) 
obţinem ap = 0 pentru i #1 şi ay = 1. Se obţine astfel un nou 
tabel de forma. : 


| LA ba bs 
pe Xa a'ia o 
ba xan O'22 i 
as % 33 0'a2 0 
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în cazul în care 1=2, k = 3, coeficienții æy fiind daţi de (2). 

Vom demonstra teorema în cazul | = 2 și k = 3, demonstraţia 
făcîndu-se analog pentru orice indici I şi k. 

Să arătăm mai întîi că dacă «s # 0, atunci vectorii aa, ba, as sînt 
liniar independenţi. 

Conform definiţiei vectorilor liniar independenţi, trebuie să 
arătăm că dacă 


Aa + Aoba + dada = 0 (3) 


atunci coeficienții M = às = ħa = 0. 
Din egalitățile (1) obținem că b, = ajg + asta H Casa, de 
unde substituind în (3) pe b, cu această expresie se obține : 


(Ar F Aa das) Ga + Aztag Ga F (àg ocaa)aa = 0. (4) 


Deoarece prin ipoteză vectorii a, da; @ sînt liniar independenți, 


rezultă că: 


M T àz% = 0, 
nas = 0, (+) 
àg F raag = 0. 


Din a Il-a ecuație, deoarece «as z 0, se deduce ^, = 0, valoare 
care înlocuită în celelalte două ecuații ne conduce la M = à, = 0, 
deci vectorii a, ba, ag sînt liniar independenţi. 

Pentru a arăta suficiența condiţiei, adică faptul că vectorii 
Gu bs, âs fiind liniar independenţi implică «s, + 0, să demonstrăm, 
prin reducere la absurd, că dacă ass = 0, atunci vectorii au, bs; ds 
sînt liniar dependenţi. 


Într-adevăr, din ecuaţiile (1) obţinem b = cs di +&s3 Gay adică. 


o ecuaţie de forma (3) : 43 a — ba + das da = 0, în care coeficientul 
lui b, este diferit; de zero, deci a,, da, aş sînt liniar dependenţi. 

Mai rămîne să arătăm că noii coeficienți a; daţi de formulele 
(2) verifică relaţiile de forma (1) asociate însă tabelului transformat, 
adică : 


, , , 
b, = ai dak an bs + sa a 


1 r 
ba = dis Qi + daa Da F ase a (1) 


Ultima, egalitate fiind verificată, să arătăm de exemplu că are 
loe prima egalitate. 
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Pentru aceasta, să luăm ultima egalitate din (1): 
ba = dis a + das Ga + aa da. 


Din această egalitate, deducem exprimarea vectorului a, în funcţie 
de vectorii a, ba şi a: 


Substituind această valoare în prima egalitate din (1) deducem 
exprimarea vectorului b, numai în funcție de a, ba Şi aa: 


Da = an af ai Ga + aa ds = aa da F ta 


a IL « 

13 33 

| a + bs —- je F&a ds. 
Cos X23 das 


Regrupind termenii obţinem : 


1302 ` A Casa 
b, du a“ + ba + | aaa A 


azs a i Cos 


Prin identificarea, coeficienţilor găsim noile valori : 


Ajg 
Pan 13 %21 
Luci e 1 E 
aa 
a 
, 2: 
an =; 
X23 
&33 01 
Ası = da 
zh 


care coincid cu valorile date đe (2) pentru 


i=1,3; j=1; k=3; l=2. 


În mod asemănător se deduce și a doua egalitate din (1). 
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Aplicaţii la rezolvarea sistemelor de ecuaţii liniare şi la inver- 
sarea matricelor 


Vom arăta în continuare cum poate fi aplicată metoda substi- 
tuției la rezolvarea sistemelor de ecuații liniare şi la inversarea 
matricelor. 


Să considerăm mai întîi un sistem de trei ecuații liniare cu trei 
necunoscute de forma : . 


Aaa + hY + das? = bu 


aat -F AY F ao = ba (5) 


at + aay F aa = ba. ` 


În cazul coeficienţilor a; î, j = 1, 2, 3, se observă că primul indice 
indică linia, iar al II-lea, indică coloana în care se găsește coeficientul 
respectiv în sistem. Deci coeficienții necunoscutelor v, y, z, formează 
o matrice : 


Gia Aiz Aig 


4 = Aaa 


da 023 


Aaa @za (a 


Elementele situate într-o aceeaşi coloană a matricei A formează 
fiecare cite un vector cu trei componente. Aceşti vectori îi notăm : 


da = (tu Ga da) 
da = (n, Qaz az) 
ag = (Qis, Qes; Asa), ă 
iar vectorul termenilor liberi îl notăm : 
b = (bus ba, ba) 
Cu aceste notații, sistemul (5) se serie : 
oa + Ya, + Zaz = b. (6) 
Într-adevăr, aceasta este o egalitate între doi vectori cu cîte trei 


componente, iar doi vectori sînt egali prin definiție cînd compo- 
nentele lor de acelaşi ordin sînt respectiv egale. Egalind respectiv 
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primele, cele de a doua şi cele de a treia componente ale vectorilor 
din cei doi membri ai egalităţii (6), se obţin tocmai cele trei ecuaţii 
ale sistemului (5). Rezolvarea sistemului (5) revine deci la a exprima 
vectorul termenilor liberi b, ca o combinaţie liniară de vectorii 
formați cu elementele situate pe cîte o coloană a matricei A a coefi- 
cienţilor sistemului. 

Coeficienții acestei combinaţii liniare vor fi tocmai soluţia siste- 
mului (5), 

Pentru a rezolva sistemul (5) vom pleca de la tabelul : 


| a ag a |b 
4 Cu fig Tis b 
H Ca Aag aa ba 
U3 iy ij 033 ba 


unde, Ug, Ug sînt; vectorii unitari din R? : 
u = (1, 0, 0), 
u = (0, 1, 0), 
us = (0, 0, 1), 


care sînt liniar independenţi, aşa cum cere metoda substituției, iar 
vectorii d, aa, ag și b se exprimă în funcție de acești vectori unitari 
cu coeficienții care apar în coloana respectivă din tabel, de 
exemplu : 

b = bity -H bata + bata. 


Să observăm că matricea care apare în acest tabel este tocmai ma- 
tricea A a sistemului. 

Încercînd să substituim pe rînd vectorii unitari cu vectorii 
Qi, do Qş (nu contează în ce ordine), folosind de fiecare dată un pivot 
nenul şi aplicînd formulele (2) de transformare a tabelului, putem 
ajunge la situaţia cînd vectorii Ur, 22; Ug au fost substituiţi toţi cu 
vectorii a, da, ds, obţinînd tabelul: 


| ti a3 a | b 
a 1 o o | & 
ag 0 1 0 | 8 
az 0 0 1 | B 


În acest caz b= Ba, + Pata + fsas, conform definiției unui 
astfel de tabel. 

Deci coeficienții unei astfel de exprimări sînt tocmai soluţia 
sistemului (5), adică 


& = Pu Y = Ba Z = Ps- 
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Dacă făceam substituția într-o altă ordine în tabel, prima coloană 
conţinea vectorii a, 2, & într-o altă ordine, însă întotdeauna, valoarea, 
numerică a soluţiei æ se găseşte în linia Iui a, valoarea lui y se găsește 
în linia lui a, iar valoarea lui z în linia lui as. 

Se poate să ajungem însă la un tabel în care substituţia nu poate 
continua pînă la capăt, deci pentru care nu mai găsim pivoţi nenuli. 
Să presupunem că am ajuns la situaţia : 


| a d a |b 
& 1 0 cs Bi 
tı 0 1 Ces Pa 
GA 0 0 9 a 


Nu-l mai putem substitui pe 43 CU ag, deoarece în linia lui u, toate 
elementele tabelului sînt nule. 

În această situaţie, sînt posibile două cazuri : 

1) «æ = 0 şi în acest caz o soluţie a sistemului este : 


2 = Bu Y = Bp2=0. 


Se poate arăta* că în acest caz sistemul este compatibil şi 
nedeterminat, adică admite o infinitate de soluţii. 

2) ax 0 şi în acest caz sistemul este incompatibil, deoarece 
vectorul b nu se poate exprima ca o combinaţie liniară de vectorii 
coloană a; Ga, da, adică sistemul (5) nu are soluţie. y 

În mod analog se tratează şi cazul cînd substituţia nu mai 
poate continua după ce am aplicat o singură iterație, adică după 
substituirea lui 4, cu a, am obținut un tabel de forma : 


l m al a; | b 
Si 1 kia ĉia Pı 
GA 0 0 
uz 0 0 0 Y 


În acest caz, dacă « = y = 0, sistemul este compatibil şi nedeter- 
minat, o soluție "particulară fiind e = Bi y = 0, z =0. 

În caz contrar, adică dacă măcar unul dintre numerele æ şi 'Y 
este nenul, atunci sistemul este incompatibil. 

Folosind. metoda substituţiei sîntem deci conduși la una din 
cele trei situații posibile : sistem compatibil şi determinat, cînd obținem 
soluția unică a sistemului, sistem compatibil și nedeterminat, cînd 
obținem o soluţie particulară a sistemului și sistem incompatibil. 


* Demonstrația acestei proprietăţi se face ținind seama că în acest caz rangul 
matricei coeficienţilor este doi. 
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Vom vedea pe un exemplu că această metodă este chiar metoda 
substituţiei unei necunoscute dintr-o ecuaţie în celelalte ecuaţii 
ale sistemului. i 


Avantajul metodei constă într-o scriere mai compactă a coefi- 


cienţilor ecuaţiilor obținute prin substituţii, fără a mai serie şi necu- 
noscutele respective. 

Ordinea, de efectuare a operaţiilor este următoarea, : 

1) Se alege un pivot nenul. 

2) Se transformă toate elementele tabelului, inclusiv coloana 
corespunzătoare lui b, elemente care sînt situate pe linii şi coloane 
diferite de linia şi coloana pivotului, după regula dreptunghiului. 

3) Se împart elementele liniei pivotului la pivot. 

4) Se înlocuiese elementele situate pe coloana pivotului prin 


zero, mai puţin elementul situat pe linia pivotului, care a fost; înlocuit; . 


la 3) prin 1. 

Aceste transformări se pot face pe un același tabel în această 
ordine, prin ştergerea vechilor elemente şi înlocuirea lor cu noile 
elemente calculate la 2), 3) şi 4). Acesta estemodulîn care operează 
an calculator electronic. 


Exemplul 1. 
Să considerăm sistemul: 
3p+y—38a=—2 
— æ+ y+ 5z=4 (7) 
2æ — Ty — z= — 1 
Vectorii formați cu elementele situate pe coloanele matricei A a 
coeficienţilor sînt în acest caz a = (3, —1, 2), a = (1, 1,—7) şi 


as = (—3, 5, —1). 
Se formează tabelul : 


| a az a | b 
a E 1 o 
u, | ~ 1 5 4 
LA 2 —7 —1 —1 
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Alegînd ca pivot numărul 3 încadrat în tabel cu regulile 2) — 4) 


se obţine noul tabel: 


a, CA CA b 
TEE 
3 

|4 14 
a E EN 
23 7 
op sa ES 


Alegînd ca nou pivot numărul 4/3 şi făcînd transformările 2) — 4} 
se obţine tabelul următor : 


a, CA LA b 
a 1 0 —2 |—1/2 
aa o 1 3 |72 
DR 0 0 |24 | 49/2 


Alegînd ca pivot numărul 24 obţinem tabelul final : 


bei 


da 


a | 


b 


& 
a 


0 


1 
0 


0 
0 
1 


37/24 
21148 
49/48 


aa 
deunde se citește soluția sistemului: œ = 37/24, y = 21/48, g = 49/48. 
“Dacă împärțim prima ecuaţie a sistemului (7) la 3, se obține 


Ş 1 A 
a+ Š y— z=, jar dacă il inlocuim pe = — >y +2 + 2 în 
3 


celelalte două ecuații ale sistemului, obținem sistemul echivalent : 


EREE E 
cir zi 3 
14 

—y + 4z= — (72) 
3 
— Bz 
3 


ai cărui coeficienţi apar în cel de al doilea, tabel. Dacă se împarte 


acum a II-a ecuaţie prin z şi se substituie y din a II-a ecuație 


a sistemului (7”) în prima şi a III-a ecuație a lui (7) se obţine sistemul 
echivalent (7%): 


a — 22 = — 1/2 E 
y + 3a = 1/2 r) 
242 = 49/2 


95 


ai cărui coeficienţi şi termeni liberi apar în cel de al III-lea tabel. 

Împărţind în sfirșit a III-a, ecuaţie din (7) cu 24 şi substituind 
valoarea lui z astfel găsită în primele două ecuații ale lui (7”) obţinem 
sistemul echivalent 


w = 37/24 
y = 21/48 (°) 
2 = 49/48 


de unde se citesc soluțiile şi care are drept coeficienți şi termeni liberi 
elementele ultimului tabel obținut prin metoda substituției. 


Exemplul 2. 


Fie sistemul: 


æ +y tz =2 
— &+yt3z=4 (8) 

2% — 2z = 0 

Formăm ca mai înainte tabelul: 
a LA a; b 
m a f2 
DA -1 1 3 4 
DA 2 o —2 0 


Alegind ca pivot numărul 1 încadrat în tabel se obține prin” 


substituție noul tabel : 


CA üg LA b 
a 1 ES 1 
A 0 |2] 4 6 
DA 0 --2 —4 |—4 


Substituind pe æ, în locul lui t, se obţine noul tabel: 


LA LA da b 
a 1 o —i |— 
az 0 1 2 3 
uz 0 0 0 2 


Procedeul substituţiei nu mai poate continua deoarece pe linia 
Tui u şi coloanele Mi a, 2, as se găsesc zerouri. Sistemul este incom- 
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palibil, deoarece în linia Iui u, și coloana lui b se găseşte numărul 2, 
care este diferit: de zero. 

În cazul inversării unei matrice A = (@y)ij.ig SC cere să se 
determine o matrice B = (by sosea astfel încit AB = E, unde B 
este matricea unitate cu trei linii şi trei coloane : 


1 9 ð 
E=[0 i 0 
0 0 1 


Notînd vectori linie ai matricei E cu w =, 0, 0), ua = (0, 1, 0), 
ua = (0, 0, 1) şi ținiud seama de definiţia. înmulţirii matriceale, 
se obține : 


ba | ba by ije o 


ba | ba Pa = loli o 


ba boa ba 0|o 1 


tbir + haba F Gaba = 1, 


anbi F toba F Qaaba = 0, (9) 
ab F Gabr + agba = 0; 
prin egalarea elementelor primei coloane din produsul AB gi din 
matricea P. 
'Pinînd seama că elementele primei coloane a matricei A 
formează vectorul : 


a = (Gus; Ca; Gah 


-gi usemănător 


i 


Qa == (Min, Aro, aa), 
ag = (0133 aa Gas), 
egalităţile (9) reprezintă tocmai egalitatea vectorială : 


bnt F- baz + bas = th 


T- e. 2101 97 


a a ei 


scrisă pe componente. Se obţine în mod analog, că 
bioty F bastz + Basta = Uz 
biat, -F baza + basla = tgs 


egalităţi care ne arată că elementele matricei inverse B = A”! sînt 
tocmai coeficienții exprimării vectorilor unitari t, s, ug, că nişte 
combinații liniare de vectorii formați cu elementele situate în 
coloanele matricei A. 

Deci pentru a calcula inversa matricei A (dacă aceasta există), 
procedăm astfel : 

Plecăm de la tabelul : 


| a a a | mw da Up 
(Ei au da aya 1 9 0 
LA au azy azg 9 1 D 
Us au da a33 o 0 i 


şi căutăm să substituim vectorii unitari 44, 12, %a Serigi în prima 
coloană cu vectorii a, Qa, 4. 

Dacă această substituire se tace pină la capăt, se va obţine 
tabelul : 


| a A ds | H w Us 
a 1 D o bu ba Dia 
CA 0 1 è bat baz bog 
& 0 0 1 da baa baz 


care conține în partea dreaptă tocmai matricea B = A71 Dacă 
această substituire nu se poate face pină la capăt (obținîind o linie 
numai cu zerouri în dreptul lui 4, >, ds), atunci matricea A este 
singulară (nu are inversă), deci det A = 0. 


Baemplul 3. 
Să inversăm matricea : 


1 2 —1` 
4A=|30 a 
1 1 4 
prin metoda substituţiei. 
Formăm tabelul : 
| a LA d | m LA LA 
t h] 2 
a 3 0 
Ea 1 1 
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Luind ca pivot numărul 1 încadrat în tabel obţinem : 


a LA [A în Ua uz 
a 1 2 i 1 o o 
t |—6] 5 | -3 1 0 
tlg 0 —1 5 ict 0 1 


a LA CAS LA Us i 
a 1 o 2/3 0 18 o 
aa o 1 —516 | 42 -16 o0 
A 0 0 |256| |-a2 = 1 


Alegînd ca pivot numărul 25/6 se obţine tabelul final : 


LA LA a | m LA LA 
a 1 o 0 [2125 9/25 —4/25 
a o 1 0- |25  —1/5 15 
az o 0 1 i—8/25 —1/25 6/25 
Din acest tabel deducem matricea inversă, : 
2/25 9/25 — 4/25 
A™ =| ap —1/5 1/5 
—3/25 —1/25 6/25 
Exemplul 4. 
Să se inverseze matricea : 
1 0 1 
A=|—1 2 4 
0 2 5 
Se formează tabelul : 
| a LA d | u LA ig 


Alegînd ca pivot numărul 1 încadrat se obține : 


a az Hy Hi Uz us 
a 1 o 1 1 [ o 
ua o 2] . 5 1 1 0 
u o 2 5 0 0 1 

' Substituind pe uz cu a, găsim tabelul: 

d az 1 u ls 
a 1 o 1 o 0 
as o 1 5/2 12 12 0 
tlg 0 0 0 —1 —1 1 


pentru care substituţia nu mai poate continua, deoarece elementele 
situate pe linia lui up şi în coloanele lui a, aa, aa Sint nule. 

Deci această matrice A nu este inversabilă. 

Într-adevăr, matricea A are Qeterminantul egal cu zero, deoarece 
ultima linie este egală cu suma primelor două linii. 


Notă: Metodele expuse se pot utiliza pentru cazul general al rezol- 
vării unui sistem de m ecuaţii liniare cu n necunoscute și la 
inversarea unei matrice cu n linii şi n coloane, cu elemente 


numere reale. 


Baeroiţii 


Folosind metoda substituţiei, să se rezolve următoarele sisteme 


de ecuaţii liniare : 


1. æ + W= 


—s -+ y=5 
R: («= —1, y= 4) 
3. æ -+ 3y = — 1 


—4g — öy = 0 
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2. 


4 


32 —y=t 

6a — 2y =8 

R: Sistem nedeterminat 
æ—y=3 


4g — 4y = 1 


R: Sistem incompatibil 


5. s+y—22=0 
—g — y+ Tz=7 


R: Sistem nedeterminat 


11. &— y— z=1 
3s +- Ty — z=4 
6s + 4y — 4 =" 

13. @+ y+ z+ t=4 
2% — y+ z— i=1 
—g + 2y + 3z — 3t =1 

p— y+Ta— t=6 


6. 3s —2dy=a4 
s+ 2=3 
y— 2=2 


8. 2% + 3y — 2=0 
—t+ y+ 22=1 
x+ 4y + 2=2 


R: Sistem incompatibil 


19. 2—y—2=3 
- 4 +y—z=5 
—3g —y+z=—_1 
R: @ = 4, y = — 5, z = 8) 
12. 2% + 3y + z=0 
— t y— z=3 


—2æ + Ty — 3z = 1 


R= ekr L tel 


2% -+ 2y + =4 
— æ — y+38i=3 


14, —3a y+ z—i=5 


15. a y+.z— t=i1i 16. Ta — y— 2421 =—10 


2% + y—382=3 —3g + 4y + õ5t=1 
g y + 22 — 2t = 4 2g — y— z+ t=0 
3e — W+ z—4t=9 8æ + y —3z+ 9% =4 


R: Sistem nedeterminat R: Sistem incompatibil 


Folosind metoda substituţiei, să se inverseze următoarele matrice : 


17, 1 =F 18. 1 4 19. ( =] 
(, :) 5 20 i 3 4 
R: 5/8 1/8 
E ia) 
20. 1 1 4 21. 0 1 —3 22. (1 2 3 
2 0 3 —1 2 4 4 5 $ 
—1 ł 2 0 1 0 5 T 9 
R: 342 —1i —3/2 R: Matrice singulară 
742 —3 —5j2 
—i 1 1 x 
23. (a 12 24. 1 1 1.1 
4 —4 3 —1 1 1 1 
2 —2 7 1 — 1 1 
1 1 —l 1 
R: Matrice singulară 1 EE INA 3 
4 4 44 
= N 41 
R: 4 4 4 4 
LaLa 
4 4 4 4 
-i111 
4 4 4 4 
25. 1 2 0 = 26. (1 4 — 
2 1 0 3 3 2 1 
—1 0 1 1 0 1 5 
1 3 2 4 


R: Matrite singulară 
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27. Să se evalueze numărul de operaţii la rezolvarea unui sistem de 
trej ecuaţii algebrice liniare cu trei necunoscute prin metoda substi- 
tuției. 


R: 12 adunări şi scăderi, 12 înmulţiri, 18 împărțiri, pentru sistem 
compatibil şi determinat. 


28. Să se evalueze numărul Ge operații la inversarea unei matrice 
cu trei linii şi trei coloane prin metoda substituţiei. 


R: 12 adunări şi scăderi, 18 înmulţiri, 27 împărțiri pentru o matrice 
nesingulară. ` 


H.5. Metoda simplex pentru rezolvarea problemelor 
de programare liniară 


Fentru calculul unei soluții optime a problemelor de programare 
liniară vom prezenta un algoritm sistematice care conduce la soluții 
posibile din ce în ce mai bune, pînă, cînd se ajunge la o soluție optimă, 
sau la concluzia, că problema nu admite un optim finit. 

Să considerăm următoarea problemă de programare liniară : 


Problema 1 (forma canonică) : 

Să se determine vectorul u = (2, Y, 2) ER? cu componente numere 
reale nenegative £, Y, z > 0 astfel încât să fie verificate restricțiile 
următoare : 

ant + aY + ns? S bys 
aat -F loy F Aoa? S bas (1) 
aat F Aay H lg? S bas l 


unde coeficienții Qu, i, J = 1, 2, 3 sînt, numere reale, iar bi, ba, ba 
sînt numere reale nenegative*, cu condiția ca funcţia obiectiv ce -+ 
+ Ext -+ ez să fie masimă. : 

- Ne putem rezuma la problema de maximizare, deoarece problema 
de minimizare se reduce la aceasta. Într-adevăr, 


min (ag +F Cy + 032) MAX (— E — CY — Cg). (2) 


* Cazul în care b, ba, b, sînt numere oarecare nu va fi tratat aici. 


Justificarea acestei egalităţi se face în felul următor : ` 

Presupunînd că expresia cw -+ CY + oa poate lua valori nega- 
tive, atunci minimul ei se obține pentru acele valori (x, y, 2) pentru 
care 0% + eY -+ eg < 0, cînd ags -+ ey + e este maximă în 
valoare absolută, adică : 


min (040 + ey + e) = — max |as + ey + e] = 
= — max (— 68 — Cs — 632). 
În caz contrar, expresia cw + Cy -+ eaz ia numai valori pozitive, 


deci — 04% — ey — 6? < O şi deci maximul acestei expresii — 4% — 
— cy — 6z se obține pentru acele valori (x, y, 2) pentru care expresia 


E 039 F ee = | — cat — ey — cae | este minimă. În acest caz, 
putem serie de asemenea : ; 
min (04 -+ CY H- 092) = — MAX (— 0U — êY — 032), 


deci egalitatea (2) este demonstrată. 

Pentru a rezolva această problemă cu algoritmul simplex, vom 
introduce mai întîi variabilele de compensare sau variabilele-ecart 
z, J, Z > 0, care transformă restricțiile (1) în egalitäți : 


ant + aY + dat + d = by „aa 
aat F Goa + das HY = ba (3) 
lat + agy F Agg +E = bs ; < 


Am transformat astfel problema inițială de programare liniară în 
problema 2 (forma standard): . X 

Să se determine numerele reale nenegative e, Y, Z, Ü, Yı Z >Ô, 
care verifică restricțiile (3) și pentru care functia obiectiv 0.20 + CY +6 
este maximă. 

Dacă (25 Yu 2» By Yy 2) este o soluție optimă pentru această 
problemă de programare liniară, atunci (® Yy 2) este o soluție 
optimă a problemei inițiale de programare liniară, 

Într-adevăr, deoarece (Vn Yı, 2u Bu Yy 2) verifică egâlităţile (3), 
rezultă că (d, Wu 2%) verifică inegalitățile (1), deoarece 4, > 0, 
Jı > 0, 2 > 0, şi în plus expresia Crt, + Cata F 03%, este maximă. 
Problema de programare liniară caracterizată de restricţiile (3) 
admite o soluţie nenegativă și anume : 


bu Ya cat bas 2 => ba. 


2 > 0, y =0, a =0, li 
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Algoritmul simplex porneşte cu această soluţie posibilă, numită soluție 
pg de bază, pe care o îmbunătăţeşte succesiv, făcînd calculele 
în tabelul : É 


t GA CA u, LA n, | b 
Hy an aig 013 1 0 o |. 
Hg Ca azg aag 0 1 9 bo 
LA az CAN a33 0 0 1 LA 
ei Ca — 6 0 9 9 0 
Tabelul 1 


În acest tabel am notat prin a, dz, aş vectorii formați cu elemen- 
tele situate pe coloanele matricei A a coeficienţilor æ care apar în 
restricţiile problemei inițiale : 


Cu Aa Wa 
A= |an ün a23 
Ca Tia Qas 


adică a, = (Gn, Ga Oz), de = (Qio, Qoz, Ag) ȘI do = (G13 ea tss). 
Vectorii 4, a, Ug sint vectorii unitari din R3, 


u, = (1, 0, 0), u= (0, 1, 0), «= (0, 0, 1 


care sînt liniar independenţi, iar vectorul b = (b,, ba, ba). Primele 
trei linii ale tabelului sînt alcătuite după regula introdusă la metoda 
substituţiei : vectorii a, Ca, @g, Un, Ua, Ug Se exprimă ca nişte combi- 
naţii liniare de vectorii care apar în prima coloană şi anume th, ta 
ua, cu coeficienţi care sînt; serişi în coloana vectorului corespunzător : 


a, = Os + Gata + Aaa 
a= Canta + Calle + Agta 


a3 = Casta -+ Caatla + Casta 


th = Lua + Otta + Ottis, 
ta = Ota F Lug + Otia, 


(4) 


Ug = Oty + 0ua + Itiz; 
b = byu, + Bata + bats. 


Aceste egalităţi se verifică folosind definițiile operaţiilor vectoriale. 
Astfel : 


a = an (1, 0, 0) -+ aa (0, 1, 0) + asa (0, 0, 1) = (a 0, 0) -+ 
+ (0, an, 0) + (0, 0, an) = (äu, ass €x) Celelalte egalități din (4) 
verificindu-se analog. Numerele din ultima linie a tabelului se cal- 
culează tinind seama de funcţia obiectiv a problemei 2 care este : 


OL + CY + C = Cu + Cy + Coe + 05 + 0y + 0z. (5) 


Numerele din ultima linie a tabelului se calculează după regula 
următoare : 

Pe coloana lui a, se înmulțește a cu coeficientul lui & din 
fúncția obiectiv (5) (care este zero), se adună cu a înmulţit cu 
coeficientul lui y din funcţia, obiectiv (care este zero), iar rezultatul 
se adună cu ay înmulţit cu coeficientul lui z din funcţia obiectiv 
(care este zero). Din rezultatul obținut se scade coeficientul lui z 
din funcția obiectiv pentru problema de programare liniară 2, care 
este c,. Calculele se fac analog pentru celelalte coloane ale tabelului. 

Deci putem serie : 


coeficientul lui 2 


| | RT lui 2 RE | 


coeficientul lui şi. 
coeficientul lui ¥_ 


— Cp = Ian + Oaza + Dag -a | 
_ coeficientul lui y 
— ia = Oye H Oaza H Oase ta 
__eoeficientul lui z 

eg |” + |” + Oaa —c3 

0 = 0-1 +-0:0+ 0-0 —0 

0 = 0:0 4-0:1+0-0 —0 

0 = 0:040:0+01 —0 


Pentru a înțelege mai bine modul de calcul al elementelor 
situate în ultima linie a tabelului 1 în dreptul coloanelor lui 4, az, 
lg, ti Ua Us, Să stabilim următoarea corespondenţă, : a, = V, day, 
ag Z, mă, tao, wg 2. Elementele care se scad vor fi mereu 
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aceleași și anume coeficienţii lui v, y, z, &, y, 2 din funcţia obiectiv, 
adică C, Cz; 03 0, 0, 0 

Elementele cu care se înmulțesc numerele 4, depind de vectorii 
care se găsesc în prima coloană. De exemplu, dacă în prima coloană 
a tabelului în loc de t, tp, tg S-ar fi găsit de exemplu d; Uz da 
prin corespondenţa stabilită mai sus, trebuia să înmulțim numerele 


a; respectiv cu coeficienţii lui æ, y, z din funcţia obiectiv, care sînt, 


respectiv c, 0, 63. În acest caz, în ultima linie a tabelului am îi avut 
armătoarele numere : £ 


da >| Cana F Otar + Cotar | — 6r 
Mot REA 
de Cilio F Oloo F Callaz Cz 
+ T $ 
da —> | Callas -+ Oaza +-.0aaa — 63 
oa t 
d, >| 6-1 + 0-0 + 63-0 —0 
ES PE e 
d; >| 00 + 0-1 + 3-0 —.0 
Î t t 
de >| 0 +00 + ezl . —0 
coeficientul lui z 


|___coeticientul lui y 
coeficientul lui œ 


| E 


Mai departe vom nota prin di, do, da, da, ds, de sumele încadrate 
în figura anterioară. 

Valoarea 0 din ultima linie și ultima coloană a tabelului 1 
reprezintă valoarea funcției obiectiv os + Cay + Ca? pentru soluţia 
posibilă iniţială a problemei de programare liniară, care este (0, 0, 
0, bi, ba, ba), adică; valoarea funcţiei liniare ce ~- Cay + caz + 02 + 
+ 0y + 02 pentru g = y = z = 0, adică zero. În continuare algo- 
ritmul urmăreşte să înlocuiască vectorii unitari 24, Us, Ug Care apar 
în prima coloană a tabelului prin vectorii a; da 23, după anumite 
reguli care vor fi expuse mai departe, pînă se va ajunge fie la o 
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soluţie optimă a problemei de programare liniară, fie la concluzia 
că problema nu are maxim finit, adică funcția liniară o -+ cay + Caz 
poate lua valori oricît de mari pentru soluţiile posibile ale problemei 
de programare liniară. Să presupunem că am înlocuit vectorii «4 şi 
üa din baza lui R? prin vectorii a, Şi aa, obţinînd după doi paşi de 
aplicare a algoritmului simplex tabelul 2 : . 


Li de 03 uy Ua EA b 
be 1 0 93 Zia is l o Bı 
az 0 1 B23 ča az 0 Ba 
us 0 0 | ag | Xa as 1 Bs 

0 0 da — cs de de 0 v 

Tabelul 2 


unde vectorii a, 3; ug din R? sînt; liniar independenți. 
Pentru acest tabel regula de formare este aceeași cu a tabelului 
iniţial, adică există, relaţiile : 
a = L-as + 0-a + Otia 
a = 0-a, + L-a + 0-u 
da = Argl, F %asda + astia 
Ur = Qut F oala F Asilis, : d (6) 
Ua = asi F asia F asa i 
Ug = 0-a, + 0 -aa + L tig; 
b= Bi + Pato + Pats. 


Elementele din ultima linie a tabelului se calculează astfel : 


0 La H 0e H00 e. 
0 = 0-0, + 1e + 0-0 — êy îi: 
dg — Ca = ra C1 F dosta E aas: 0, 


da = cela + aab + si 0 —0, 


dg = asa H aofa + ass:0 — 0, w A K 
0 = 0.0, +0 +100, : 
Vi = Bau + Pala + 0-03. 
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Valoarea +, a funcţiei obiectiv se calculează ţinind seama că 
soluţia nenegativă corespunzătoare noii baze (a, az, %3) este (Bi 
„Ba > 0; 0, 0, fa), deoarece din (6) obţinem: b = fa. + aaa -+ O -as + 
F O- F Ots + Bata şi ea se citeşte din tabel, pinînd seama de 
vectorii care intră în noua bază. 2, este valoarea funcției obiectiv 
pentru această nouă soluție nenegativă de bază, iar pe ultima linie 
“a-“tăbelului diferenţele d; — e, corespunzătoare vectorilor din bază 
sînt nule. i i 

Desfășurarea algoritmului simplex se bazează pe următoarele 
proprietăţi ale unui astfel de tabel : 


Propoziția I1.5.1. 


Dacă toate diferenţele d, — c, care apar în ultima linie a tabelului 
sînt astfel încât d; — e, > 0, atunci soluția nenegativă de bază care 
apare în ultimă coloană a tabehilui este optimă, adică realizează mavimul 
functiei obiectiv ow + Cy + Cg. 

În cazul nostru, să arătăm că dacă da — €s > 0, da > 0, ds > 0, 
atunci soluția-posibilă (Bı, Ba, 0, 0, 0, Bs) a problemei de progràmare 
liniară 2 este optimă. Pentru problema de programare liniară 1 am 
văzut că se poate obține imediat o soluție optimă, luînd în consi- 
derare primele 3 componente ale vectorului soluție optimă pentru 
problema 2, adică luînd e = Buy Y = fu Z = 0. 

Trebuie să arătăm că pentru oricare altă soluție posibilă a pro- 
blemei 2, de forma (£1, Yy, 2y Ün Vi Zi) CU Vy Yo Zu Ün Yy % > 0 
care. verifică restricţiile (3), valoarea funcţiei obiectiv, adică cu + 
+ Cada + Caz este mai mică sau cel. mult egală cu valoarea funcţiei 
“obiectiv pentru soluția nenegativă de bază (fu; Pa 0, 0, 0, fa), adică 
cu c Pi + Ceba Faptul că vectorul (2, Yy 215 fie Yı 2) verifică (3) 
se scrie vectorial sub forma următoare : 


£m + Vaz T Zig + th + Jat Za = b. (8) 


u; prin egalităţile următoare, deduse din (6): . 


Us = Qly T Aola T aa 


„Va = sh T, Solz T Aaga. 
Înlocuind în (8) pe ax, t, 2 dați de aceste egalitäți, obținem : 
i (2 + zas + Bta + Yit) -H (Ya F Zaza F Erdo F 


+ Fiaos)ta + (2 T 21%a3 + Žids F Yta = b. (9) 
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Însă din (6) mai obţinem că: 


Bia + Bata + Batia = b. (10) 


Deoarece vectorii dı, da; #3 sînt liniar independenţi, rezultă că expri- 
marea vectorului b ca o combinație liniară de acești vectori este 
unică, şi deci comparînd (9) cu (10) deducem egalitățile : 

ZF Aag F Bita + Yrtis = Bo 

Yı T Zoe H Pioa + Yita = Bar i (1) 


Zi T Zita F Brdas F Yita = Pa 


Valoarea funcției obiectiv v, se poate calcula, tinind seama de egali- 
tățile (11) astfel : A $ e egali 


Di = Bala + Balo = (81 + Zas F aaa + Yrs) + 
F (Va -H Ziaa F Braoa + Va a5)0a = 10, + Vu0a + 
-H (agt, H atasa) H Fy (AraC, H degto) + 


H Yalt F a02) = W6, A Vica + ida + Dida + Vida 


tinind seama de definiția numerelor da, da, ds, date de (7). Diferenţa 
dintre valorile funcției obiectiv_pentru soluțiile nenegàtive (8, Ba 
0, 0, 0, Pa) Si (8u Yi Z Fry Yı 21) este egală cu: 


Bila F Vila F Als + Bida + Fads — (C161 + Y1Ca A- 2163) =zu(da — 
— 03) + da + Yad; > 0, deoarece conform ipoteżei d, — 03 > 0, 
d > 0, ds > 0, z > 0, & > 0, J, > 0, deci (B, Ba 0, 0,0, Ba) 


pompis care se citește în ultima coloană æ tabelului este optimă. 
e.c.t.d. 

Această proprietate este adevărată, oricare ar îi vectorii din 
pază şi ne dă un criteriu de recunoaştere a unei: soluţii optime pentru 
problema 2 (deci şi pentru problema 1): cînd toate diferenţele 
d, — 0, care apar în ultima linie a tabelului construit prin metoda 
substituţiei sînt nenegative, soluţia nenegativă de bază obținută în 
acel moment este optimă, adică realizează maximul funeţiei obiectiv 
a problemei de programare liniară 1 sau 2. 
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în caz contrar, cel puţin una dintre diferenţele d; — e; < 0. În 
acest caz să presupunem de exemplu că = 3 pentru tabloul nostru. 


Propoziția I1.5.2. 


Dacă d — c, < 0 şi toți coeficienții de pe coloana lwi az sînt 
mai mici sau egali cu zero, adică œs S 0, aa < 0, ass < 0, atunei 
problema de programare liniară 2 nu are optim finit, adică funcția 
obiectiv poate lua valori oricât de mari (proprietate valabilă şi peniru 
problema 1). 

Pentru a arăta acest lucru să luăm două dintre relaţiile (6): 


b = Bia, + Pata + Bus 
Qs = Oas F agla + Agata 


Înmulţină egalitatea a doua cu A > 0 şi scăzind-o membru cu membru 
din prima obținem: 


b = (Bu — doza) + (Ba — Aaaa) + (Ba — Aos)ua E Ada (12) 


Egalitatea (12) ne arată că vectorul cu componente nene- 
gaitive : (Bı — Aaaa Ba — Adas, M 0, O, Bs — Aass) este o soluţie 
posibilă a problemei de programare liniară 2, deoarece verifică 
ecuaţiile (3) scrise vectorial sub forma (12) pentru zi = Bi — a> 
> 0, Yı = Bz — Ačo > 0, =A >0, Zi=0, J1=0, = Ba — Mag > 0. 

Valoarea funcției obiectiv pentru această soluție posibilă 
este cu( fi. — Agia) + Ca(Ba — Adas) + ACa =t Brt CBot Mea — Crag — 
— 02023) = Cbi F CaBa F Mos — da). Deoarece Cs — d = — (da — 
— c3) > 0 conform ipotezei şi A > 0, se vede că valoarea funcției 
obiectiv poate fi făcută oricît de mare în mulțimea soluțiilor po- 
sibile ale problemei 2, luîndu-l pe A suficient de mare. Deci în acest 
caz funcția obiectiv nu admite un maximum. ¢.c.t.d. 

Cazul complementar este următorul: există un indice 4 
pentru care a > 0, deci mulţimea indicilor I; = ţi | aa > 0) este 
nevidă. În acest caz continuarea algoritmului se bazează pe pro- 
prietatea următoare : 


Propoziția I1.5.3. 
Dacă d, — c, < 0 și mulţimea Is £ Ø, atunci să presupunem 
că min fi = Pa adică minimul expresiilor B după a >0 se reali- 
iE, Qis Xag Ciz 
zează peniru indicele i = 3. (Acest indice care realizează minimul 
poate să nu fie unic.) 


111 


În acest caz substituind în tabel pe u, din prima coloană leu 
a se obține în coloana lui b o soluție, nenegativă de bază pentru 
care valoarea functiei obiectiv crește. 
i Pentru a demonstra această proprietate, să presupunem că, 
tabelul obţinut prin metoda substituţiei cu pivotul «s, încadrat în 
tabelul 2 este următorul : 


& az a3 u LA LA b 
a E 0 0 ia ajs A pi 
de A 9 Lo 25 Zz Ba 
2s, d o 1 Gas 435 ag Ba 
0 0 o di di d vi 
Tabelul 3 


unde elementele tabelului 3 se obțin din elementele tabelului 2 


după regula dreptunghiului, alegînd ca pivot al transformării 
elementul «s, adică: - 


, 43 &: pe 
13 Xa 13 &, 
ai= aa , a= as — 5, 
a, 
33 C2 
, LI : r 
PE = &is Ba i 
xie vi B= B > a24 
aa as 
oa Oa, x a, 
app M Si, cepe = apg — TP 085, 
X33 ag 
2 "ma gt [> A 
a26 > B2 = Ba — 5 Pa, ag = e 
as ss Ua 
r LA , 1 , 
a835 > a3 = —, Ba Ba . 
sa Oas X33 


Primele trei zerouri. din ultima linie s-au obținut ţinind seama că 
diferențele d, — c, pentru vectorii a, aa, a, din bază sînt nule 
El tele di y Par $ 1 2 3 z y) 
lementele di, ds, da şi v se calculează după aceleaşi reguli ca 
pină acum, adică : 
a EDINE y E 
d= 4u F aba H asta — 0, di = aise, + azsa + 
7 i - d 

F asss — 0, de = aas0u F seca + aşs0s — 0, vi = 


= Prea + Bala + Bsta. 
112 


Se poate arăta că şi aceste elemente ale tabelului se pot calcula tot 
cu regula dreptunghiului, de exemplu : 


Sos (da — Ca) Maai Ba (ds — 63) , 


Asg C33 


di =d 


ss Xa 


BA : r . Oa Q: Aagi, 
Într-adevăr, d4 = &146, + 02403 -H 346 = (cu — =a] e, + (ou — a) & + 


Oas EN Zy 
p e= teal H Costa — 7 (Cta Ca H Cog Ca — Ca) = da — — (d — Ca) ete. 
Xss iá as š as 


Să arătăm că Bi Bz 83 sînt numere nenegative: 


1 = 8. — sa fa deci, condiţia Bj > O se transcrie B, > i Ba, care este îndeplinită 
33 a3 


dacă o SS0, deoarece B, > 0, B; > sa > 0. În caz contrar, putem împărţi cu ca 


0, a, 
în ambii membri și obţinem EL > bB, inegalitate care este îndeplinită, deoarece 


13 as 
min Pi Be. La îel se arată că Bz > 0, iar 65 = Pa >20. 
Gis > 0 Zig Oas Ca 


Deoarece b = Bia, + Bia + Bias, rezultă că soluţia (Bi, Bz 
8i, 0, 0, 0) este tot o soluţie nenegaiivă de bază a problemei de pro- 
gramare liniară 2. E 

: Să observăm că deoarece «s, # 0, vectorii din prima; coloană a 
tabelului obținut prin substituție şi anume a; 2, ag, sint de asemenea 
"linia independenţi, ceea ce justifică denumirea de soluție nenegativă 
de bază. ` 

:.Să mai demonstrăm că prin această substituție valoarea funcţiei 
obiectiv crește ` . ran S mn 
Balda: — cs), 


“sa > 
Deoarece d, — Cs < 0, Ba > 0 şi aga >O rezultă că vi > d, adică 
valoarea funcţiei obiectiv se îmbunătăţeşte. Dacă fs >0 rezultă 
vi >v, şi valoarea funcției obiectiv crește strict. e.c.t.d. i 
Dacă pentru tabelul astfel obținut di, de şi de sînt nenegative, 
atunci e = fi, y = Bz, 2 = Ps este o soluţie optimă a problemei 
de programare liniară 1, conform propoziției 1.5.1. În caz contrar 


Conform, regulii dreptunghiului, vi = v; — 


- se alege un coeficient negativ (de exemplu cel mai mic număr negativ) 
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8 — e. 2161 


dintre d, d; şi de şi se continuă aplicarea, algoritmului în modul 
următor : dacă pe o coloană toţi coeficienţii ai; sînt negativi sau 
nuli, problema de programare liniară nu admite o soluţie optimă, 
conform propoziției II.5.2. Dacă, pe coloana celui mai mic număr 
negativ dintre di, di, d; există un coeficient pozitiv, se schimbă 
baza după regula dată de propoziţia II.5.3. După un număr finit 
de pași se va ajunge fie la o soluţie optimă a problemei de progra- 
mare liniară, fie la concluzia că problema nu admite o soluție optimă, 
ea putînd lua valori oricît de mari în mulțimea soluțiilor posibile. 

Am presupus, făcînd această afirmație, că bazele nu se repetă, 
și deoarece numărul lor este finit (în cazul nostru vectorii din fiecare 
combinaţie din ce! mult C} = 20 combinații de cîte trei vectori 
dintre a,, da, as, Ur, tta, U sînt liniar independenți), rezultă că algoritmul 
are un număr finit de pași. 

Dacă toate componentele Bus Bas Ba sînt nenule (deci strict pozi- 
tive), atunci soluţia, de bază (Br Pa O, 0, 0, Ba) se numește nede- 
generată. 

În caz contrar, dacă măcar una, dintre componentele $,, B, sau 
Pa» corespunzătoare vectorilor au, e Şi ta din baza lui R$, este nulă, 
atunci soluţia de bază obţinută se numeşte degenerată. 

Să observăm că valoarea funcţiei obiectiv la o iterajie creşte 
strict dacă soluţia este nedegeneraiă, deoarece în acest caz Ba >0, 
tinind seama, de relaţia dintre vi şi a. 

Dacă valoarea funcției obiectiv nu se modifică în cursul cîtorva, 
iterații succesive, este posibil să revenim la una, din bazele prin 
care am trecut deja și atunci procesul continuă la infinit fără a 
conduce la soluţie. 

Această situație se numeşte ciclare. Am văzut că această situație 
poate apărea numai atunci cînd problema de programare liniară admite 
soluţii de bază degenerate. ; 

Pentru o problemă; nedegenerată, convergența algoritmului sim- 
plex, adică faptul că conduce la soluţie într-un număr finit de iterații, 
este asigurată. 

Pentru o problemă degenerată, adică avînd cel puțin o soluție 
de bază degenerată, este posibilă în principiu ciclarea. 

Subliniem faptul că degenerarea nu implică neapărat ciclare. 
Deşi multe probleme practice sînt; degenerate, nici una nu a cicat 
pină acum şi exemplele de ciclare au fost construite cu dificultate. 
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Pentru înlăturarea ciclării se poate utiliza fie tehnica de perturbare, 
fie metoda lexicografică (vezi de exemplu cartea : M. Dragomirescu, 
M. Malița, Programarea pătratică, Întroducere în programarea 
convexă, editura, Științifică, Bucureşti, 1968). Următoarea, problemă 
ne va furniza o interpretare geometrică, a, noţiunii de degenerare : 


max (s + 4y). 


Fig. 1.5.1. 


Maximul se atinge în vîrful A, adică soluţia optimă este e = 2, 
y = 5 sau, pentru problema adusă la forma standard (prin intro- 
ducerea, variabilelor — ecart 4, J, 3, î în cele patru inegalităţi) : 


æ = 2, y = 5, 5 = 0, y = 16, z= 0, = 0. 


Acest program este degenerat, deoarece are. mai puțin de 4 
componente strict pozitive. De asemenea, soluția de bază corespun- 
zătoare viriului B este : 


æ = 0, y=1,3=0,y=0,z=8,i=2 


și ea este degenerată, > 
Prin ambele vîrfuri, A şi B trec câte trei din dreptele care defi- 
nesc tronsonul convex al soluţiilor posibile, deși două ar fi fost sufi- 
ciente pentru definirea, vîrfurilor respective. 
Foniomenul acesta este general : orice soluţie de bază degenerată 
a unei probleme cu două variabile sub forma, standard corespunde, 
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în problema sub forma canonică, unui punct aflat la intersecția a 
mai mult de două dintre dreptele care definesc tronsonul (incluzind 
pe cele de forma æ = 0 sau y = 0). 

Proprietățile anunțate se pot transcrie pentru un număr oarecare 
de variabile şi de restricții ale problemei de programare liniară, iar 
demonstrațiile se fac pe baza aceloraşi principii şi în cazul general. 

Vom da, în continuare cîteva exemple de aplicare a algorit- 
mului simplex. 


Exemplul 1. 


Să considerăm următoarea problemă de programare liniară : 
Să se determine w, y e R astfel încît : 


sy >o 
2# — 3y < 4 
5w +2% <12 


max (3 + y). 


Reprezentînd grafic în planul sOy mulţimea convexă a soluțiilor 
posibile ale acestei probleme de programare liniară gäsim mulţimea 
punctelor interioare și laturile patrulaterului cu vârfurile în . punctele 


de coordonate (0, 0), (0, 6), (2, o (na) 11.52). 


ea 
Fig. IL.5.2. i 
Calculînd valorile funcției obiectiv 3% + y în cele patru vîrfuri 


ale patrulaterului, găsim că în virtul a această funcţie ia 
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ien aere a 


136 


valoarea maximă, adică max (32 + y) = œ 7,158 în mulţimea, 


19 


soluţiilor posibile ale problemei de programare liniară dată. - 


Deci soluţia, acestei probleme 


de s =#, y San 
. 19 19 


Să rezolvăm aceeaşi problemă 


de programare liniară este dată 


ă cu algoritmul simples, trans- 


formînd mai întîi inegalităţile în egalităţi prin i i 
j galit galităţi prin introducerea variabi- 
lelor de compensare &, y > 0. Obţinem problema : să se determine 


&, yY æ, y > 0 astfel încât: 


28 — 3y + 


5 + 2y - 


i 
“e 


zi] 
Il 


4. 
=12, 


„max (32 + y). 


Pentru această nouă problemă matricea coeficienţilor este urmă- 


toarea :. 
pea ( 
5 i 


vectorul b = (4, 12) iar vectorul e 


—3 10 


2:0 aj . 
= (3, 1, 0, 0). 


- Tabelul cu care plecăm este următorul : 


ay LA LA LA b 
m | 2] -3 1 0 4 
LA 5 2 o 1 12 

—3 —1 0 0 0 


şi care corespunde soluției nenegative de bază : e = 0, y =0,5 = 4, 
y = 12, vectorii a şi ap fiind a, = (2, 5), ax = (— 3, 2) iar, u, = (1,0), 
ta = (0, 1). Elementele de pe ultima linie a tabelului s-au calculat 


astfel : 


— 3 = 2-03 + 5-04 — 6y 


1 3-63 + 2-04 — Ca, 


0 = 1-03 + 0-64 — 63, 


0 = 0:03 + 1-6; — Gas 
0 = 0-3 +0-1+4-0 -+12.0 (valoarea func- 


ției obiectiv pentru soluția nenega; 


tivă de bază (0, 0, 4, 12)). 
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Deoarece în ultima linie există numere negative (în afara 
coloanei lui b), rezultă că soluția (0, 0, 4, 12) nu este optimă. 
Cea mai mică valoare negativă din ultima linie fiind —3, com- 


parăm rapoartele PŁ după au >0, care în cazul nostru sint: 
Qir 

ign, Cel mai mic fiind 2, rezultă că schimbăm baza, înlocuindu-l 

pe m GU Gu, deci alegind ca pivot pe 2, pe care l-am încadrat în 

dreptunghi. Tabelul obținut cu metoda substituției este următorul : 


a, GA LA u, | b 
a, E ee lea 00|. a 
2 2 
19 5 
tlg 0 a T73 í 2 
o d 3 ol] 
2 2 


z 3 iE : P IA si 
Deoarece în ultima linie există un coeficient negativ şi anume — Fi 


(care nu apare în coloana lui b), rezultă că soluția nenegativă de 
bază obţinută şi anume (2, 0, 0, 2) nu este optimă pentru pro- 
blema de programare liniară modificată (prin introducerea variabi- 
lelor de compensare), deci ne uităm în coloana lui a în care se găseşte 
acest coeficient negativ. În această coloană există un coeficient 


ri 19 ; x : ; ; 
pozitiv şi aame deci putem continua aplicarea algoritmului, 


alegînd drept pivot pe = . Se obţine tabelul următor : 


LA LA u u | b 
a 1 Ge e se, E ||| „04 
19 19 19 
5 
a, 0 EIN E 
19 19 19 
o o L n| 008 
19 19 19 


Deoarece în ultima linie diferențele d, — e, (i = 1, 2, 3, 4) sint toate 


nenegative, rezultă că soluția de bază astfel obținută ia 00); 
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este 'soluţie optimă a problemei modificate. Suprimi i 

gi _ emei n e uprimind com- 
pondiitele corespunzătoare lui g, y, obţinem soluția optimă a 
problemei de programare liniară inițială, şi anume : 


s=%, pd. 
19 19 
Tot din ultima linie, în dreptul coloanei lui b găsim valoarea maximă 


a funcţiei obiectiv, şi anume, 
Etemplul 2. 
Să se determine v, y, ze R, astfel incit : 
2, Y, Z > 0, 
æ + 2y — z < 5, 
s — 4y +zsl, 
: max (32 + 2y -+ 42). 
Tntrođucînd variabilele de compensare £, y > 0 obținem problema : 
Li Y, 2 g, yO 
g + 2y —z +4 ġ=5 
s —4y +e +ī=l 
f max (32 + 2y -+ 4z) 
Tabelul cu care plecăm este următorul : 


a A az a LA A b 

i 1 2 -1 i o 5 
uz 1 —4 E o 1 1 
3 2 4 o o 0 


Deoarece există, diferenţe d, — 6 (i = 1, 2, 3, 4, 5) negative, aplică 
algoritmul simplex alegînd ca pivot numărul 1 în ti Satra 
pg ea a a p încadrat în dreptunghi 


| 


a , A a, A u | b 
LA 2 —2 0 1 1 6 
% l -4 1 0 i la 

1 —18 o o 4 | 4 


f = i loana - lui 
Deoarece în acest tabel dz;— 02= — 18 < 0, iar pe coloa 
æ toţi coeficienţii « (i = 1, 2) sînt negativi, rezultă că problima 
nu are optim finit, funcţia obiectiv 32 + 2y A de putind lua valori 
oricît de mari în mulţimea soluţiilor posibile, 


Exemplul 3. 
"Să se determine æ, y > 0 astfel înċit : 
~e +y<si 
+ 2y s4 
2. 2... 
max (Tæ -+ 4y) 
ic î i iilor ibile ale 
ntind grafic în planul 20y mulţimea, soluţiilor posi e 
Enea Ma de programare liniară obținem mulțimea convene 
din figura TI.5.3 cu virfurile de coordonate (0, 0), (2,0), (2,1), (F) 


i . Funcția- obiectiv îşi atinge maximul în vîrful (2,1). > - 
a pc a aplica algoritmul simplex introducem. variabilele de 
compensare î, y, 3 > 0 şi transformăm inegalităţile în egalităţi : 


—g+y+g=i 
s+ 2y+y=4 
a+2=3, 


funcția obiectiv rămînînd neschimbată. 


PTSTA ZEIT ET AN 


Tabelul de la care plecăm este următorul : 


| a az u, EA LEA b 

Ra —1 1 1 0 0 1 
BEA ES 2 6 1 o 4 
ua +] o 0 0 1 2 
-7 —4 0 0 0 0 


Pe ultima, linie alegem cea mai mică diferenţă d, — ce; negativă, şi 
anume — 7. Pe coloana lui —7 alegem coeficienţii pozitivi şi com- 

i 2 4 3 E a D ; 
parăm rapoartele Z și T> alegîndu-l pe cel mai mic Ii Pivotul 


transformării va, fi numărul 1 încadrat în dreptunghi. Se obţine 
tabelul : | . 


a a ZA LA Us. b 
2 0 1 1 0 1 | 3 
23 o SEI] o 1 —1 2 
a, 1 g o -0 1 2 
o -4 o o 7 14 


Pe ultima linie a tabelului mai există o singură diferență negativă 
(—4). În coloana, corespunzătoare lui —4 selectăm coeficienţii pozi- 
tivi şi alegem pivotul comparând rapoartele 3. şi = . Minimul fiind 
Zi Tuăim pe 2 ca pivot şi obținem tabelul : 


t üa i 


LA us | o 

a. o o 1 EEN 2 2 
qe 2 2 

a o 1 o: ES EREN 1 
2 2 

a, 1 .0 0 o 1 2 

o o o 2 5 las 


Soluţia, nenegativă de bază (2, 1, 2, 0, 0) este soluţie optimă pentru 
problema de programare liniară cu variabile de compensație, deoarece 
toate numerele scrise în ultima linie pe primele 5 coloane sînt 
nenegaiive. j i 


Scrierea componentelor s-a făcut tinind seama de ordinea 
vectorilor din baza scrisă în stînga tabelului, componentele corespun- 
zătoare vectorilor ce nu fac parte din bază fiind nule. Pentru a 
obţine soluţia optimă a problemei inițiale neglijăm variabilele de 
compensație, obţinînd æ = 2, y = 1. Valoarea optimă a funcției 
obiectiv se citeşte din tabel din coloana lui b, ea fiind egală cu 18. 


Observaţie. 


Se poate arăta, că soluţiile nenegative de bază ale unei pro- 
bleme de programare liniară cu restricţiile de forma egalităţilor 
corespund unor vîrfuri ale mulțimii convexe a soluţiilor posibile. 
Algoritmul simplex, lerînd numai cu soluții nenegative de bază, 
constă deci din trecerea de la un virf la alt virf al mulţimii convexe 
a soluţiilor posibile ale problemei cu mărirea valorii funcţiei obiectiv, 
pină cînd se ajunge la un optim, dacă, acesta există. 

Dacă problema de programare liniară are forma standard, 
adică forma : 


ant + Gaz) + Ms? = by 
Goa F Gaz) + Aaa = ba, 


aat + aga + asa = ba, (12) 


S, Y, 2 > 0, 
max (ou + ey + 0a), 


vom indica cum poate fi obținută o soluție nenegativă de bază: 


iniţială, cu care poate fi inițiată aplicarea, algoritmului simplex. 
Metoda, care va fi expusă se numeşte metoda celor două faze. 
Putem presupune că în egalităţile problemei (12) avem b > 0, 
ba > O şi ba > 0, înmulţind în ambii membri ai restricţiilor cu b; < 0 
prin — 1. ; 
Adăugăm apoi în fiecare ecuație cîte o variabilă artificială a, 
Yu 2, Şi vom considera problema : 


ant + Aaa + Gas? F 2 = biy 


ant F aY F Az + Yı = bzs (13) 


ayt + Gas - Gas? +% = bas 
2, Y, 220; ty Yva >0, 


max (— da — Yı — 2). 
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i 
| 
| 
| 
| 
] 


Este clar că această nouă problemă de programare liniară ; 
soluție nenegativă de bază, şi anume : A E 0, sira pe iată 
„Putem deci aplica algoritmul simplex pentru această problemă, 
plecînd cu vectorii unitari ty, tz, da, în prima coloană a tabelului” 
Deoarece 2, Yy în > 0, rezultă æ, + Yı t 2 > 0 şi deci max (—a să 
— a — 2) = — min (m, + Y + 2) < 0. Sînt posibile două cazuri : 
1) Dacă max (—®, — Yı — 2) <0 în condiţiile (13), atunci 
problema sub formă standard (12) nu admite soluții nenegative 
Într-adevăr, dacă (12) ar admite o soluție nenegativă (30 ap 20), 
atunci (13) ar admite o soluţie nenegativă de forma (2%, yo 2, 0 
0, 0) şi deci max (—®, — Yı — 2) = 0 (pentru m, =y, = 2 2 ). 
i 2) Dacă max (—a — Yı — 21) = 0 în condiţiile (13), atunci am 
găsit o soluţie nenegaiivă de bază (a%, y°, 20, Œi, 98, 24) pentru pro- 
blema (13) cu —a0 — y} — 2 = 0. Deoarece a, 4, 2? sint nene- 
gative, rezultă cu necesitate că ad =y} = 20 =—0 şi se vede că 
soluţia (2, Y°’, 20, 0, 0, 0) verifică ccuațiile problemei (12) și este 
nenegativä, deci am obținut o soluție nenegativă de bază pentru 
problema (12) sub forma standard. 

„Deoarece toți vectorii 2, ws, ttg, corespunzători variabilelor 
artificiale au fost eliminați din bază, putem incepe faza a doua, care 
constă din rezolvarea propriu-zisă a problemei (12) cu ajutorul algo- 
Saul, simplex, pornind de la soluţia nenegativă de bază obţinută 
i aza Intu. 


, 


Exemplu. 


Să considerăm următoarea problemë de programare liniară în 
formulare standard : 


2g -+ y — z = 3, 
æ + 3y —-ti=4, 
2, Y, z, t> 0, 
max (— æ — 6y). 
În faza întîi vom transforma problema în următoarea : 
2g + y~z 


g + 3y —t +y = 4, 
z, Y 2% i Fis Y > 6, 


ra = 3, 


max (— 2 — Yı). 


Pentru această nouă problemă matricea coeficienţilor este urmă- 
toarea : : 


2 223 
af? 1 1 
1 3 0 


0 1 0 
—1 0 1 


Vectorul b = (3, 4), iar vectorul e = (0, 0, 0, 0, —1, —1). 
Tabelul de la care plecăm este următorul : 


«a LA LA LA LA EA 
DA 2 a —1 o 1 0 3 
IA 1 |3| o —1 0 1 4 
—3 —4 1 1 0 0 —7 


Deoarece în ultima linie există numere negative (în afara coloanei 
lvi b), rezultă că soluţia (0, 0, 0, 0, 3, 4) nu este optimă pentru 
problema auxiliară. ; e AAA N 

Cea mai uucå valoare negativă din ultima linie fiind —4, com- 
parăm rapoartele după a > 0, găsind min fil = pi 

Cig 

Îl vom alege ca pivot pe 3, introducînd pe a, în bază în locul 

lui ta Obținem tabelul : 


| a da a3 i i Ha b 

mo o| 5/3 | 0 —1 1/3 1 —1B8 5/3 

EA 1/3 1 0 —1/3 0 1/3 41/3 
—5/3 0 1 —1/3 0 43 | —5/3 


În ultima linie a tabelului există acum 2 coeficienți negativi. 


- Cel mai mic este — 5/3, iar din coloana sa îl alegem ca pivot pe 


5/3, cu regula cunoscută. Obţinem tabelul final al fazei întii : 


| a de 23 fi U u | b 
a 1 0 —315 15 3]5 —1/5 1 
a | 0 1 2/5 —215 —15 215 1 
| 0 o o. o 1 1 0 


rai 


Deoarece în ultima linie toate numerele sint; nenegative, f întăi 
s-a încheiat, obţinînd soluţia nenegativă de bază a problemei inițiale : 
(1, 1, 0, 0). Plecînd de la această soluţie, trecem la faza a doua 
aplicînd algoritmul simplex pentru tabelul : i 


a, aa i; CA b 
a „d 0 —315 1/5 | 1 
a, o 1 [ars] -235 i 
0 o —3/5 41/5 | -7 


Se observă că elementele —3/5 şi 11/5 ale ultimei linii au trebuit 
să fie recalculate, restul coeficienţilor fiind identici cu cei din tabelul 
obținut la sfârșitul fazei întîi. 

k Alegînd ca pivot pe 1/5 încadrat într-un dreptunghi, obținem 
tabelul : Ta 


a; da i EA b 
a, 1 3 0 i 4. 
a 0 5 1 —2 5 

0 3- 0 $ —4 


Deoarece toate numerele din ultima linie (exceptind coloana lui 
b) sînt nenegative, rezultă că şi faza a doua este încheiată, obținind 
soluția optimă : œ = 4, y = 0, z = 5, t = 0. 
Problema inițială era echivalentă cu problema : 
22 +y > 3, 
g + 3y > 4 
v > 0, y >00, 
min (x + 6y). 
Pentru aceasta, mulțimea soluțiilor posibile este reprezentată in 


figura II.5.4. 


Fig. 11.5-4 


Se observă că minimul funcției obiectiv se atinge în punctul 
Cit, 0), iar aplicarea algoritmului simplex a constat din deplasarea 
din B(1, 1) în C(4, 0), unde funcția obiectiv îşi atinge valoarea 
minimă în mulțimea soluţiilor posibile ale problemei de programare 
liniară, 


Il. 6. Probleme 


1. Să se transcrie matematice următoarea problemă de progra- 
mate liniară : 

Într-o fabrică funcționează două procese de producţie, notate 
Pi; Pa care produc substanțele $,, Sa, Sa, astfel încît în unitatea, 
de timp (ora) P, produce a, unităţi din 8, @a unităţi din Sa, dus 
unităţi din S, iar P, produce a, unităţi din S, @z unităţi din 
Sz, Gaza unităţi din Sz. Prin plan fabrica trebuie să producă cel puţin 
b, unităţi din 8, d, unități din Sa, bg unităţi din Sa. Se ştie că 
procesul P, consumă în unitatea de timp ce, lei (materii prime şi 
manoperă), iar P, consumă în unitatea de timp 0, lei. Procesele de 
producţie fiind liniare (adică procesul P, (i = 1, 2) funcţionind t 
ore consumă te; lei), să se determine un plan de producţie optim, 
care să asigure îndeplinirea: planului cu cheltuieli de producţie 
minime. 

2, Să se rezolve următoarea: problemă de programare lini iară : 


2>0,y>0, 
25 +y>5 
a+4y >6 
min (2 + 6y), 


reprezentind grafic mulțimea convexă a soluţiilor posibile în planul 


z0y. 
R: (6, 0). 
3. Să se rezolve următoarele probleme de programare liniară, 


folosind reprezentarea. grafică a mulţimii soluţiilor posibile : 


a) 4,5 2 — 3,2 y < 4,2 
3e +7,25 y < 2,5, 
2, Y > 0, 
max (—34 + 2,6 y). 
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b) —2g + 8y < 5, e) 2g + 3y > 4, 
3s + 4y < 10, gts, 
s—yx<si, i 4 y > 0, 

æ, y z0; i max (52 — y). 
max (22 + 4y). R: (7, 0) 
f 35 ) 
R:|— a 
8 32 i 
d) &—y <2, e) ix +y s3 
æ, y > 0, p+ysl 
max (22 + Ty). æ+ 3y >i 
R: Nu există maximum. £, Y > 0, 


max {tæ — 39). 


4. Se ştie că fiecare animal dintr-o fermă are nevoie de o canti- 
tate minimă de principii nutritive în unitatea de timp (pe zi), care 
depinde de specie, vîrstă, scopul urmărit în alimentaţie ete. 

Principiile nutritive se află în diferite proporţii în produsele ce 
compun rația furajeră, 

Folosind datele din tabelul de mai jos, să se determine canti- 
tatea æ din alimentul A, şi cantitatea y din alimentul As, expri- 
mate în kg, ce trebuie să intre în compoziţia raţiei furajere a unui 
anima] astfel încît costul ei să fie minim, folosind metoda grafică. 


| Conţinutul în principii 


Denumirea princi- nutritive al alimentelor Cantităţile 


viilor nutritive te) minime pre- 
A dz 
P, 0,1 0 | 0,4 
P, o 0,1 0,6 
| P. P 0,1 0,2 2 
0,2 0,1 1,7 
Cosa (lei/kg) 2,4 0,8 ra 
E 3 
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R: Modelul matematic al problemei este următorul : 
`æ > 0, y >O, 

0,1 > 0,4, 

0,1y > 0,6, 

012 + 0,24 > 2, 

0,2% + 01y > 1,7, 

min (2,4 æ + 0,8 y). 
Folosind metoda grafică se găseşte æ = 4, y = 9. 

5. Să se rezolve următoarea problemă de transport: Din A 

trebuie transportate 2500 t produse și din B trebuie transportate 


1400 t la depozitele D, și Da. D, poate primi 2800 t, D, poate primi 
1500 t, iar matricea costurilor transportului este următoarea. : 


D, Da 
A 80 leii 40 leijt 
B 50 lei/t 35 lei/t 


Să se determine planul de transport optim. 


R: Probiema se transcrie matematic sub forma : 


e > 0,y>0,2>0,i>0, 


æ + y = 2500, 
a + t = 1400, 
æ+ 2 < 2800, 


y +t s 1500, 
min (80% -+ 40y + 502 + 35t). 


Mulțimea soluţiilor posibile din planul Oz este desenată în figura 


1.6.1. Se găseşte soluţia optimă : 
= 1000 £, y = 1500 î, z= 1400 t t= 0 tone. 


(1000,1800) 


Fig. II. 6-1 
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6. Să se rezolve următoarele probleme de programare liniară 
cu algoritmul simplex (dacă: este cazul, problemele de minimizare 
se vor transforma în probleme de maximizare, transformind de 
asemenea toate inegalitățile în < ): 


a) s— ysl, b) æ —y +z > — 1, 
2y -+z < 4, 2y — 2 <S 4, 
my 2>0, g, Y, 220, 
max (e — y + 32). . min (— 2g +y — eg). 
R: (1, 0, 4) 5 R: Nu există minimum 
€) gA+2y —2s4, d) s+y—z<1, 
—2g + y +32 <2, 2% + y + 32 <3, 
s + 3y -+z <4, 2—y+ast, 
t, Y, Z> 0, &, Y, z > 0, 
max (2 + 2y + 32). max (æ +y +2). 
l R: (2, 0, 2) n:( TE +) 
2 2 
e) z+ysi, f) 2% +y < 3, 
y+asl, 2A+2s3, 
gz+azsl, 2z -+ æ < 3, 
g, Y, Z > 0, %, Y, 2 > 0, 
max (æ +y +2) max (æ +y +2) 
n(4 Ay +) R: (4,11) 
2 2 2 : 


3a+y+2e+ăisil, 
2, Y, Z, t> 0, 


max (æ + 2y + 32 + t). 
i R: (0,1; 0, 0) 
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7. Scriind variabila æ = z — t eu z, t > 0 să se rezolve urmă- 
toarea problemă de programare liniară cu algoritmul simplex, apoi 
să se verifice rezultatul folosind reprezentarea grafică : 


y > 0, 
22 +y < 4$, 
—e +y <4, 


max (—22 + y). 
R: (— 4, 0) 


8. Folosind metoda de la problema anterioară, să se rezolve 


următoarea problemă : 
By y > 0, 


a+y-thsi 
—e -Hy +2 <4, 
max (x + 23y — z2). 
R: Nu există maximum 
9. Să se rezolve problema : 
æ> 2, y> 3 z0, 
z+y+as6, 
2% + 3y — z < 14, 
max (® + 2y + 32) 


făcînd o schimbare de variabile (2! = æ — 2, y’ = y — 3). 
R: (2, 3, 1) 


10. Să se rezolve problema : 
g >z 6 yz —3 zzl, 
æ +y +2 sð, 
< 5, 
22 + 3y — 2 < 2, 


max (2% + 2). 
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11. Să se rezolve problema : 


2, Y, 2>0, 


max (æ +y — z2}. 


da 
Nils —.o0 
7 Te 


12. Introducind variabile-ecart şi utilizind metoda celor 
două faze, să se rezolve următoarea problemă de programare liniară : 


e—y sl 
zi+y>l, 
æ — ysl, 
a, y >o, 
min (22 + 3y) 
și să se arate că soluția ei optimă este degenerată. 
R: (1, 0) 


13. Să se arate că intersecția a două mulțimi convexe este o 
mulţime convexă. Să se dea un exemplu din care să rezulte că 
această proprietate nu mai este adevărată pentru reuniunea a două 
mulţimi convexe. 
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